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THERMODYNAMICS  OF  RAMJET 
FLOW  PROCESSES 

I.  INTROOUCTION 


In  principle,  the  ramjet  is  an  air-breat aiiiK  jei.  propulsion 
heat  engine  which  depends  for  its  continued  operation  on  ram  compres¬ 
sion  of  the  entering  air,  since  It  contains  no  moving  parts  to  effect 
this  process.  Figure  1  shows  a  schematic  diagram  of  a  conventional 
ramjet.  However,  the  constructional  simplicity  of  the  fixed  geometry 


ramjet  somewhat  belles  the  control  complexity  necessary  to  ensure  suc¬ 
cessful  operation  over  a  wide  range  of  flight  speeds  and  altitudes,  as 
will  be  emphasized  by  the  various  chapters  of  this  volume.  The  ramjet 
Is  primarily  a  supersonic  propulsion  device,  since  at  subsonic  speeds 
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lis  efficiency  of  uperation  is  very  low.  In  addition,  because  of  its 
Innorcnt  dependence  cn  r?.m  compression,  the  ramjet  cannot  produce 
usable  thrust  to  <  r..a(>le  it  to  start  mo^'ing  under  its  own  power,  and 
it  suffers  therefore  from  the  disadvantage  of  needing  to  be  boosted 
to  flight  speed  by  other  means,  such  as  by  a  rocket  or  turbojet.  Thus 
it  cannot  lam  the  sole  power  plant  of  an  aircraft.  At  present  ram¬ 
jets  operate  at  Mach  numbor.s  up  to  about  3,  but  rt  is  confidently  ex- 
poctcii  that  future  developments  of  hypersonic  units  may  raise  this 
figure  to  as  high  as  7  or  above.  Such  a  hypersonic  ramjet  is  as  yet 
in  the  destgn  stage,  but  the  Inaications  are  that  it  will  look  very 
different  from  designs  now  current.  The  principles  of  propulsion, 
however,  will  remain  much  the  same. 

The  performance  of  a  ramjet  engine  is  analyzed  by  first  making 
a  force  balance  at  the  exit  and  inlet  in  order  to  find  the  net  thrust. 
This  is  generally  accomplished  by  determining  the  flow  by  the  use  of 
one-dimensional  gas  dynamic  theory  at  both  stations  and  then  integrat¬ 
ing  the  forces  acting  on  corresponding  areas.  The  net  thrust  is  gen¬ 
erally  converted  into  a  dlmeiisioules.o  coefficient  (the  thrust  coeffi¬ 
cient)  by  dividing  it  by  the  product  of  an  engine  reference  area  and  the 
freestream  dynamic  pressure.  I’erformance  is  also  a  function  of  com¬ 
pression,  expansion,  and  temperature  ratios.  The  use  of  a  fixed  exit 
design  for  an  engine  eliminates  the  variable  of  expansion  ratio  and 
makes  possible  a  convenient  plot  e'  performance  as  a  function  of  com¬ 
pression  ratio  as  derived  from  the  flight  Mach  number  and  the  engine 
temperature  function  referred  to  later  (this  last  quantity  is  not  a 
dimensionless  ratio).  Analytical  estimation  of  the  performance  of  a 
ramjet  engine  and  the  determination  of  its  optimum  design  cannot 
readily  be  accomplished.  The  flow  of  air  through  an  engine  is  set  by 
the  parameters  of  heut  addition  to  the  combustion  chamber  and  the  exit 
nozzle  throat  area,  and  the  equations  relating  exi+  and  inlet  condi¬ 
tions  are  usually  of  a  complex  Implicit  nature.  Accordingly,  calcu¬ 
lations  initially  set  out  to  relate  conditions  at  the  combustion 
chamber  inlet  to  those  at  the  exit  nozzle  throat,  and  diffuser  r 
formance  characteristics  are  then  used  to  determine  the  diffuser  in¬ 
let  conditions  on  the  requirement  that  the  flow  at  the  diffuser  exit 
match  that  at  the  combustion  chamber  inlet. 
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A  I  '.'riaiii  ar.icunt  of  Introductory  thermodynamic  theory  is  pre- 
■^entcd  ii  ‘'Cf',  ions  2  r.nd  3  of  this  TG  since  the  application  of 
this  raalerial  is  discusseo  in  Section  3  onwards.  Fuel  and  flow 
controls  are  discussed  in  more  detail  in  TG  370-5  and  -6  as  fuel  sys¬ 
tem  design  considerations  are  .an  outgrowth  of  propulsion  requlrenients ; 
such  devices  when  properly  designed  are  not  subject  to  'feedback'  from 
tile  engine  proper.  Other  aspects  of  component  performance  will  be 
discussed  in  more  detail  in  TG  370-4  on  Diffusers,  in  TG  370-7  through 
-12  on  Combustion  Systems  and  in  TG  370-13  on  Exit  Nozzles.  In  gen¬ 
eral,  this  TG  will  attempt  to  relate  theory  to  practice,  and  engine 
performance  to  that  of  its  components.  ('unfortunately,  however,  much 
otherwise  suitable  material  has  had  to  be  withheld  as  It  has  not  yet 
beer,  publishrii  in  the  open  literatui  e . )  Consistent  energy  units  are 
u.sed  tin  ougliuut ,  su  thu ..  Uic  vyi  the  conversion  factor  J  (the  mechani¬ 
cal  equivalent  of  boat)  is  not  needed.  The  symbols  used  are  in  accord 
with  the  A.S.A.  List  of  Preferred  Symbols  as  far  as  possible. 
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A  i))  ’i  !  .uu'.  I,  i!ii  1  mil  A.'vo lopmont  of  useful  gas  flow  rela- 
!  :  i'-.  1 lure  i;  'senttti  t'oth  for  f  ami  1  lar  i  7.ai  ion  ot  ttie  novicr, 

>-  .1  re*  -e-'e.' ,  Such  flow  is  first  considered  from  the 

-  i;i.  nsio'.a  1  .i-peci,  i.e..  with  constant  properties  across  a  plane 
■;j3.  iiuic’.i  lar  to  tta  ciircction  of  flow.  (More  strictly,  one -dimer. - 
I'hil  i'n  .v  tw  '.iuuight  of  as  t.iat  in  which  the  component  of  velo- 

I  .t;  .\t  1  iptit  to  the  direction  of  the  average  velocity  at  any 

•  os.,  section  may  !k'  .issumod  to  be  neg’igible.)  Vorw  rigorous  treat- 
..  a  l  ol  one-dineiisior.al  flow  dues  not  appear  justified  as  real  flows 
la  vi'r  trulv  oiii-d i nt'iis ion.o  1  .  but  are  treated  as  such  onlv  to  make 
an.il..-;s  traetable.  Kclevant  flow  phenomena  that  do  not  fall  into 
i.t  o;;.  iiiit  its  luna  1  eatef'orv  are  then  nonsidered  in  due  course. 

Ill  go.icT,' 1  ihe  treatment  is  concerned  with  ideal  nonviscous 
ilows,  Kince  these  cctnstitute  a  large  part  ot  the  subject  matter  of 
:)viid  mechanics,  and  does  not  refer  to  the  nonuniformities  introduced 
liy ,  lor  caamplc,  the  growth  ol  the  more  slowly  moving  boundary  layers 
; loh  form  adjacent  to  a  solid  surface,  and  in  which  the  fluid  /elo- 
city  is  progressively  reduced  to  zero.  (It  is  usually  possible  to 
.ioal  separately  with  such  viscous  effects,  and  then  fit  them  into  the 
1,1.1  in  flow  field.)  Certain  boundary  layer  effects  will,  however,  be 
reierred  to  l.ater  In  connection  with  the  subject  of  flow  through  noz¬ 
zles  . 


2.1  Basic  Thermodynamic  Relationships 

Any  nna'ysis  of  the  working  fluid  requires  consideration  of  the 
lul lowing  fundamental  laws: 


1.  Conservation  of  Matter:  p  A  V 


constant . 


iTHMis  iwii.i  i  Gas 


2.  Newton's  Laws  ol'  Motion; 

(a)  Matter  at  rest  remains  at  rest  or  in  uniform  motion  un¬ 
less  acted  upon  hy  external  forces. 

(b)  The  rate  of  change  of  momentum  of  a  body  is  proportional 
to  the  applied  force  and  takes  place  in  the  direction  of 
that  force.  This  definition  is  applicable  to  a  fluid. 

(c)  To  every  action  there  is  an  equal  and  opposite  reaction. 
(This  implies  that  forces  never  occur  singly,  but  that 
every  force  on  a  body  is  exerted  by  apnlying  an 

and  opposite  force  to  another  body.) 

3.  The  First  Law  of  Thermodynamics:  All  energy  is  conserved; 
when  energy  in  one  form  disappears,  an  equal  amount  of  energy 
reappears  in  .mother  form. 

4.  The  Second  Law  of  Thermodynamics:  Heat  will  not  of  its  own 
accord  flow  from  a  cold  to  a  hot  body.  Heat  will  therefore 
flow  only  where  there  is  a  favorable  temperature  differen¬ 
tial.  This  law  governs  the  behavior  of  all  heat  engines 
since  these  depend  for  their  operation  on  the  flow  of  heat 
into  and  oat  of  a  working  fluid. 

5.  (a)  The  perfect  gas  law 

F  p  R  T,  and 

(b)  The  adiabatic  (and  Isentropic)  law 

Y 

p/p  =  constant  . 

Both  the  perfect  gas  law  and  the  adiabatic  law  are  idealized 
representations  of  the  observed  behavior  of  gases.  The  departure  from 
the  perfect  gas  law  is  often  taken  care  of  by  expressing 

P  =7 

riTT 

where  Z  is  a  compressibility  factor  which  may  differ  from  unity.  R 
is  used  here  in  preference  to  the  universal  gas  constant  3!  (to  which 
it  is  related  by  the  “"nation 


R 


weight  ) 
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si!ico  pi'oso'ii  applications  deal  essentially  with  the  same  fluid.  For 
the  iMso  of  nonisentroplc  flow  of  gases,  the  adiabatic  law  is  con- 
sideri'd  to  be  a  particular  case  of  the  polytropic  law 


P/p"  “  constant, 


where  n  may  have  various  values.  When  n<y  ,  heat  flow  takes  place 
into  the  fluid  during  an  expansion  process,  and  when  n>  y  ,  heat  flow 
takes  place  away  from  it.  If  n  -  1,  the  equation  represents  an  iso¬ 
thermal  How  case. 

The  main  problem  to  be  considered  in  ramjet  flow  analysis  is 
that  of  tracing  the  one-dlmenslonal  gas  stream  In  terms  of  iiS  state 
parameters  through  various  flow  transformations.  If  gas  initially  at 
a  state  A  is  transformed  into  a  new  state  B,  its  final  state  will 
be  a  function  both  of  its  Initial  state  and  of  the  process  through 
which  it  Will  have  passed.  Some  processes  are  easy  to  follow  because 
one  can  relate  the  initial  and  final  stream  states  by  equating  the 
difference  of  state  to  the  prescribed  change;  if  there  is  no  change  of 
state  this  difference  is  zero.  Analysis  therefore  depends  upon  the 
determination  of  the  stream  states  for  such  processes  at  various 
points . 


The  obvious  intrinsic  properties  of  a  gas  are  mass,  volume,  and 
pressure  or  temperature  (since  these  latter  may  be  related  to  each 
other  by  the  perfect  gas  law).  If  the  gas  is  in  motion,  it  possesses 
kinetic  energy;  it  also  possesses  potential  energy  by  virtue  of  its 
pressure  and  volu.me,  but  its  gravitational  energy  of  position  is  nor¬ 
mally  negligible.  Four  Independent  variables  a^c  thus  required  to 
desc.'ibe  completely  the  state  of  the  working  fluid— the  physical  con¬ 
stants  uf  the  gas  apart.  The  two  familiar  thermodynamic  variables, 
pressure  and  temperature,  may  be  replaced  by  point  stream  properties 
which  are  uniquely  related  to  them.  One  may  thus  establish  new  ther¬ 
modynamic  state  properties  by  manipulation  of  the  equations  stating 
the  basic  laws.  For  example,  the  First  Law  of  Thermodynamics  states 
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that  in  general  a  quantity  Q  of  heat  added  to  unit  Mass  of  a  gas  in¬ 
creases  the  Internal  thermal  energy  u  and  perforu  mechanical  work 
W.  This  Is  expressed  In  differential  form  as 


dQ  -  dll  dW , 


Mechanical  work  may  be  expressed  in  of  force  and  distance  as 


dW 

-  (Forci' ) 

But 

Fn  (» 

-PA. 

dW 

-  P  A  d. 

whence 

dQ 

-  du  -  P 

For  Ideal  gases  the  Internal  energy  u  is  a  lunction  o;  temperature 
only,  yielding  the  relationship 

dQ  -  c^  dT  ♦  P  dv.  (2) 


As  It  stands,  this  differential  equation  is  not  directly  Ictegrable, 
since  It  has  not  yet  been  established  what  portion  of  the  added  heat 
Increases  the  Internal  energy  of  the  body  and  what  portion  Is  con¬ 
verted  to  mechanical  work.  It  may  be  shown  from  the  Second  latw,  how¬ 
ever,  that  the  right  hand  side  of  the  equation  may  be  made  an  'exact' 
differential  by  dividing  by  the  temperature,  l.e.. 


(3) 


This  equation  defines  a  quantity  s,  called  entropy,  which  Is  thus 
shown  to  be  a  function  of  the  thermodynamic  state  of  the  system,  and 
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tt'iiich  In  fact  determines  i^.ethcr  or  not  the  flow  oe  reKarded  as 
reversible  Ols  -  0).  A  reversible  adiabatic  process  Is  therefore 
iscr.li'upic .  Conversely,  an  irreversible  adiabatic  process  is  non- 
Isent  roplo . 

If  the  h'^at  addition  be  specified  to  occur  at  constant  pressure, 
itj.  •;  i )  Is  directly  Integrable  since  P  will  not  then  be  a  function  of 

V  : 


P  constant 


H  “  n  +  Pv. 


(4) 


The  quantity  H  so  defined  Is  known  as  enthalpy,  which  thus  repre¬ 
sents  the  heat  absorbed  by  a  body  at  constant  pressure.  By  defini¬ 
tion,  Cp  ”  (^Jp  ‘  *  **‘®  Independent  functions 

involving  T,  P,  and  v,  they  also  serve  to  define  the  thersodynaaiic 
state  of  a  systee.  Tlic  above  relations  are  essentially  theraos-catlc 
ones,  since  they  presuppose  that  the  systee  has  attained  tberaodynaalc 
equlllbriua.  At  very  high  flow  velocities,  this  is  not  necessarily 
the  case. 

Many  tberwodynaalo  calculations  can  be  conslitorably  slaplified 
if  It  Is  assumed  that  the  heat  capacities  c^  and  c^  reaain  constant  as 
the  temperature  varies,  which  would  be  so  if  all  the  internal  energy 
of  the  gas  existed  In  the  fom  of  translational  energy  of  the  BOle- 
cules.  In  combustion  processes,  however,  temperatures  are  high  enough 
for  the  Internal  energy  to  exist  In  the  rotational  and  vibrational 

temperature  is  therefore  of 
that  y  should  be  given  by 


forms  The  variation  of  y  (-  c  /c„)  with 

'  p  V 

some  importance,  and  klnet'^c  theory  shows 


y  -  1  . 
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ro 


whtTP  X  is  th'j  numbor  <»f  dctjreos  of  freedom  ov  the  gas  molecule  ■>  3 
for  monaiomu'  gases,  5  for  dlitloinlc  gases,  an<l  6  for  polyatomic  gases. 
The  corresponding  theoretical  values  of  y  are  thus  1.67.  1.40  and 
1.33.  >’or  the  case  of  air,  the  value  of  1.  M)  is  very  nearly  correct 
at  lo*  temperatures  and  pressures,  the  value  for  0^  being  slightly 
lower  and  that  for  N2  slightly  greater.  (The  value  for  combustion 
products  differs  somewhat.)  The  effect  of  pressure  on  the  value  of 
y,  particularly  within  the  relatively  low  pressure  range  at  which  the 
ram.jct  engine  operates,  i.e.,  less  than  about  200  pda,  is  sma'".  Fig¬ 
ure  2  shows  the  variation  of  y  with  temperature  and  pressure  (Ref. 

1)!  The  value  of  1.40  is  a  convenient  one.  in  that  Its  use  results  in 


TEMS€»*TURE  (‘Rl 


Fiq.  2  Vofiotion  of  )  (e^  c, )  with  Temptrofurt  ond  Prttturt 


the  exponents  common  In  various  Mach  number  functions  referred  to  be¬ 
low  having  whole  number  values,  for  example,  | 

I 


5, 


and 


y  +  1 
2Ty  -  1) 


3. 


y  ■  1.40  has  therefore  also  been  Incorporated  as  a  basic  parameter  Into 
the  International  Standard  Atmosphere.  .'.Ithough  for  the  sake  of  simpli¬ 
fying  analyses  many  relatively  simple  empirical  equations  have  been 
proposed  to  express  c^  and  as  functions  of  temperature,  these  equa¬ 
tions  can  at  best  yield  only  average  values  over  certain  ranges  of  tem¬ 
perature.  It  Is  therefore  preferable  to  use  actual  data,  particularly 
as  the  burden  of  handling  this  material  can  be  alleviated  to  a  large 


References  may  be  found  on  page 


173. 
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e.'itent  by  the  use  of  computers  with  large  storage  capacities.  (In 
practice,  a  computer  would  store  all  the  equilibrium  constants,  vibra¬ 
tional  properties,  and  moments  of  Inertia,  and  then  utilize  a  program 
to  calculate  the  thermodynamic  constants  aS  needed.) 

2.2  The  Speed  of  Sound 

This  IS  a  parameter  of  major  significance  in  defining  the  nature 
of  high  speed  fluid  flow.  A  derivation  is  accordingly  presented  here, 
the  starting  points  being  Newton‘8  Second  Law  and  the  Law  of  Conserva¬ 
tion  of  Matter.  Consider  first  a  mo«inc  coordinate  system  which  has  the 
same  velocity  as  a  weak  shock  or  sound  wave  (a  small  pressure  discontin¬ 
uity  moving  at  the  speed  of  sound)  so  that  the  wave  becomes  stationary 
with  respect  to  the  moving  coordinates.  If  the  wave  is  taken  as  being 
propagated  from  right  to  left  with  respect  to  the  gaa  at  rest,  then  in 
tr.c  Biu.  ir.g  r-ccrdir.ate  system  there  Is  a  flow  of  gaa  from  left  to  rlt.'l 
through  the  wave.  Since  this  wave  is  considered  one-dimensional,  l.e., 
planar,  it  extends  theoretically  to  Infinity  In  a  direction  normal  to 
that  of  propagation.  Under  these  conditiona  there  ia  thus  no  change 
along  the  wave  normal  to  the  flow  direction;  furthermore,  a  stream 
tube  entering  from  the  left  must  retain  the  same  cross  section  .'.nd 
the  mass  rate  of  flow  through  it  remains  constant,  even  though  velocity, 
density,  temperature,  and  pressure  say  change  as  the  tube  passes  through 
the  wave  front,  as  in  Fig.  3. 
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i3y  Newton’s  Thira  Law.  the  i.-rci'S  acting  front  the  left  on  the 
sount!  wave  must  balance  those  acting  from  the  right.  The  reactions  on 
the  two  sides  are  therefore  respectively  the  sums  of  the  pressure 
forees  .iiiu  the  laics  of  change  ot  momentum; 

"l  ^1  ^  Sr  (""l  =  A.  e  (n,^  V2) 

tor  steady  flow. 

Vj  and  arc  const.oit  with  time  at  fixed  locations  and  their  deriva¬ 
tives  are  thus  zero. 

dm .  um„ 

■  •  ^1  *  ''1  ”3?  “  ^2^2  ^  '^2  “dT  • 

The  mass  rates  of  flow,  dm/dt ,  are  given  by  the  Equation  of  Continuity, 
or  of  Conservation  of  Matter; 

dm^  dmg 

“3t  ■“  ^1  Pi  '^l  -31  ‘  ^2  P2  ''2  • 

^1  ^1  ^l^S  Pi  '^l^  '  ^2  **2  ''^2^^2  P2 

2 

or  A(P  r  p  V  )  =  constant. 

This  latter  will  be  referred  to  as  the  momentum  equation.  The  quantity 
2 

A(P  +  p  V  )  is  constant  along  any  stream  tube,  in  the  absence  of  wall 
forces,  and  can  be  differentiated  to  determine  how  P,  p,  and  V  vary 
as  the  flow  passes  through  the  sound  wave.  For  a  constant  area  strean, 
tube , 

A(dP  +  2  V  p  dV  +  dp)  “  0,  (5) 

but 

A  p  V  =  constant. 
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Differentiating  the  last  expression, 


A(p  dV  +  V  dp)  -  0 


Inserting  Eq.  (6;  Into  Eq  (5), 


AldP  -  2  dp  +  dp)  -  0  , 


and  solving  for  V  , 


,2  dP 


If  the  flow  velocity  is  high  enough  for  there  to  be  negligible  heat 
conduction  from  the  wave  front,  the  flow  will  be  adiabatic,  for  which 

p 

—  -  -  constant , 


dP  -  y  ^  dp  , 


So  that 


which  for  ideal  gases  gives 


V  -  Ky  R  T 

-  a,  the  speed  of  sound.  (8) 

y  is  here  the  'instantaneous'  value  corresponding  to  the  temperature  T. 
Also , 

2  dP 
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Thih  equation  ciofino^  the  speed  at  which  a  amplitude  disturbance 

will  propaK-ate  itself  through  the  medium.  This  sonic  velocity  will  be 
toui'd  to  he  a  eonvonient  quantity  for  simplifying  analyses  since  the 
iiach  number  M  is  defined  In  terms  thereof.  It  may  again  be  emphasised 
that  the  equation  .a  =  yV~R  T  is  only  valid  for  a  condition  of  thermo¬ 
dynamic  cnuilibrium,  which  may  not  always  be  a  realistic  assumption 
■vhen  combustion  gases  are  being  considered.  If  a  chemical  lag  occurs, 
a  a.s  given  above  must  be  regarded  as  a  Icwei-  liiuiiing  value.  A  higher 
limiting  value,  determined  on  the  assumption  that  t  lo 

is  'frozen' ,  i.e.,  its  chemical  composition  is  unable  to  react  to 
ch.Tiiging  state  parameters.  This  distinction  may  be  imporlaiil  if  the 
flow  velocity  is  changing  rapidly,  as  in  a  nozzle. 

Sound  waves  arc  very  weak  types  of  shock  waves,  the  theory  of 
which  is  dealt  witli  in  greater  detail  below.  The  speed  of  sound  is 
thus  the  rate  at  which  only  the  smallest  pressure  disturbances  travel. 
True  shock  waves,  in  which  the  pressure  ratio  across  the  discontinuity 
may  be  much  greater  than  unity,  may  on  the  other  hand  travel  at  speeds 
of  several  times  that  of  sound.  In  this  category  fall  explosion  waves 
and  certain  types  of  detonation  waves  in  combustion  chambers. 

In  view  of  the  definition  of  M  in  terms  of  a,  it  is  of  inter¬ 
est  to  consider  how  the  value  of  a  is  affected  by  a  variation  in  tem¬ 
perature,  as  both  y  and  R  are  temperature-dependent.  In  the  case 
of  R,  where  R  *  vt>p  -  Cy) .  this  difference  is  almost  constant  with  T 
up  to  temperatures  of  several  thousand  degrees  Fahrenheit,  f  is  there¬ 
fore  the  parameter  most  affected  by  the  variation  of  T,  and  in  the 
atmosphere,  since  the  temperature  varies  in  a  complicated  way  with 
altitude,  so  also  does  the  speed  of  sound,  the  latter  ranging  in  value 
from  somewhat  less  than  lOCO  fps  in  thp  stratosphere  to  over  1400  fps 

^  at  400,000  feet. 

1 

2.3  The  Mach  Number 

^  The  nature  of  the  wave  system  representing  small  disturbances 

traveling  with  the  speed  of  sound  is  dependent  on  the  absolute  (free- 
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sti-'  an'.)  velocity  oi  'he  source  of  the  disturbances  .in  relation  to  the 
:-:peod  at  which  the  disiurhances  travel,  or,  more  succinctly,  on  its 
Mac)i  nu.Tibei-  M.  a  parameter  which  has  already  been  referred  to.  In 
Fi^r.  1(a),  let  sutccssive  positions  of  the  souj'ce  of  a  disturbance 


(a)  (b) 

Fig.  4  Mach  Patterns  Set  Up  by  Small  Pressure  Disturbances  Moving 
at  Subsonic  and  Supersonic  Velocities 


'traveling  at  velocity  V  <  a)  be  A,  B,  and  C.  When  the  point  C  has 
been  reached  after  a  time  2  t,  the  front  of  the  wave  emitted  when  the 
source  was  at  A  will  have  traveled  a  distance  2  a  t;  that  emitted 
from  position  B  will  only  have  traveled  a  distance  at,  i.e.,  the 
source  can  never  overtake  the  wave  representing  the  'advance  warning'. 
In  Fig.  4(h),  however,  when  V  >  a,  it  is  seen  that  there  exists  a  coni¬ 
cal  envelop**  to  the  wave  fronts  (which  in  three-dimensional  flow  are 
spherical).  Supersonic  flow  is  thus  characterized  by  the  existence  of 
a  'zone  of  silence'  external  to  the  envelope,  in  which  disturbances 
representing  pressure  signals  cannot  be  received.  If,  for  example, 
the  source  of  disturbances  is  stationary  and  the  fluid  is  moving  in 
the  opposite  direction,  i.e.,  from  right  to  left,  supersonic  flow  will 
continue  in  the  zone  of  silence  undisturbed  by  the  source  whose  pres¬ 
ence  is  only  felt  downstream  of  the  'Mach  cone'.  The  envelope  shown 
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IS  known  as  a  \la<'h  wave  (in  either  two-  or  three-dimensional  flow), 
and  it  is  eiin::'  tinit 


fM  -  —  —  . 

V  t  M 

II  is  here  the  Mach  an^le  If  the  source  of  a  disturbance  moves  at 
the  speed  of  sound,  fj.  v  2 . 

The  wave  conf itturat ioii  shown  in  Fig.  4(b)  will  remain  similar 
even  when  the  pressure  disturbance  is  laige,  except  that  the  velocity 
of  propagation  will  now  be  greater  than  a,  and  the  half-cone  angle 
will  thus  be  larger  than  sin~^  1/M.  Such  a  wave  is  known  as  an  oblique 
shock  wave.  Since  it  is  traveling  faster,  a  stronger  shock  wave  will 
thus  overtake  a  weaker  one,  and  by  a  similar  process,  a  series  of  weak 
waves  will  coalesce  to  form  a  stronger  one,  as  will  be  the  case  if  the 
boundary  has  a  continuous  concave  curvature. 

It  is  convenioiit  to  express  many  of  the  common  thermodynamic 
flow  relations  in  the  form  of  functions  of  the  Mach  number,  and  various 
relations  involving  M  are  therefore  derived  in  a  subsequent  section, 
and  plotted  in  the  Appendix. 


When 

M 

< 

1, 

the 

flow  is  said 

to  be  subsonic 

when 

M 

1 ,  it 

is 

transonic , 

and 

when 

M 

> 

1, 

it 

is 

said  to  be 

supersonic , 

hypersonic  flow  referring  to  the  case  when  M  >  5  approximately.  This 

latter  definition  is  merely  one  of  convenience,  since  the  ranges  in 
which  supersonic  and  hypersonic  flow  relations  are  valid  will  of  neces¬ 
sity  overlap  (Ref.  2).  The  Mach  number  may  also  be  regarded  as  an  in¬ 
dication  of  the  ratio  of  the  ordered  kinetic  energy  of  flow  of  a  fluid 
to  the  random  kinetic  energy  oi  molecular  motio.i. 
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2.4  Isftitropic  Compressible  Flow 

The  simplest  type  of  flow  to  consider  is  steady  flow  without 
t)aiisfer  of  heat  or  work  across  the  boundaries  of  a  duct.  The  i.OuJi- 
licjns  at  any  two  Stations  (1)  and  (2)  may  be  related  by  noting  the 
equality  of  the  energy  states;  the  energy  conserved  is  partly  thermal 
(internal)  energy,  partly  nressuve  (potential)  energy,  and  partly 
liinetic  energy.  Writing  expressions  for  these  terms,  one  has  for 
unit  raa.ss ,  neglecting  gravitational  effects, 

2  2 

c  c  T+—  ”  constant.  (9) 

p  2  V  p  2 


For  stagnation  conditions  (reduction  of  kinetic  energy  to  zero). 


c  T, 

V  t  p, 


c  T  4  —  1“  TT— 
V  P  2 


From  the  isentropic  law.  using  the  perfect  gas  relationship. 


y-i 


Let  this  ratio  be  K. 


Si. V 


M  V>  F.  T  1  ir  I  hi’  (  Ifivt  ol  a  compross ible  medium, 


•)  R  T 


“  ^  r,  1 

1 )  y  R  T  ^  “  K 


IK  -  11  . 


r  -  1 

i  . 


From  F.q .  ill),  the  mas.s  flow  per  unit  area  is  given  by 


(12) 


(13) 


Th'i  mass  f  lov?  cannot,  howeve*- ,  be  made  to  increase  Indefinitely  with 
iti-reasing  P^.  for  the  flow  velocity  cannot  be  acceler.ated  beyond  the 
so, lie  value  by  a.  simple  contracting  section. 


For  a  maximum,  therefore, 


I-  2  y+l-] 


(14) 


'or  air  (y  >•  1.40)  this  ratio  has  the  value  of  0,328,  and  for  combus¬ 
tion  bases  it  is  about  0.548.  The  critical  velocity  corresponding  to 
this  pressure  ratio  may  be  found  by  substitution  of  F.q.  (14)  in  Eq . 
(11)  : 


V 

c 


2y 

y  +  1 


R  T., 


y  R  T 


a . 


(15' 
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rims  a  i;}a.\j:num  (ijischar(;o  occurs  when  reaches  a  value  (1/0.528)P  for 
ail  ,  tlie  Mach  iiutiber  beint:  then  equal  to  unity.  This  is  known  as  the 
'(hokiiv;'  condition,  i.e..  an  Increase  of  upstream  pressure  will  not 
further  Increase  the  velocity  bevend  the  local  speed  of  sound  [al- 
ilK'ugh  the  mass  flow  can  still  increase,  because  of  the  presence  of 
llm  density  tei-.m  in  Kq .  (l.'l)'  .  The  location  where  this  occurs  is 
known  as  a  'throat,'  from  which  small  pressure  disturbances  cannot 
tr.".v.  1  upstream,  as  previously  demonstrated  for  the  case  of  the  zone 
Ilf  si  lenco  . 

Summarizing;,  tie  pressure,  density,  and  temperature  ratios  may 
!)  "  r  i  I  t  on  as 


>■ 


(16) 


1 


T 


(18) 


It  is  interestiiiK  to  note  that  the  total  or  stagnation  temperature  T^^ 
is  independent  of  the  method  by  which  the  gas  is  brought  to  rest,  pro¬ 
vided  no  external  work  is  done  and  no  heat  is  transferred  across  the 
boundary.  The  relationship  in  Eq .  (IG)  is  therefore  true  whether  or 
not  the  flow  is  isentropic. 


These  relations  are  of  great  value  in  determining  the  nature 
of  tiie  flow  in  a  variable  area  duct.  From  Bernouilli 's  equation  in 
dif  f  nrenti.al  form, 


V 


dV  + 


P 


0 


Onc-Dinu‘iis  lor.a  1  G.th 


Suhst'.  ImI  nu  from  .i‘ 


cIP  (i/>  . 


a" 

P 


-  V  av 


-  1)2  dV  since 


Wheii,  r).iioforc.  M  <  1.  the  velocity  will  Increase  faster  than  the 
density  can  decrease,  so  that  tlie  area  must  decrease  for  the  flow  to 
continue  lu  fill  the  passage.  When  M  >  1,  the  density  will  decrease 
faster  than  the  velocity  will  increase,  so  that  the  area  must  increase. 
Nozzles  to  produce  supersonic  flow  are  constructed  so  as  to  make  use 
of  this  principle,  and  are  thus  convergent-divergent.  To  find  the 
area  relationship  in  terns  of  M  for  isentropic  flow,  one  has 


D  A  V 


•  *  • 

p  A  V  , 


where  the  asterisk  refers  to  the  throat  condition  (at  which  U  -  1),  so 


.  *  „*  *  * 

A  ^  p  V  _  p  a 

A*  p  V  p  aM 


Furthermore 


p  Pt  P 


K  -If 


y -  1  =  i+y 


Y  P  p 


=  — r  .  Since 
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Hcncc 


A 


1_ 

It 


(1&) 


The  corresponding  (Isentropic)  pressure  may  be  dortved  fro*  Eqs .  (!G) 
and  (19) : 


(20) 


The  relationship  between  area,  pressure  and  Mach  number  may  likewise 
be  obtained  from  Bernoullli’s  equation  and  the  continuity  relation  in 
differential  form: 


V  dV 


0 


and 


dp  dV 

-i  *  -? 


dA 


0  . 


Now 


since 


dA 

T 


dp  dV 
p  -  -V 


~  ^  ’  since 

P  p'T 


dV 


dP 
■  pV 


(21) 


P  V 


2 


pM^y  RT  -  PyM^. 


Since  matter  is  conserved  in  a  flow  process, 


m  -  ^1  ^1  ^1  “  ^2  ^2  ^2  ’ 
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V  F.  T'  '^.wi  V  -  Ma. 


ri;  "i  ‘1  -S 


rr:  ‘’2  ''2  ^2 


S i lu’c 


S  Pi  «1 


Y 

ITT^ 


V  y  R  T, 


^2  P? 


“2  Vr“ 


ITTJ 


(22) 


Equation  (22)  is  In  an  inconvenient  form  since  the  temperatures 
arc  'static'  values,  as  would  be  determined  by  an  'observer'  moving 
i  tl'.  the  stream.  In  attempting  to  measure  temperature,  a  sensing  de¬ 
vice  would  normally  bring  the  flow  to  rest  and  would  thus  sense  the 
'total,'  or  stagnation,  value.  Combining  Eqs .  (18)  and  (22),  there¬ 
fore,  to  get  over  this  difficulty. 


Ai  Pj 

1 

^  „  -^2  ^2  ‘*2  1 

L  R  J 

■  R  J 

(23) 


Suffixes  for  are  unnecessary  here  since  remains  constant  in 
adiabatic  flow. 


Thus  far,  flow  quantities  have  been  considered  on  a  mass  basis. 
For  flow  on  a  weight  basis. 


(24) 
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1  tii'fi.iCu  tn  t  ie  sytr'.';.!  *  and  lo  read  'w  eirclc.  ^  Weight  flow  may 
t  nus  be  expressed  I'v 


! 


1 

I 


I 


■ 

I 


w 

Replacing  the  static  pressure 
.  ('I'.e  has 


APw 

vTi 


(25) 


P  In  Etj,  (25)  by  the  total  pressure 


-  [<P.  P  )w].  (26) 

\t;  I  *  -1 


Equation  (26)  may  be  used  to  relate  constant  weight  flow  between  any 

two  points  (1)  and  (2)  for  two  values  P  and  P  of  total  pr“=sure 

^1  ^2 

which  may  be  dissimilar.  Thus, 


w 


1 


^  ^t. 


(27) 


For  constancy  of  entropy  and  enthalpy,  the  total  pressure  as  well  as 
the  total  temperature  will  be  the  same  at  S+ition  2  as  at  Station  1. 

Thu*5 , 


(P/'P.  )w 

J  1 


(2«) 


The  function  j(P/P.^)wj  may  also  be  expressed  in  terms  of  the  Mach  num¬ 
ber  : 


^he  symbol  *  here  replaces  the  symbol  S,  erroneously  listed  as  a  mass 
flow  function  in  Ref.  22. 


2 


I 


i.)i 1  G:’.s 


i 

i 


i 

i 

i 


«2  ry‘ 

Mj  LK^J 


(29) 


It  A2  bo  replaoocj  by  A.  as  at  f  h*»  throa^  of  a  sui>ersonic  nozzie. 


(30) 


It  follows  that 


[(p  p^)S]m 


01* 


[(i  Pt)w] 


.Kid  11  is  also  obvious  that 


(31) 


(A/A  )j  (A/A  )2 


Several  of  the  above  relations  are  plotted  (for  two  repi-^sentative 
values  of  y)  in  terras  of  U  in  the  Appendix.  These  area  relations 
may  be  used  to  find  the  areas  corresponding  to  different  values  of  M . 
both  sub-  and  supersonic.  When  flow  is  assumed  to  be  isentropic,  no 
change  of  total  pressure  through  Ihe  system  is  considered.  [Static 
pressures  may  be  computed  at  any  point  by  using  Eq.  (16).  1  Since 


TP7P777 


TP7pp7 


(32) 


so  that  (P/P^ )  at  both  Mj  and  Mg  may  be  found  from  t.ables  or  curves 
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siil.s  *  i  • t ;uto  V-n  .  (J2).  For  4>xamplt‘,  if  a  flow  acce  Ici^itos 
i :  ;in  M  D .  1  to  M  ^  1  (i  tho  value  ol  (P  )  chanttus  from  0.993  to 

<1.528  Ik.iir, 


[o'^3  J  ‘'j 


M-0 . 1 


(P  is  .1  cenllnuously  decrcasiiit’  and  s inKle-valueii  runction  oi  »iach 
num'.icr  . 

When  the  flow  is  adiabatic  but  nonisentropic ,  Eqs .  (16)  and 
(dl)  lio  Icr.jjer  apply,  but  Eq .  (27)  may  be  used.  The  amount  of  energy 
degradation  (entropy  increase)  can  be  specified  by  means  of  the  total 
pressure  ratio  (P  ,'P  ).  Using  Eq.  (27),  an  exit  Mach  number  and 

^  t  . 

4.  1 

area  may  be  related  to  an  Initial  Mach  number  and  area  as  follows: 


A 


2 


(P/P 


^2  ^1 


(33) 


This  relationship  may  be  solved  for  either  the  unknown  area  or  Mach 
number  function  as  required,  and  It  Is  useful  for  the  calculation  of 
ramjet  inlet  capture  area  ratios  for  imperfect  diffusers — where  the 
total  pressure  recovery  ratio,  exit  Mach  number,  and  freestre-am  Mach 
num'/er  are  known.  In  this  case,  the  equation  is  set  up  for  Stations 
0  and  2  (see  Fig.  1),  and  is  expressed  in  terms  of  area  ratios: 


A 


7  ‘’to 


[_(i*. 


.'p. 


)Sl 


Solving  for  (A^/A^), 


Aq  ^  /  t2'\/A2\  [(P/P^)w]2 

^  \^/\V  [(P/P^)81(, 


(34) 


All  tho  terms  on  ttie  rii^ht  l.ind  jiclo  are  known.  Alternatively,  if 
••ill  tile  terms  oxeejii  (!>  ,P  )  are  measurable  quantities,  this  rela- 

lion  may  bo  used  to  de..e’  ine  the  total  pressure  recovery. 


2.5  The  Normal  Shock 

One-dimensional  gas  dynamics  is  adequate  to  describe  case  of 
noiiisentropic  flow  represented  by  a  discontinuous  rise  of  ^  essure , 
density,  temperature,  and  entropy  known  as  the  normal  shock  wave.  This 
vise  of  pressure  is  of  particular  importance  in  the  design  of  a  ramjet 
diffuser  which  is  designed  to  produce  compression  of  the  entering  air 
without  the  use  of  moving  parts. 

Nov'  from  Euler’s  equation 


V  dV 


P 


0, 


and  the  isentropic  equation 


P  ■=  p  ^  k  , 

which  can  be  differentiated  to 

y-1 

dP  =  K  T  p'  dp  , 


there  results 


^bjk> 


y  dp 


the  stagnation  value.  This  is  the  energy  equation  in  compressible  flow 
and  may  be  applied  to  either  side  of  a  shock  wave. 


2ti 
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Co.isiJoi  now  ti.“  flow  in  a  duct  through  a  normal  bhock  front,  using 
subscripts  1  and  2  for  the  irunediate  upstream  and  downstream  loca- 
t ions  (Fig.  o )  . 


Then  for  the  two  sides  of  the  shock--which  may  be  considered  infinitely 
thin — using  the  continuity  and  momentum  relations, 


1  2 


From  the  first  two  parts  of  Eq .  (35), 
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2  2  /k'Vk*'' 

«1  “2 


•  •> 
K  M“ 


‘1  K 


Sdlviir;  for  tjy  subat  itutiiig  for  K.j^ ,  K2 ,  and  K 


1  +  __2_ 


2YMj^  -  Y  +  1 


^^2  _  a* 
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The  variation  of  (V2/Vj^)  and  (Mg/Mj^)  with  is  shown  in  Fig.  6.  Pres¬ 
sure,  density,  and  temperature  relations  across  a  normal  shock  wave  may 
be  found  as  follows: 


P2  '^1 

Pi  ^ 


*1  K7  • 


From  the  momentum  relationship 


^2  -  Pi  -  Pi  Vi'-PaVg"  . 


and  since 
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\>  ' 
1 
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P,  ’  ■  P, 
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2 
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y  a 


y  -  1 
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-  n 

+  X  ) 


1)  . 


Hence 


P 
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2 
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2. 

4 
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(Mf 
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K*  +  y(Mj  -  1) 

- s - 

K 


2yM^  -  (y  -  1) 

~~ 


(40) 
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'.iiU'  tli;vi  ,iS  M.— O'  (5\.,,  Pj>  •».  O'.  ,  l>ut  " — “‘I-  (which  has  the 

■■  '  ’  IS"  !  r  ;  i  '}  i  , 


These  relations  arc  shown  In  Figs.  7  ami  S.  Figxu'e  5  shows  a  ccspari- 
ooa  between  the  curve  relating  pressure  and  density  ratios  for  Isen- 
tropic  flow,  and  Xhal  tor  flow  across  a  normal  shocK  wave 


I  2  3  4  R 

M, 

Fig.  7  Kormol  Shock  Functions  ( Pj/P, )  ond  (pj/pi ),  y  =  1.4 


At  very  high  Mach  numbers,  the  static  temperature  rise  of  the 
^as  passing  tiii'uugh  a  shock  wave  may  give  rise  to  considerable  ioniza¬ 
tion,  which  may  be  felt  ahead  of  the  wave  if  the  ionized  particles  are 
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iM.  ; : n;  t  *!•.  i  t'J.itu  :'  ^  aii'i  M,,  .o'ro^;s  a  normal  shock  wave 
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(42) 


It  is  ulBo  a  simpio  natter  to  lieternine  the  ratio  of  the  stagnation 
pressure  downstream  of  a  shock  to  the  upstream  static  pressure.  These 
(luiiititios  represent  readings  conveniently  available  from  the  use  of 
an  impact  pro’oe  or  i>i»ot  tube  inserted  into  a  supersonic  stream.  As 
■saown  tn  lig.  10,  the  detached  bow  wave  has  the  form  of  a  normal  shock 
i mi:i>  d iaie  1  y  ahead  of  the  probe. 
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FiQ.  10  Pi-of  Tube  ir.  'lupersonic  Stream 
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-  K*  1|2 

2  K* 

2(2y  -  y  +  1) 

2y  -  y  +  1 
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(43) 


This  relation  is  known  as  the  Rayleigh  pitot  tube  formula  and  Is  plot¬ 
ted  in  Fig.  A-12  in  the  Appendix. 


It  may  be  noted  that  a  velocity  transition  from  subsonic  to 
supersonic  through  a  normal  shock  wave — though  permitted  by  the  above 
consideration--is  not  possible,  as  this  would  involve  a  contradiction 
of  the  Second  Law  of  Thermodynamics.  The  total  pressure  relatlcnchip 
shows  that  a  loss  occurs  during  deceleration  (as  actually  observed), 
whereas  flow  in  the  reverse  direction  would  imply  an  increase  in  total 
pressure  and  thereby  a  reduction  of  entropy. 


2.6  Steady  Dlabatlc  Flow 

The  foregoing  adiabatic  analysis  obviously  cannot  refer  to  the 
case  in  which  flow  takes  place  with  a  change  in  enthalpy,  such  as 
occurs  through  the  combustion  chamber  of  a  ra;ajet,  or  in  a  duct  where 
fj iction  effects  are  not  inappreciable.  This  may  conveniently  be 
handled  by  one  of  the  methods  relating  to  dlabatlc  processes,  i.e., 
those  in  which  heat  Is  added  to  the  system  either  from  within  or  across 
the  walls,  as  developed  hy  Shapiro  et  al^.  (Ref.  3)  and  ctheis. 

From  the  momentum  equation,  taking  subscripts  (1)  and  tk)  be¬ 
fore  and  after  the  point  of  heat  addition. 
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a:\d  from  the  perfect  gas  law»  P  -  p  R  T,  and  ¥  -  ^  ij  “  -  p" 
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FroB  the  perfect  gas  law, 
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Likewise 


■^1  •‘l  M^(l  ^  y  Kj 

^^*^2  m|(i  +  y  Kg 


Pj  Vg  M^(l  ^  y  M^) 

^2  ''l  Mj(l  +  r  Mg 


(47) 


(48) 


1 

Since  p  -  K  (where  the  den8lt>  assumes  its  stagnation  value 

if  the  flow  is  brought  adiabatlcally  to  rest), 


and 


Mgd  +  y  Mj) 

- - - 

Mj(i  +  y  Mg) 


(49) 


^tj  Mgd  +  y  Mj) 

'’tg  Mjd  +  y  M^)  \’^2/ 


(50) 


Since  the  stagnation  temperature  changes  with  heat  addition,  a 
more  convenient  reference  condition  is  one  associated  with  a  particular 
value  of  the  mass  flow.  This  is  the  condition  at  which  'thermal  chok- 

4i  4i  4i 

iag’  (M  “  1)  occurs.  There  is  thus  but  one  value  of  T  ,  P  and  p  for 
any  particular  flow. 


The  corresuondiiu'  pressure,  temperature  and  density  equations 
are  thus  found  by  substitution  of  M  -  1  into  Eqs .  (44)  to  (50): 
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From  Kq.  (48), 
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(52) 


(53) 


(54) 
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Theae  functions  are  plotted  In  Figs.  A-16  and  A-17  in  the  Appendix 
for  two  values  of  y  and  they  are  also  ta))ulated  ir.  various  texts 
listed  In  the  Bibliography  (Cashel  and  Jennlnga,  Zucrow,  and  Shapiro) 
and  in  Ref.  22.  For  Itiltlally  subsonic  conditions,  Addition  of  *'eat 
to  the  flow  causes  P,  p^,  p  and  P^  to  fall,  and  V,  and  T  to  rise. 
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.  ti.  t  A,  ,1,',  11!'  uhi  !i  iu:U  w  1 1  lulr;iwii  from  the  fluid.  This 

J'i  l  1  ;p!'  I.ae  !'■  <a  .ip;'l!cci  i  ii  the  ‘  aero  I  lu’ I'mopressof  '  ileve  i  i)y 

hiiar.ii-ii  Ohl.  n.  lu'i'ip  i  ti.jri-ti'd  into  a  sub.soiiii'  stream 

llowi!,;  la  a  iiu>  t  ; i  a  di  r  to  iiuroasc  the  stagnalioti  pressure.  It 

* 

1  !  !  lx  Cl,  a.it  le-'d  that  an  .aiumalv  occurs  in  the  'ui'i'es  tor  (T.  T  ) 

and  (p  p)  111  that  ihe.se  ratios  reacii  a  raar  mum  at  a  value  somewhat 
It'ss  than  unity.  It  is  easily  .seen  from  d  i  f  f  oreti  t  ia  t  i  hfi  Eq  .  (53)  that 
this  m.i.'.imum  oi-  urs  at  M  1  %  >  ,  i.e..  at  M  (1.845  for  >'  ”  7/5,  and 
at  M  i)..S82  lor  }  ‘.'.7.  For  1'%/  <  M  s  1,  therefore,  where 

1.029  •  (T  T*)  •  1.  lurthor  addition  of  heat  to  the  fluid  brings 
about  an  increase  in  kinetic  energy  to  compensate  for  the  drop  in  the 
value  t.'l  (T  T  ).  The  behavior  of  (p  p)  is  similar. 

Since  with  flow  initially  subsonic,  it  is  clear  from  the  ther¬ 
mal  choking  effect  that  the  addition  of  moie  heat  to  the  fluid  cannot 
then  result  in  a  further  increase  of  velocity,  the  assumed  mass  flow 
condition  at  the  entrance  to  the  heat  transfer  region  becomes  incom¬ 
patible  with  that  at  the  thermal  throat,  and  the  mass  flow  must  then 
decre.'>se.  This  reduction  in  mass  flow  is  analogous  to  the  pressure 
Choking  effect  resulting  from  a  progressive  reduction  in  the  size  of 
a  geomecric  throat.  The  situation  in  initially  supersonic  flow  may  be 
illustrated  by  a  comparison  with  the  case  of  flow  through  a  duct  with 
multiple  (geometric)  throats,  as  shown  in  Fig.  11,  which  is  an  exten- 
siv'Jii  ot  the  nozzle  condition  discussed  later  with  reference  to  Fig.  18. 
The  adjustment  to  sonic  conditions  at  the  thermal  throat  is  accom- 
jjlished  by  the  flow  first  becoming  subsonic  through  a  normal  shock;  it 
may  then  accelerate  again  if  further  heating  occurs  (as  by  friction) 
up  to  the  sonic  point,  finally  decelerating  through  a  normal  shock. 

At  least  in  theory,  shock-free  deceleration  of  supersonic  flow  to  Mach 
1  by  heat  addition  is  possible,  but  if  friction  effects  are  present, 
it  is  clear  that  the  length  of  .a  duct  is  an  important  variable,  since 
the  effect  of  friction  is  to  bring  about  a  progressive  addition  of 
heat  to  the  fluid,  with  consequent  change  of  Mad.  numher- 
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Fig.  11  Analogy  af  Diobotic  Flow  of  Air  with  Ihot  in  o  Duct  Having 
Multiple  ((..eomelric)  Throats,  Showing  Unstable  Nature  of 
Decelerating  Flow 


It  may  be  noted  that  the  creation  of  a  therm.:  '  throat  is  not 
normally  aimed  at  in  the  design  cf  ramjet  combustion  chambers,  since 
the  function  of  accelerating  the  flow  to  sonic  velocity  is  reserved 
to  the  propulsion  nozzle . 

For  no  output  of  shaft  work,  the  energy  e.;uation  (per  unit  mass) 
is  as  follows: 


!  I,  t  ■  era  t  1  iiL  . 
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1  2 


9 

-  V  .  “i 


or,  in  terms  ol  stagnation  enthalpy. 


«1 


2 


T.  )  , 

‘1 


'■■•here 


When  the  continuity  anh  momentum  equations  apply,  then  for  a  constant 
area  duct  the  .r.aas  velocity  G  is  given  by 

G  =  pv  , 

and  therefore 

2 

P  +  3  /p  -  constant.  (57) 

Likewise,  for  constant  stagnation  temperature,  the  energy  equa¬ 
tion  may  be  written 


yZ 

H  +  — j  =  constant.  (58) 

The  effects  are  shown  graphically  on  a  diagram  plotted  to 
ordinates  of  temperature  (or  enthalpy)  and  entropy,  such  as  Fig.  12. 

The  lines  on  the  diagram  which  represent  Eqs .  (57)  and  (58)  are  known 
an  Rayleigh  and  Fanno  lines  respectively.  The  points  A  and  B  of 
maximum  entropy  for  the  two  curves  correspond  to  the  choking  condition, 
which  can  thus  be  approached  along  either  the  subsonic  (upper)  or 
supersoni-c  (lower)  bro.iches.  A  process  ot  thermal  choking  brought 
about  Dy  heat  addition  to  initially  subsonic  flow  is  shown  by  follow¬ 
ing  along  the  line  DB;  heat  addition  to  supersonic  flow  with  decelera¬ 
tion  through  a  shock  wave  (with  corresponding  Increase  of  entropy)  ard 


One  Dimensional 


sut'iit  «  U-T’;uii>n  til  I  111’  (’liol’.irK  condition  is  illustrated  by  the 
curvi,'  CRB.  riic  piunl  oi  n.ixunuTi.  cnlhalpy  on  the  Rayleigh  line  corre- 
ti  the  value  .M  1  ■,  I 


h 


Fig.  12  Fonno  ond  Rayleigh  Lines 


Some  practical  outcomes  of  these  considerations  as  related  to 
supersonic  combnsT.i on  arc  ef erred  to  later. 


2.7  Thrust  and  Drag 

Further  use  will  now  be  made  of  the  momentum  equation  as  derived 
initially  in  the  section  on  the  speed  of  sound.  In  the  form  given, 

2 

P  •*  nV  =  constant, 

■a  constant  area  duct  was  assumed.  The  momentum  char.»4es  that  the  inter¬ 
nal  flow  through  a  ramjet  undergoes  due  to  both  area  variation  and  heat 


aaviition  ,m  r.iaiio  us**  *>1  in  Iht*  (leiitirat icii  of  thrust.  Consider  first 
1.....  1  .•  t  or  11*  1  thiMsi  IS  sot  op  uy  tJuid  flow  in  a  duct  of  variable  area 

(lio.  13). 


On  any  cross-sectional  area  A,  the  force  acting  is  composed  of 
the  static  pressure  over  the  area  p’us  the  rate  of  transport  of  momen¬ 
tum  ,  i  .e  .  , 


P  A  +  V 


dm 

at  ’ 


or,  using  the  mass  flow  relation  rh  =  p  V  A, 

If  -  A(P  +  pV^) 

=  AP(l  +  y  M^)  for  a  perfect  gas. 


A  function  f  is  also  defined  by 

f  =  P(1  +  y  M^)  . 


(59) 


This  function  is  shown  plotted  in  the  Appendix.  If  i.s  here  known  as 

the  stream  thrust  (or  -ninulse  function).  It  is  clear  that  the  stream 
thrust  remains  constant  ir  a  duct  of  constant  area  (such  as,  for  ex¬ 
ample,  certain  types  of  ramjet  combustion  chamb-.cs)  and  thus  corre¬ 
sponds  to  tlie  Rayleigh  line  previously  referred  to.  In  a  variable 
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.1  I  .1  iiUC  I 


hv 


;•  .  ai- !  1  lurre  oxi-rtei*  h;  the  fluid  Oi'.  the  duct  is 

■  inlei'iia!  •■nd  nctweon  Stations  (1)  and  (2)  is 


F 


1  n  L 


th(V2  -  Vj)  ^  P2  ^2  -  Pi 


This  force  acts  in  an  opposite  direction  to  that  of  the  flow.  The  in¬ 
ternal  thrust  is  thus  equal  to  the  difference  of  the  stream  functions. 
It  may  be  noted  that  the  relative  values  of  Vg  and  will  reflect  not 
only  the  difference  of  cross-sectional  area,  but  also  the  difference 
In  enthalpy  of  the  gas  at  the  two  stations.  Considering  the  ramjet 
engine  as  a  whole,  the  entering  fluid  will  be  air,  while  that  leaving 
will  represent  the  result  of  nixing  air  and  fuel.  An  additional  in¬ 
ternal  thrust  term  due  to  the  injection  and  acceleration  of  fuel  is 
therefore  given  by  The  magnitude  of  the  external  thrust  (or 

drag)  is  given  by  the  resolution  in  the  direction  of  flow  of  the 
forces  acting  on  the  projected  area  Ap- 

^ext  '  -  A2)  -  Po(Ap  -  A,)  =  -A^). 


The  sum  of  the  internal  and  external  thrusts  is  known  as  the  gross 
thrust  of  the  ramjet  engine: 


"e  '^e  "  \(Pe  '  ^0^  " 


*0  '^C 


=  -  "0  -  -  Aq) 

=  A^  Ped  +  -/m2)  -  Pq  A^(i  +  ^  y  m^)  , 


1 


1  ni* 


■  11.1  ■ 


IT.I'- 


i!  t  .!■  prcssui'i  .It  till  exit  Iroin  tl.o  nozzle-  differs  froia  the  free- 

•sti  i  in  V  UK- .  The  iiu  rease  of  momontuin  due  to  fuel  addition  may  be 

i.iktn  .uiuunl  of  by  i  n  t  roduc  i  :iK  the  factor  (1  y).  tihore  y  is  the 

!  I  e  1  -.1 11'  ra  tin: 

^  y>  *0  '^e  *  Ae^^’e  '  Pq^  '  *0  '"o  ' 

The  net  thre.st  is  loss  than  this  by  the  amount  of  t!.:-  ."xteriial  fric¬ 
tion  draft  force  -  tf  the  internal  friction  drag  is  con.sidered  to 

!'<-  sma  1 1  : 


F  -  F  -  D 
^net  g  ‘^ext 


Under  steady  flight  conditions, 


ext  ■ 


Fj^g^  is  observed  to  depend  strongly  on  the  ambient  pressure  and  the 
freestream  Mach  number  for  a  fixed  engine  geometry.  A  more  constant 
parameter  is  the  dimensionless  coefficient  obtained  by  dividing 
by  the  product  of  the  engine  reference  area  (nominal  cross  section) 
and  the  freestream  dynamic  pressure  Qq. 


C 


F 


(61) 


This  quantity  varies  only  slowly  with  Mach  number  and  ambient  pressure 
level.  Since  it  is  Independent  of  engine  size,  it  is  a  useful  measure 
of  comparison  between  engines  of  different  size  under  different  flight 
conditions.  The  expression  for  the  thrust  coefficient  may  be  trans¬ 
formed  by  the  introduction  of  the  mass  flow  equation  and  the  perfect 
gas  law.  Since 
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sinr»-  (or  af)v  urhitraiv  area  A  , 
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A  different  comparison  of  engine  performance  that  may  be  men¬ 
tioned  here  is  given  by  the  quotient  of  the  net  thrust  and  the  fuel 
weight  flow: 


1 


f 


(62) 


This  quantity  is  known  as  the  fuel  specific  impulse.  High  values  for 
signify  efficient  engines,  and  values  of  from  1000  to  2000  seconds 
are  obtainable.  The  reciprocal  of  is  clearly  the  specific  fuel 
consumption;  this  latter  is  often  stated  in  units  of  reciprocal  hours. 

It  is  also  useful  to  express  stream  thrust  in  terms  of  other 
flow  properties; 


'f 


Substituting  for  w, 


(1  +  Y 


M  g 


M^) 


:f  =  * 


R 

Ylf 


(63) 


r(.lb'  rcl.U'.on  will  ‘Mive  a>.  .1  ntasuif  iil  j<>t  reaction  if  the  total 
t  t '■  re  atnl  Mach  nur.lier  are  known  as  well  as  the  rate  of  flow. 

Since,  however,  it  is  difficult  to  measure  elevated  stagnation  tem¬ 
peratures  in  flowing  gases,  it  is  therefore  customary  to  evade  this 
bv  measuring  instei.d  the  sonic  Jet  reaction  if  which  would  be  that 
t'_lt  at  a  choked  nozzle.  The  quotient  of  the  Jet  reaction  if  invo 
the  weight  rate  of  air  flow  is  then  defined  as  the  air  specific 
impulse  S  • 

d 

1-  * 


S  is  a  parameter  widely  used  to  compare  the  performance  of  a  ramjet 
onplne  imjp'i  different  conditions.  Substituting  for  ‘i  the  expression 

in  Eq.  (63).  for  M  ”  1,  gives 


,/K 


Since  w  ”  +  y)t  then  from  Eq.  (64), 


Sa  =  (1  + 


V^t\l2R  (1  -i  y) 


Combining  Eqs .  (63)  and  (64), 


if  =  w  s 

“  “  M  V2(l  +  y)  K 


*a  Sa  • 


(64) 


(65) 

(66) 

(67) 


This  Mach  number  function  0(M),  according  to  Rudnlck  (Ref.  5),  thus 
relates  the  stream  thrust  of  a  Jet  of  Mach  number  M  to  the  sci. ic 


Cne-Dlaiensional 


v.iluf,  w,  ,  '>(  t!i<‘  lhrii<^t  (at  Iho  tctal  *  nniptraturc )  .  By  defi- 

I'lltlil.  .  t'K  voloi  c  . 


0{ll) 


'( 

-T 

‘t 


(68) 


This  function  is  also  shown  plotted  in  the  Appendix. 

It  ffisy  be  noted  that  both  vdeceleration  of  supeido.iiv  flow 
through  a  normal  shock  and  heat  addition  to  tlor  in  a  combustion 
chamber  of  constant  area  are  cases  that  Illustrate  constant  stream 
thrust.  Between  Stations  (1)  and  (2),  Eq .  (67)  gives 


for  constant  stagnation  temperature. 


*a  Sa  0(“)l 


*  0(1I)_ 

a  a  z 


whence 


0(ll)j  -  0(14)2  ,  (69) 

In  flow  with  heat  addittnn,  ditterent  values  of  S  must  be  used  corre¬ 
sponding  to  the  appropriate  inlet  and  exit  stagnation  temperatures. 
Using  subscript  (2)  to  denote  conditions  at  the  combustion  chamber  en¬ 
trance,  and  (b)  for  conditions  after  combustion,  then  for  y  -  0  and 

y  “  7/5,  and  using  £q.  (65),  then  for  T  In  °R,  and  ♦  In  Lb/sec, 

t  a 

0(M)2  _ 

0(1!).  -  2.382  -8 — .  (70) 

‘‘b  ^  2 

The  ratio  of  the  Mach  functions  is  thus  Inversely  proportional  to 

(S  /  i/TV  ) ,  parameter  which  will  frequently  be  reieried  to  oelow. 

^  ^2 

In  practice,  the  above  relationship  needs  to  be  moa..fled  by  the  in¬ 
clusion  ct  a  flameholder  drag  term  to  enable  it  to  apply  to  a  real 


'’"••-lUnn'ns  ;  T.a  1  iiai-. 


combust iun  chamber  if  combustion  is  stabilized  on  such  a  device.  The 
lirai?  force  acts  in  ••n  ttpposite  direction  to  tliat  of  the  flow  through 
the  romhustion  chamber  and  is  hence  to  be  deducted  fro*  the  inlet 
stream  llirust  when  the  exit  value  is  being  determined: 


-  ‘f  -  D.  . 

b  o 


Si'.rh  a  drag  term  is  generally  computed  in  the  form  of  a  coeffjelnnt 

Cjj  ,  obtained  as  before  by  dividing  by  the  product  of  the  combus- 

b 

tlon  chamber  area  into  the  combustion  chamber  inlet  stream  dynamic 
pressure  q„  ; 


2  **2  Ag  ^  P2  Mg 


This  drag  coefficient  is  not  strictly  constant  except  for  flows  at 
very  low  speeds  which  may  be  treated  as  incompressible. 


^2  ^  ^2  ^2 

*a  =  2.382  0(U)^ - ^ - .  (72) 


may  be  expressed  in  terras  of  P„,  M„ ,  and  T.  in  order  to  eliminate 
a  z  z  tg 

it  from  the  preceding  equation; 


^2  ^2  *2  o  *2  \  ^  ^2  “2 

\  *2  ^2  ^2  ^^**>2 - T - 

^  ^2 


whence 


X  ^  “2  • 

0(M).  -  2.382  0(M)„ - -  . 

2  2  «2 
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Alternatively,  using  Eq .  (72) 


K' 


1  -  (C„^  -  2)  r-f 


This  equation  for  covbustor  flow  cannot  be  further  Simplified 

unless  Cp  takes  on  such  values  as  0,  1,  or  2.  If  Cp  “  0,  Eq.  (70) 
b  b 

results.  If  Cp  -  1,  the  numerator  of  the  right  hand  side  of  Eq.  (74) 

^  2 

reduces  to  (1  +  y  H2/2) ,  which  is  approximately  equal  to(P^  ^^2^’ 

2  ^ 
higher  powers  o^  Mg  are  neglected,  so  that 


[s,  /  V'tT'  ~  -r - - — JTT 

L  *b  ^  (P/P^)  9 


If  Cq  -  2,  the  numerator  of  the  right  band  side  becoaes  unity,  and 

b 


For  each  of  these  three  cases  one  may  select  any  two  of  the  pararaeters 

(S^  /  y/T^ ) ,  Mp,  and  Mg.  in  order  to  solve  for  the  third  quantity 
h  2 

directly. 

The  total  pressure  after  combustion  may  thus  be  related  to  the 
inlet  total  pressure  by  means  of  the  continuity  of  either  mass  flow  or 
stream  thrust.  It  is  therefore  possible  to  write  an  expression  for 
strean  thrust  in  terms  of  Mach  number  and  total  pressure  across  the 
region  of  combustion.  Since 


(f/P.)  -  (1/P)  •  (P/P.)  , 

'2  2  ^2 


then  for  constancy  of  stream  thrust, 


!* '-  * '  ’i'l''’  *’''sI‘?5-'4 

i  1  ^  f  ^■'‘m 


M  P  ),, 

I’  p  ,  »  *“*  /''/«•• 

'!,  'p  !\‘r  • 

h 

i  '  ‘'ii'l  t  f  inpi  ■  c.'i  t  u  rc  aitor  (M^mlui.s  i  i  on  may  be  ob- 

i!'iini  l.otal  pi'essiire  aiul  tempei'ature  i  ti  a  manner 
'■  eM'il  lo  obtain  tile  valiK?  ol  itic  flow  after  a  normal 
■l>  'Plan  t  i  t  i  es  are  .seldom  of  latere.sL.  The  flow  after 
eomp 1 1.' t.e*  I y  clef  inc.'ci  by  a  knowlodtic^  of  the  parametei’s  of 
'.Ola  I  t  ('m])era  lure  (or  S,^ )  ,  and  mass  flow  (or  total  pre.i- 
I'r  I '.ssn  re  ,  Maeh  number,  and  area  of  the  exl'.aust  stream 
I  no/.a.io  niay  tie  cnalculated  from  these  parameters  as  pj-e- 

l  I.  e  < !  , 


best  available  copy 
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iwo  ISIMI  f.'SlfjtJAl,  SUPrffSGNIC  WAVi;  SYSTEMS 


;■  '  ni.  ili'.lit  VI' 1 1 II' 1  t  1  IS  ,  a  rain.li.'l  is  adcuriitaiilM.l  Uy  a 
i  I  n,  !  Ml  •ainl-all  rlna'l;  wave  sysli';:;,  ai-i  lilillcalecl  111  Fit-',  1. 

'  '  .  .1',  i:,,  1-11!  -I'l'  1)1'  il  ,i  Si' iissi.' M.  uiulcr  llii,'  lu.'acljiin  of  c:io  -  U  i  - 

‘  .  I'T  I  \ai;',p  1 1' ,  I  iiiw  t  lu'iuiiili  U;('  ol)  1  1  qun  bo'v  wave  sc-l  up  by 

■ . :Mal  I'  ai)  .tclual  ratiiji’t  d  i  f  lusei’  is  a  complicated 

'■  '  >'|  I  -d  1  mi  nn- 1 1 '.na  1  las'-  that  is  diasc  r  i  iK'ti  l;y  a  nonlinear  differential 
'I'n'ii".'.  SiMiii'iiiit  analysis  will,  however,  be  pi’escnted  to  illus- 
iiai'  'A'bai  II,  liinili.'il  tiy  olilique  shock  waves  and  by  Prandt  1 -Moyf'i' 

'  I'.i  in .  I  ora- ,  'inri'  'oaiii'  oi  Hu'sc  are  fundamental  to  the  internal  aero- 
I  >’ 11.1 '  i  1  IS  I)!  till  I'aMji'i  I'liylne,  and  to  the  case  of  'external'  eom- 
i'.ist  ion  as  li  1  sc  ussi 'll  Ivj'iefly  below. 


■'i  .  1  SIk'cIa  Waves 

■rile  s’liiploi'  case  of  oblique  .shod:  waves  ironorated,  e.g.,  by  a 
.'.I'liae  n  a  s'aperson  ic  si  rcaiii,  is  hero  considered  to  Illustrate  the 

nr  iiic  lilies  involved.  .As  before,  the  first  assumption  made  is  that  the 

wave  is  i  ii  f  i  n  i  t,  e  1  y  thin.  It  is  also  necessary  that  the  wave  be 

'  a  1  i  at:  bed  ■  to  the  point  of  change  of  se  CviOTi  uwC  boundary  and  hence 

ti'.ai  it  is  plane.  This  is  because  a  rigorous  theory  for  detached  flows 
tioos  not  exist.  A  conical  case  gives  rise  to  curved  waves,  and  this  is 
also  true  of  a  detached  wave  which  results,  for  example,  when  the  angle 
'  in  Fig.  11  i;--  large  at  low  Mach  numbers.  In  general,  flow  downstream 
d.  :'.n  oblique  wave  is  supersonic,  but  it  is  possible  to  obtain  subsonic 
'.  1  ow  itsst  bciore  6  becomes  large  enough  to  cause  detachment. 


It  is  convenient  to  commence  by  considering  how  the  flow  is 
iiiriicd  by  the  action  of  a  wetlge  of  angle  d  projecting  into  the  stream. 
I'r.i  inciin.it)  on  U  of  the  oblique  •.v;'.ve  sol  Uf)  by  a  wedge  Is  greater 
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of  a  shock  wave  Is  greater  than  ihal  of  a  Mach  wave.  The  flow  vectors 
and  V.,,  representing  the  flow  on  entering  and  leaving  the  shock 
wave,  are  shown  resolved  into  normal  and  tangential  components  along 
the  wave.  V2  is  acttially  obtained  by  vectorial  addition  of  a  tangen¬ 
tial  component  V_  (‘  V_  )  to  the  normal  velocity  V„  .  The  relation 

*2  *1  ”2 

between  V„  and  V„  is  that  for  flow  passing  through  a  normal  shock 

"1  "2 

so  that  V.  is  Guhs-inic  These  quantities  may  also  be  related  by 

"2 

means  of  the  continuity  equation  normal  to  the  wave  front, 


'1  'N, 


2  ’N„ 


(V« 


dV  ^  f' 
N  '■ 


dp  ) 


whence 


V„  dp 
^1 


-  pj  dVjj,  so  that  dVjj 


V  ^ 

\  Pi 


Application  of  the  momentum  equation  parallel  to  the  wave  front  gives 


(^1  ’^Tj 


-  Cpj  +  dpXVj,^  +  clVjj)(V.j.^  +  dV.j,)  , 
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Fion  iho  momentum  equation,  for  the  normal  component  of  the  shock, 


P,  ♦  Pi  Vn  ^  P->  *  Po  V, 


2  >^2  Ng  ’ 


so  that 


Vv"  -  V  2 
N ,  No 


(P„  -  P  )  (^  4 
^  1  Pj  P2 


Cl)  r.l)  1  nin;4  Eqs  .  (7.;)  and  (7E).  one  has 
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Tlicse  are  well-known  Rankine-liugoniot  relations. 
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IntroUuc  inti  the  Mach  numbers  on  the  two  sides  of  the  wave, 


=  V  sin  9  =  u,  M  sin  0  (tiS) 

‘1  "  -  ^ 

and 

Vjj  -  sjn{0  -  6)  -  =^2  ''’2  -  <5)  .  (01) 


From  the  momentum  and  continuity  equations, 
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so  that 
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'r  +  1)  Pg  +  (>'  -  1) 
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(Y  +  1)  Pj  +  (y  -  1)  P2 


by  substitution  from  Eq .  (80). 
From  Eqs .  (83)  and  (84), 


Mj  sln^  8 


y  +  1  **2  .  y  -  1 
“  2y  ~W  ’ 


and 


sin(9  -  5) 


y  ->■  1 
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(85) 
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Standard  relationships  thus  exist  between  Mj ,  ,  *,  0,  (Pg/Pj^)  and 

(pg/Pf).  Some  of  the  most  important  are  shown  in  Fig.  15,  in  which  it 
will  be  seen  that  for  given  values  of  and  5,  two  solutions  of  0 

are  possible;  experimentally,  however,  it  is  observed  that  only  the 
lower  or  weaker  wave  solution  occurs.  .4bove  the  line  marked  0-0 

max 

the  shuck  wave  detaches.  Although  the  constant  value  of  1.40  has  been 
used  for  y  in  plcttlng  this  figure,  for  Mach  numbers  higher  than 
about  3  the  temperature  rise  through  the  shock  wave  will  necessitate 
account  being  taken  of  variable  specific  heat.  Similar  effects  occur 
in  the  case  of  flow  through  a  three-dimensional  oblique  shock  set  up 
by  a  cone  facing  into  a  supersonic  stream.  In  this  case  stream  prop¬ 
erties  are  constant  along  radii  from  the  vertex.  The  compression 
effect  of  a  cone  is  weaker  than  a  wedge,  however,  and  the  oblique 
shock  becomes  detached  at  a  lower  Mach  number. 


eo 


(b) 


Fig.  15(b)  Voriotion  of  Shock  Wove  A.igle  with  Upstreom  Moch  Number 
for  Various  Flow  Deflection  Angles 
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IL  -i.i  .,i  ill  ll.i.i  li  ciiiisidcr  briclly  the  rcliec- 

lion  and  i  n  li.  rs*.  c  t  iun  ol  oblique  nbock  waves.  I’igurc  lC(a)  shows  the 
iimple  I'a.-it  of  ideal  reflection  of  a  shook  wave  at  a  solid  boundary 
whorebs  the  dovinstreani  flow  is  turned  parallel  to  the  boundary.  A 
similar  e.ise  is  shown  in  Fig.  16(b)  for  the  intersection  of  shook  waves 
of  equal  strciij-th.  In  Fig.  16(c)  the  effect  of  the  intersection  of 
shock  waves  of  unequal  .strength  is  Shown:  pressure,  temperature,  den¬ 
sity  and  entropy  discontinuities  occur  while  the  t'.vo  flow  voloci:.y 

vectors  remain  parallel.  If  the  wedge  angt<-  5  is  fixed  and  is  de¬ 
creased,  a  point  is  reached  when  it  is  no  longer  possible  to  deflect 
the  M2  stream  through  the  angle  -  parallel  to  the  wall.  This  case, 
shewn  in  Fig.  16(d),  is  known  as  a  'Mach  reflection.'  A  similar  sur¬ 
face  of  di ecnnt i nil ity  is  Set  ap.  Th<  corresponding  situation  brought 
about  by  the  Mach  reflection  of  two  shocks  of  equal  strength  is  shown 
in  Fig.  16(e).  The  above  case,  Fig.  16(a),  of  'positive'  reflection 
at  a  solid  boundary,  may  be  compared  to  the  case  of  'negative'  reflec¬ 
tion  a*  the  free  boundary  of  a  let,  as  shown  in  Fig.  19,  which  is  the 
P-a-  .  .ever  expansion  case  dealt  with  in  the  next  section. 


3.2  Prandtl-Meyer  Expansions 

in  contradistinction  to  subsonic  flow,  supersonic  flow  may  ex¬ 
pand  around  a  sharp  corner  without  loss,  and  unlike  the  shock  phe¬ 
nomena  uescribed  above,  the  expansion  wave  is  an  Isentropic  flow  case, 
a.nd  can  be  considered  to  be  a  succession  of  Mach  waves.  The  two- 
dimensional  case  shown  in  Fig.  17  vill  now  be  referred  to.  A  turning 
angle  dv  is  shown. 

The  equation  cf  continuity  may  be  written  down  as  for  the  ob¬ 
lique  shock  case,  giving 


dVj^ 


Pi 


(87^ 
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ami  from  thf  morn.' ii turn  equation  applied  parallel  to  the  shock  front, 

V„  -  V,_  ,  (88) 

1  *2 

as  ’'ofore. 

From  the  equation  of  momentum  normal  to  the  wave  front, 

dP  =  P^-P^  -  (Pj  V^^)  -  (p^  t  dp)(V^^  H  ..  dVj^.),(89) 

whence 


o, 

•  i 


2  dP  2 

Eliminating  dv^  between  Eqs .  (87)  and  (89),  ^  ®'''” 

phasizes  that  this  is  an  isentropic  process. 

From  the  Mach  triangle, 


-  Vi  =  dV  =  dVj^  sin  /i^ , 
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Now  fron  the  energy  equation 
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sin'^  (y  -  -[r  -  4 


1?  ^  -  4' 


(92) 


Siiu'<‘  Mj  ami  V'^  are  re  lated  by  Eq  (90),  Eq .  (92)  may  be  transformed 
Uilo  the  fol 'i(.)*'lng  form: 


A  iTiareiraura  value  of  v  occurs  when  M,  ^  cd  (assuming  that  the  flow  is 

L 

still  gaseous )  ,  i  .  c  .  . 


=  130.5 


for 


In  practice,  flow  will  not  expand  around  a  sharp  corner  to  anything 
like  this  value,  and  boundary  layer  separation  limits  v  to  about  30  , 
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.'i.3  Nozzle  Flow  and  Jet  Structure 

The  IsentroDtc  equations  for  variable  area  duct  flow,  togatuer 
with  the  shock  wave  equations  derived  above,  are  also  useful  for 
Illustrating  various  phcnosena  occurring  in  nozzles.  For  example,  it 
wa.";  shown  in  Eq .  (21)  that 


I.C.,  the  change  of  pressure  with  area  depends  on  the  Mach  number  of 
the  flow.  When  M  <  1 ,  dP  has  the  same  sign  as  dA.  but  when  M  >  1,  dP 
has  the  opposite  sign,  the  change  occurring  at  a  throat  where  M  -  1. 

A  nozzle  attached  to  a  reservoir  may  thus  consist  of  a  subsonic  por¬ 
tion  (converging)  and  a  diverging  portion  which  will  run  partially  or 
wholly  supersonic  If  the  reservoir  pressure  is  sufficient  to  cause 
choking  at  the  throat.  As  may  be  seen  In  the  two-dimensional  diagram. 
Fig.  18,  if  the  pressure  P  as  obtained  at  the  throat  during  flow  is 
initially  such  that  (P/P^)  >  0,528,  the  critical  pressure  ratio  for 
air,  the  nozzle  corresponds  to  a  subsonic  venturi,  and  the  pressure 
variation  follows  a  line  above  ABC.  The  pressure  downstream  of  the 
nozzle  then  corresponds  to  a  certain  exit  pressure  ratio  (Pg/P^)- 
When  the  pressure  ratio  (P/P^)  -  0.528,  two  flow  solutions  are  theo¬ 
retically  possible,  the  pressure  variation  lying  along  either  ABC  or 
ABD,  for  flow  In  the  diverging  portion  to  be  either  completely  subsoiilu 
or  supersonic.  Which  one  is  followed  depends  on  the  value  of  the  down¬ 
stream  pressure  ratio  (P^/P^),  and  acceleration  of  the  flow  beyond  the 
subsonic  condition  may  only  be  brought  about  by  reduction  of  ^Pg/P^) 
below  the  value  at  C,  by  varying  either  component  of  this  ratio. 

(If  P^  remains  constant,  the  mass  flow  through  the  taroat  cannot  then 
be  Increased.)' 

It  may  be  noted  that  the  surface  In  the  flow  at  the  throat 
where  M  becomes  unity  Is  not  plane,  but  Is  slightly  convex  down- 


■  '  ;  -  :  !-  '  ■  . . .  k  >  ■  -  . ,  If.,  i.f.  ii.'ua  1  flow  is  an  ovor- 

'  1 !  i  !  I  <  .1 1  1  oi,  .  flu-  ''tno.iih  no^ilo  ciuitonr  ''hown  in  FIk-  IS  is  not  tht- 
|  ,11  1  ril  lor  acof  11  1  at  inu  the  flow  lo  suporsonic  vf- 

I  '.'i'.s',  •.asi'  •■lio'o.n  lo  illustrate  the  various  flow  con- 

liiln  :s.  luc  ill  s  i  a  luii  on  I  hi-  iu-incip’o  of  avoia.ince  of 

"h.icl';  lossi  s  by  sia  ci'ss i  vo  oaiii'o  '  I  a  t  ion  rif  compression  and  expansion 
w:i\'"s  ,  an^i  by  t,-'  '*  vli  oif,  !  aoarly  ono-'J imens lona  1  flow  can  be  obtained 
:i  1  111.’  ;  .1  r  1 1  ! ‘.'I ]  exit  section  of  ihe  nozzle  s’uowii. 

For  supirsonie  flow  to  occur  downstream  of  the  sonic  throat, 
the  pressure  and  Maeb  inimbcr  reached  are  now  fixed  by  the  expansion 
ratio  (A  A  )  Since  for  supersonic  flow  to  occur  throughout  the  noz¬ 
zle.  the  exit  ratio  must  correspond  tc  the  point  li;  for  all  values  of 
a  Pp,  the  solution  will  be  represented  by  supersonic  followed 

by  subsonic  flow  in  the  divergint;  portion  of  the  nozzle.  In  the  ab¬ 
sence  of  a  converging  section  to  act  as  an  isentropic  diffuser  to  bring 
about  this  transition,  the  change  must  occur  nonisentropically  through 
a  shock  wave,  its  position  being  completely  determined  by  the  pressure 
ratio  across  the  shock  with  the  downstream  subsonic  pressure  rise  fol¬ 
lowing  one  of  the  paths  BEFG,  in  order  that  P  ”  Pq  ,  say,  I  f  P^  “  Pjj , 
the  shoe!:  in  its  simplest  form  will  be  a  normal  shock  positioned  at 
the-  nozzle  exit.  As  reduced  further,  recempression  of  the 

flow  will  occur  through  oblique  abocks,  either  by  Mach  reflection 
(P^,  near  Pjj)  or  simple  reflection  (P^  near  Pp).  These  are  represented 
by  the  overexpanded  cases  shown  in  Figs.  19(a)  and  19(b)  in  somewhat 
idealized  form.  The  supersonic  core  of  the  whole  jet  gradually 
diminishes  in  cross  section  as  the  subsonic  jet  boundary  extends  due 
to  entrainment  of  the  surroundin;  atmosphere,  when  present.  A=  is 
’educed  to  Pp ,  the  external  shocks  become  weaker,  but  never  entirely 
disappear,  even  with  a  parallel  jet  emerging  from  the  exit  when  P 

© 

Pp.  Since  Pp  is  fixed  by  the  reservoir  pressure,  further  reduction  of 
P^  will  result  In  the  ."ernation  of  expansion  waves  at  the  exit,  as 
shown  in  Fig.  19(c).  This  is  the  underexpanded  case. 

So  far  no  mention  has  been  made  of  viscous  effects,  but  in 
actual  flow  the  pressure  rise  through  the  shock  wave  may  give  rise  co 
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Fig.  18  Flow  with  Normol  Shock  In  Air,  M  i  it  Lovol  Nozrle 

boundary  layer  separation  on  the  walls  of  the  nozzle,  with  resulting 
flow  downstream  as  shown  in  the  upper  diagram  of  Fig.  IS.  As  an  ex¬ 
ample  of  this,  It  may  be  noted  that  In  the  rocket  (straight-sided) 
type  of  nozzle  referred  to  below,  separation  of  the  boundary  layer 
will  occur  when  the  nozzle  exit  pressure  Is  approximately  0.46  of  the 
ambient  pressure,  whatever  the  divergence  angle  may  be.  More  compli¬ 
cated  shock  wave-boundary  layer  Interactions  are  also  possible  but 
need  not  be  considered  In  detail  In  this  treatment.  A  typical  case 
Is  the  bifurcation  Os  the  shock  at  the  nozzle  wall  (Mach  I’of lection) . 
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An  actual  supersonic  let  emerKing  from  a  nozzle  may  saow  con¬ 
siderable  fine  structure,  but  it  consists  basically  of  a  number  of 
'cells'  representing  cyclic  flow  transitions  which  persist  until  they 
arc  damped  out  by  viscous  interaction  with  the  surrounding  atmosphere. 
The  Jet,  after  it  leaves  the  nozzle,  exhibits  successive  negative,  or 
opposite,  reflections  of  shocks  or  expansion  waves  as  required  to  turn 
the  flow  at  a  free  {.constant  pressure)  boundary,  as  shown  In  Fig.  19. 
The  successive  subsonic  cores  behind  the  Mach  reflections  In  the  over- 
expanded  case,  Fig.  19(a),  are  gradually  energized  by  turbulent  mixing, 
so  that  the  surfaces  of  discontinuity  die  out  rapidly. 

It  is  convenient  to  illustrate  such  phenomena  with  reference  to 
two-dimensional  nozzles  and  Jets,  since  the  flow  patterns  are  thereby 
simplified;  for  example,  the  Jet  boundaries  are  initially  straight 
lines.  In  the  more  common  axlsymmetrlcal  arrangement,  these  boundaries 
become  curved.  Such  a  convergent  nozzle  perating  under  high  reservoir 
pressure  gives  rise  to  a  heavily  underexpanded  type  of  Jet ,  as  shown  iri 
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the  three-dimensional  case  of  Fig.  20,  and  uuch  a  Jet  becomes  subsonic 
after  one  tell.  The  form  taken  by  the  jet  under  various  pressure  con¬ 
ditions  has  been  diuoussud  by  Wilcox  e_t  a_l .  ,  (Ref.  6),  among  others. 

All  these  relations  refer,  as  Indicated  previously,  to  equilibrium 
flow  conditions. 


Up  to  now  only  completely  expanded  nozzle  profiles  have  been 
considered,  l.e.,  those  in  which  the  flow  emerges  parallel  to  the 
axis  if  the  nozzle  is  operated  under  the  correct  pressure  conditions. 
Exit  nozzles  used  with  many  ramjet  engines  are  customarily  constructed 
with  a  conical  diverging  section,  on  the  pattern  of  rocket  nozzles, 
the  shortening  representing  a  compromise  between  reduction  of  weight 
and  supersonic  friction  loss  on  the  one  hand,  and  increased  divergence 
loss  on  the  other.  The  most  suitable  expansion  angle  will  thus  be 
somewhat  greater  than  that  at  which  the  thrust  is  a  maximum.  This 
simplified  contour  is  clearly  an  advantage  if  the  engine  is  designed 
for  a  limited  life.  The  Jet  emerging  from  a  conical  nozzle  Is  differ¬ 
ent  from  that  previously  considered  in  that  divergence  imparts  a  radial 
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componLi.  I  of  volocilv  to  t!ie  jel  ,  which  vlll  then  continue  to  expand 
..  ....!  i'l  lilt-  w'K're  recoinpi  es»s»  ton  pets  in. 


J.4  ''Supersonic"  Combustion 

In  view  of  tne  above,  it  is  of  Interest  to  refer  briefly  to  the 
modifying  effect  of  heat  addition  (by  burning)  to  a  supersonic  etrpnii< 
n  nr  frcc  jct.  (Th Is  caso  is  of  some  Importance  in  ad- 

vancea  ramjet  design,  since  supersonic  flow  in  the  diffuser  and  com¬ 
bustion  chamber  could  be  a  desirable  mode  of  operatl.n.)  A  theoreti¬ 
cal  ariilysis  of  such  a  concept  of  burning  in  a  supersonic  stream  In 
ramjet  engines  has  been  published  by  Weber  and  Mackay  (Ref.  7).  It 
has  beei.  previously  shown  that  one  of  the  immediate  effects  of  such 
heat  addition  is  the  reduction  of  the  Mach  number,  either  through  a 
shock  system  or  otherwise.  Three  recent  demonstrations  of  combustion 
in  a  supersonic  stream  are  shown  in  Figs.  21  to  23  (Refs.  8  through  10). 


Fig.  21  Strong  Standing  Delonotion  Wove 
( Two-Dimensional  Mach 
Reflected  Shock) 

(  T cketl  from  Rrt  ?.  ) 


rig.  22  Co.r.busfion  in  Open  Jet  from 
'Jn''ere»ponded  Axisymmetric 
Noxxie 

(Token  hot"  Ref.  9.) 
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Fig  23  Cambinotion  Open-Shutter  onj  Schlieren  Flash  Pliutngrnph 
of  Flame  ond  Aasocicited  Shack-Wave  Systt .  in  Stream 
belaw  Surface  of  Wing,  Moch  Number,  2.47,  Angle  af  Attack,  2 

(Tskei  from  Rof.  ]0, ) 


In  the  first  of  these  caccr  a  “ach-reflectert  nnrmai  shock  wave  was  »et 
up  by  inserting  wedges  in  a  wind  tunnel.  The  rise  of  static  pressure 
and  temperature  downstream  of  the  normal  part  of  the  shock  system 
ignited  the  combustible  gas  mixture  (H2  -  air)  at  this  point.  The  on- 
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s(  t  lii  i  i'n‘)U‘;t  ;  >ii  '‘■MS  otisfrvoJ  to  displnfo  the  normal  shock  upstream 
(Fii:;.  21).  In  the  second  Instance,  Fig.  22,  a  similar  situation  was 
sol  up  downstream  of  the  first  Mach  disc  occurring  in  an  underexpanded 
open  jet,  and  it  was  considered  that  the  axial  distance  between  the 
flame  front  and  the  Mach  disc  corresponded  to  that  expected  from  igni¬ 
tion  delay.  Interaction  between  the  combustion  region  and  the  wave 
system  was  also  observed  in  some  cases,  resulting  in  the  wave  system 
moving  upstream,  corresponding  to  a  reduction  in  Macn  number.  Neither 
of  these  two  cases,  however,  illustrates  true  supersonic  burning,  but 
they  are  considered  to  represent  strong  standing  detonation  waves. 

To  obtain  combustion  of  a  liquid  fuel  in  a  supersonic  stream 
without  the  formation  o  A  iiC-ji'iiiciX  shoe  k  waves.  It  might  at  first  sight 
seem  desirable  to  inject  the  fuel  into  a  supersonic  air  stream  with  a 
component  in  the  air  flow  direction  equal  to  that  of  the  air,  but  this 
is  clearly  out  of  the  question  because  of  the  enormous  pressures 
necessary  to  effect  this.  In  fact,  demonstrations  to  date  have  taken 
the  form  of  simple  injection  at  right  angles  to  the  flow,  and  Fig.  23 
typifies  a  case  of  external  burning  of  a  liquid  fuel  beneath  a  wing 
section.  Here  combustion  would  appear  to  be  fairly  generalized  in  the 
supersonic  region  downstream  of  the  oblique  shock  formed  at  .he  point 
of  injection.  These  and  other  studies  indicate  that  this  application 
of  external  supersonic  burning  may  be  important  in  ramjet  propulsion 
at  high  Mach  numbers.  A  recent  report  by  Woolard  (Ref.  11)  lisis 
further  useful  references  to  published  work  in  this  field. 


IV  RAMJET  COMPONENT  PERFORMANCE 


iho  above  gas  uynaauc  l uiisi I  i  or.s  ao'v  need  to  be  applied-- 
with  emphasis  on  component  port ox‘mai\ce--in  order  to  illustrate  the 
tancticninK  of  a  ram, jet  as  a  heat  engine.  Mention  will  be  made  of  the 
ccsign  and  The  detailed  functioninR  of  the  difluser,  combustor,  and 
exit  nozzie  and  how  these  are  affected  by  flight  environmental  condi¬ 
tions  a  Mil  the  requirements  of  thrust  and  power  output. 


4.1  The  Ram,jet  Thermodynamic  Cycle 

The  ideal  thermodynamic  operat.ing  cycle  of  the  ramjet  is  shown 
in  Pv  form  in  Fig.  24.  Strictly  sne.ahin't,  such  a  diagram,  which  here 
represents  a  continuous  flow  process,  should  not  close,  although  it  is 
conveniently  represented  in  this  form  for  purposes  of  ana’ysis.  This 
cycle,  the  Brayton  or  Joule  cycle,  consists  of  four  parts,  three  of 
which  are  associated  with  the  engine  proper,  while  the  fourth  is  assumed 
to  be  completed  in  the  external  atmosphere.  The  first  part  of  the  theo¬ 
retical  cycle  consists  of  adiabatic  compression  of  air  which  enters  at 
high  velocity  and  low  pressu'e,  as  shown  by  a — ►b .  Heat  is  then  added 
under  conditions  of  constant  pressure  P,  and  constant  momentum  (P  + 
pV^'  =1  constant).  This  operation  increases  the  entropy  of  the  gas,  and 
there  is  an  increase  o'"  velocity  as  the  gas  density  is  reduced.  This 
is  shown  by  b— ►c.  The  heated  gases  are  then  further  expanded,  trans- 
forti  ing  a  part  of  the  added  thermal  energy  into  mechanical  energy. 

This  portion  of  the  cycle  is  shown  by  c — »-d.  This  exhaust  Jet,  even 
if  expanded  to  the  freestream  value,  emerges  with  a  much  higher  tem¬ 
perature  and  velocity  tha..  the  entering  air  flow.  TLa  force  required 
to  accelerate  the  expanding  exhaust  gases  thus  reacts  against  the  in¬ 
terior  surfaces  of  the  ramjet  and  thereby  propels  tae  engine  forward. 

In  the  final  phase  of  the  heat  engine  cycle,  assumed  to  occur  in  the 
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Fig.  24  Thermodynomic  Cycle  of  a  Ramjet  Engine 


free  n tmosijhere ,  the  hot  exhaust  gases  are  cooled  by  contact  with  this 
relatively  infinite  reservoir,  with  consequent  increase  in  density  at 
constant  pressure,  as  shown  by  d — »a. 

Since  the  ramjet  engine  is  a  continuous  flow  device,  the  func¬ 
tions  of  compression  of  the  inlet  gas,  additicn  of  heat,  and  expansion 
cf  the  hot  gas  need  to  be  performed  by  separate  components;  the  de¬ 
vices  which  carry  out  t.hese  functions,  the  diffuser,  combustion  chamber, 
and  exit  nozzle,  were  shown  in  Fig.  1.  The  changes  in  gas  state  from 
Station  0  (located  in  the  free  air  stream)  to  Station  2  are  accomplished 
by  the  diffuser  and  are  characterized  by  partial  conversion  of  the 
kinetic  energy  of  the  gas  into  pressure  energy,  usually  through 
both  oblique  and  normal  shocks  followed  by  diffusion  in  the  internal 
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i  I  111  ri  ■  !>  1  .  t  ap  i3>  t ho  >i  the  center  body 

*  . .  ■  1  if.  t  'I'w.ii'  )?  'i<’  liil.  t:  this  leature  Is  characteristic 

;  1  i  i  <  -,1^!.-.  iipoi  at  Ir.^  ;'-T  vepers jiiio  Math  n>inbers  in  excess  of  about 
i  f.  V.ii  li'us  iuss“-s  reduce  the  pi  •ssurc  obtainable  at  the  end  of 
!if  fusion,  ,111(1  I -e  (.ilffu.scr  total  jiressure  recovery,  which  reflects 
tile  1  fflcioiiry  of  tlu  ipresslon  process,  decreases  as  the  flight 
Min-t!  number  is  lucre.  Aithouifh  On  the  actiia.l  Pv  diagram  (Fig.  24) 

t  h,.'  cemprcr.sion  process  may  seem  to  conform  closely  to  the  ideal  ad<p- 
batit  case,  a  l.s  diagram  would  show  that  the  process  of  shock  compres¬ 
sion  of  air  IS  .It  tended  by  some  considerable  increase  of  entropy  due 
to  the  shock  transitions  and  the  diffusion  losses  mentioned. 

■Addition  of  thermal  energy  to  the  ai:  in  the  combustion  chamber 
betM’-^  o-  Stations  2  and  b  is  generally  brought  about  by  the  oxidation 
of  .a  fuel  (although  it  could  equally  well  be  effected  by  other  chemi¬ 
cal  reactions  or  by  a  physical  heat  source  such  as  a  nuclear  reactor). 
As  previously  shown,  the  addition  of  heat  to  a  moving  air  stream  in¬ 
creases  the  velocity  of  the  stream  and  thereby  produces  an  effect 
analogous  in  some  ways  to  that  which  occurs  when  the  stream  passes 
through  a  shock  wave;  an  increase  of  entropy  together  with  a  fall  in 
total  pressure.  In  the  representation  of  the  combustion  process  on 
the  Ts  diagr.am,  if  this  .actually  followed  a  line  of  constant  pressure, 
it  could  be  represented  on  such  a  line,  as  shown  in  Fig.  12.  In 
practice,  tiamchoJder  drag  will  contribute  to  a  drop  in  pressure  dur¬ 
ing  the  combustion  process.  However,  it  seems  plausible  that  a 
Haylelgii  lino  for  constant  mass  flow  would  be  a  faithful  reproduction 
of  the  boat  ■‘.ransfer  procc-is ;  Fig.  25,  which  is  an  amplification  of 
the  previous  case,  shows  Rayleigh  lines  for  different  duct  areas  which 
thus  enable  ihc  whole  dilfasion  process  to  be  represented.  Referring 
to  Fig.  2b,  the  thormodynamlc  efficiency  17^  is  thus  given  by 
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Fig.  25  Rsmjet  Cycle  on  Ts  D'ogroni 
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If  the  propulsii  ^  .4-  ;  .  rt.  i-ccoants  for  the  ratio  of  useful  engine 
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Since  acceleration  ol  the  exhaust  gases  to  sonic  or  supersonic 
velocity  is  a  function  of  the  exhaust  nozzle,  the  formation  of  a  ther¬ 
mal  throat  at  the  exit  to  the  combustion  chamber  is  avoided.  The  com¬ 
plete  pressure  drop  in  the  nozzle  v.ill  not  be  realized,  however,  owing 
to  losses  due  to  friction  and  coclir;;.  Typ’r’al  average  gas  state®  at 
different  points  of  the  system  are  shown  in  Fig.  26.  These  diagrams 
are  reasonably  representative  for  an  engine  operating  at  a  flight  Mach 
number  of  about  2. 

The  functioning  of  the  various  engine  components  will  now  be 
discussed  in  somewhat  greater  detail.  For  convenience,  flo?*s  will  be 
discussed  with  reference  to  a  system  of  axes  chosen  to  move  with  the 
vehicle,  i.e.,  the  engine  is  assumed  to  be  stationary  wJrh  the  gases 
in  motion- 
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Fig.  26  Pressure,  Tenperulure,  Veloeify,  and  Density  at  Vaiious  Stations 
in  u  Ramjet  Engine 


(  Totten  from  Zucrow,  Aircraft  ond  Mitkiie  Propulsion, 
Vol  1,  p.  98.) 


.  Tire  D.'l  fuser 


This  is  generally  the  forwardmost  component  of  the  ramjet  engine 
and  ser”e-i  to  convert  the  high  velocity,  low  pressure  energy  of  the 
free  air  stream  into  low  velocity,  high  pressure  energy  to  meet  the  re- 
quirefiieiits  of  the  doonstrean  combustion  chamber.  Since  no  heat  is 
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t  r:>  ti'i  f  •■rrcii  ai  inis  stapp,  th**  forprcssion  from  0  to  2  is  adiabatic, 
dul  me  1  low  xs  not  isentropic^  because  of  shock,  skin  friction,  .xnd 
other  losses  due  to  turbulence  in  the  flow.  Static  pressure  and  static 
temperature  are  also  increased  between  Station?  0  and  2  by  cOiMpieBsiun, 
while  the  velocity  Is  reduced.  > low  within  the  diffuser  duct  from 
Stations  i  to  2  may  be  either  partly  subsonic  and  partly  supersonic, 
or  else  wholly  subsonic,  depending  on  the  location  of  the  normal  shock  ' 
wave  attached  to  the  center  body.  If  this  is  upstream  of  the  'rlet  1, 
the  diffuser  is  said  to  be  operating  'suber itically '  or  'with  spill¬ 
over.  '  whereas  if  it  is  downstream  of  the  inlet  1,  the  diffuser  is 
.said  to  be  operating  'supercr  itically  ’  or  ’with  shock  swallowed'.  If 
the  normal  shock  is  located  at  1,  the  condition  is  then  said  to  be 
that  with  'shock  on  the  rim’  or  'critical.'  These  definitions.  Illus¬ 
trated  in  Fifr.  27,  do  not  imply  that  the  diffuser  performance  is 
satisfactory  or  otherwise,  but  merely  indicate  the  relationship  be¬ 
tween  actual  and  design  operating  conditions.  Diffuser  total  pressure 
recovery  generally  tends  to  a  maximum  at  the  critical  condition  al¬ 
though  this  is  by  no  means  always  true. 

The  diffuser  exit  Mach  number  is  determined  by  the  combustor 
impedance  since  all  the  air  flowing  through  the  diffuser  (except  that 
diverted  to  power  auxiliaries)  continues  on  into  the  combustion  chamber 
at  the  subsonic  speed  determined  by  downstream  drag,  heat  addition, 
and  exit  nozzle  throat  area.  The  Independent  operational  variables 
determining  diffuser  performanc -■  are  thus  the  fixed  valu<;s  of  free- 
stream  and  exit  Mach  numbers.  These  variables,  together  with  that  of 
the  diffuser  total  pressure  recovery  ratio,  fix  the  two  diffuser- 
dependent  flow  vai fables,  inlet  area  capture  ratio,  and  Inlet  drag 
force  coefficient. 

Diffusers  which  utilize  either  multiple  shock  or  continuous 
(Isentropic)  supersonic  compression  experience  smaller  total  pressure 
losses  in  decelerating  supersonic  streams  to  subsonic  velocities  than 


^It  may  be  noted  that  it  is  also  possible  to  maive  use  of  a  center  body 
of  gradually  Increasing  curvature  ('isentropic  diffuser')  which  does 
not  give  rise  to  a  strong  oblique  shock  in  the  air  Inflow  region. 
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Fig.  27  Subcritical  (a),  Critical  (b),  and  Suptrcriticol  (c) 
Diffuser  Operatian 


do  those  which  effect  coapresslon  by  means  of  a  single  normal  shock. 

The  superiority  of  the  former  type  over  normal  shock  Inlets  is  small 
in  the  freestream  Hach  number  region  from  1.0  to  1.6,  but  this  superi¬ 
ority  increases  rapidly  thereafter,  and  makes  such  diffusers  essential 
at  Uach  numbers  above  2.0.  There  are  also  appreciable  total  pressure 
tn  the  diffusion  process  from  a  high  subsonic  Uach  number  down¬ 
stream  of  the  last  diffuser  shock  down  to  the  low  subsonic  Uach  number 
reqi.ired  at  the  combustion  chamber  Inlet.  The  pirtiire  is  further  com¬ 
plicated  by  the  possibility  of  the  diffuser  operating  out  of  alignment 
with  the  oncoming  alrstream.  The  presence  of  yaw  or  pitch  gives  rise 
to  an  asymmetric  shift  oi  ihe  shock  wave  pattern  generally  recultlng 
In  a  reduction  cf  the  total  pressure  recovery  and  air  capture,  particu¬ 
larly  If  flow  separation  also  occurs.  Normal  shock  diffusers  experi¬ 
ence  a  small  (and  nearly  predictable)  loss  In  mass  flow  and  pressure 
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■:  nora!  ''oiis  1  iUmoi 1 1  oiis  rcvlatln,.;  to  liibtabllity  have  been 
.!  hiy  !  (  .11  Cl.  ..lUif  .  Id)  ,  and  tlio  sub.Joct  is  also  covered  in 

'  ' ''  d  i’d-'I  of  tills  scries.  Pearce's  stability  crl- 

-1:''!.  i.v  that  the  s.Uipc  of  the  curve  of  diffuser  pressure  re- 
ersus  iiias'-^  flow  is  .loss  than  the  .slope  of  the  curve  of  combus- 
iiiiet  pre.sfuirc  versu.s  mas.s  flow.  If  this  requirement  is 
system  is  self-compcnsi'tln;:,  since  any  fall  in  diffuser  air 
li  .11  inns  , iboi.it  a  ni'oatoi'  fall  in  the  combustion  chamber  pres- 
n  llKii  wliich  occurs  in  the  available  diffuser  pressure  so  that 
i.niiiediaii'ly  comes  to  a  condition  of  .stable  equilibrium.  On 
IS,  .sui'p,'  cannot  occur  in  a  diffuser  which  provides  for  in- 
pressure  recovery  with  Increased  spillover  in  the  working 


Accordingly,  a  safe  criterion  for  stability  is  that  the  slope 
c'.irvci  of  diffusc-r  pressure  versus  flow  at  the  maximum  thrust 
point  lor  the  engine  shall  be  negative,  since  the  slope  of  the 
ol  comhuHiur  pressure  ver.sus  flow  is  then  always  positive.  The 
cr  n’-ossurc  recovery  is  ILkoly  to  fall  during  subcrltical  opora- 
s  a  rc.sult  of  (li  .Slurbanco  of  the  Inlet  Shock  pattern  arising 
nvi  bo'undnrv  layer  separation  as  a  result  of  increased  inlet 
rc,  and  fi'i.m  forward  motion  of  the  normal  shock  which  divides 
"*  suipcrsonic  and  subsonic  flow.  This  boundary  layer 
il'in  phi'r.or.ienon  can,  by  clianging  i.lie  boundary  layer  thickness., 
Ill"  i  liective  shape  of  the  i'n.let  and  thus  lead  to  a  clo.r.gr  in 
■  ick  'Wave  p.iitern,  tlieretiy  reducing  recovery.  This  effect  is 


BEST  AVAILABLE  COPY 


Humjt  t  lent  Pert 


Fig.  28  Fastax  Motion  Pictute  Soqutnct  of  Flow  Inslobilily  nt  Inlok*  of 

1.  7-Scolo  Modol  Romi#f,  Showing  On#-Holf  Cycle  of  the  Otcillolion 

Th»  timing  fight  it  flothing  ol  */  i20  ttcond  inttrvolt;  framing  rota 
of  1880  por  tocond  ot  M  -  272. 


Kt’t.*  r;' ;  1  ••  'i'  th;  !  hvA’  RsyiUilds  iiuabor  and  can  bP  minimized 

11)  ■' 'I'.ucU  by  boundai  y  layer  rontrnl  .  The  adverse  effect  of 

1  hr  ';i>(r>,Kl  ph‘!U)menen,  i.e.,  forward  motion  of  the  normal  shock,  can 
„o!U’r.,  lly  bo  -voided  by  posit  loninK  the  oblique  shock  wave  slightly 
ihead  of  th^-  diffuser  cowl  inlet  so  that  the  normal  shock  can  hn  ex- 
p-  lied  sli^;htly  without  disturbing  the  oblique  shock.  In  this  way  the 
fimall  .mount  of  supercritical  spillover  (which  occasions  a  sliKbt  In- 
crca.se  in  inlet  drag)  permits  stable  operation  of  the  dlffusei  av  its 
most  efficient  point.  The  engine  may  be  operated  with  the  diffuser  at 
the  critical  point,  and  operation  Is  very  likely  to  be  stable  in  the 
presence  of  low  level  pressure  transients  arising  from  he  combustion 
process  for  a  limited  amount  of  spillover.  The  measures  mentioned 
above  are  applied  to  diffusers  as  required  and  are  not  a  part  of  this 
discu.ssion — which  presupposes  stable  diffuser  and  combustor  operation. 

The  mass  flow  passing  through  the  area  Aq  in  the  undisturbed 
freestream  passes  into  the  diffuser  inlet,  is  decelerated  by  a  combi¬ 
nation  of  supersonic  and  subsonic  compression,  and  flows  out  through 
the  area  A^.  This  stream  has  done  no  work;  consequently  Its  enthalpy 
(considered  as  the  sum  of  its  kinetic,  potential,  and  Internal  thermal 
energies)  h;»B  not  changed,  and  the  stream  remains  at  a  constant  total 
temperature  T^  throughout  the  flow.  Since  the  compression  of  the 

high  velocity,  low  density  inlet  flow  to  a  high  density,  low  velocity 
state  is  not  perfectly  efficient,  the  entropy  will  have  Increased; 
this  increase  is  manifest  in  the  reduction  of  diffuser  exit  total  pres¬ 
sure  to  below  the  freestream  total  pressure.  The  condition  of  contin¬ 
uity  of  the  mass  flow  enable-^  one  to  relate  the  freestream  capture 
area  Aq  to  the  exit  area  fay  means  of  the  freestream  and  exit  Mach 
numbers.  Corresponding  to  the  continuity  eq.’atlon 


Pq  ^0  "  Pa  ^d  ^d  ’ 


there  exists  the  following  relationship  between  areas,  pressures,  Mach 
numbers  ,and  temperatures; 


(94) 
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Tbe  capture  area  Aq  may  now  be  expressed  as  a  function  of  the  ratios 
(Aq/Aj)  and  (A^/A^): 


Here  A^^  Is  tbe  diffuser  inlet  area,  and  (Aq/Aj^)  Is  thus  the  effective 
inlet  capture  area  ratio.  Tnua,  if  one  'ixes  the  diffuser  exit  and  the 
freestreas  Mach  numbers,  the  air  flow  capture  ratio  will  be  determined 
by  the  pressure  recovery  for  subsonic  flow  within  tbe  diffuser.  When 
the  value  of  (A^/Aj^)  attains  its  maximum  (critical  flow)  value  in  Eq. 
(96>,  it  is  then  no  longer  possible  tor  the  flow  to  Increase  further. 

If  is  increased  for  a  given  Mq  a*^  constant  air  flow,  the  pressure 
recovery  ratio  (P.  /P.  )  at  this  moment  must  become  a  function  of 

since  it  is  the  only  free  variable  for  a  fixed  value  of  (A^/A^).  Noa 
the  flow  within  what  is  normally  the  subsonic  diffuser  becomes  par¬ 
tially  supersonic  and  experlencca  a  total  pressure  loss  in  returning 
to  Kubsonlc  flow,  inis  acded  diffuser  total  pressure  loss  serves  to 

render  (P  /P  )  variable  iind  consistent  with  (A^/A, and  the  exit 
Xq  01  max 

to  inlet  diffuser  area  ratio  (A^/A^)  will  show  tbe  same  pressure  re¬ 
covery  in  this  region  of  ’swallowed  shock’  flow  regardless  of  geomet¬ 
rical  design.  The  Inlet  drag  (force  over  tbe  aiffvser  i..lat)  >  needed 
as  well  as  the  '-ctal  pressure  recovery  for  the  performance  of  a  diffuse: 


■  !•  •  IV  III  ■!  Ml. !  1 1  1<  r-  .11-  .  .  ills.'u'.  ly  ti-rnlriPil  by  oxperi- 

;■  !  I  ■  v'l  1-  JiilLi,"  r  u<  ''iyn,  altboufih  they  may  be  calculated  ex- 

{  1'  ini'  :  tl'iv  wholly  supersonic,  end  approximately  for  mixed 
-•.ipt  rs.  Ill  1,  -suijson  ic  t  I  o  v 

,^M,.puLln^;  the  total  pleasure  i.  every,  losses  attonjuut  on 
!  1(  .  pv  sapo  ihiuiiph  a  shock  may  rertuny  (le  caiCuLaled  from  the  stream 
M.ich  number  hy  supersonic  flow  t  oory  Subs.onic  losses  are  more  diffl- 
cuii  t-'  -alcu  late  since  1  Iv-w  propel  T  ie«  are  not  simply  related  to  local 

ihfluscr  duct  contours.  Inlet  drag  may  be  .calculated  from  the  magnl- 

•iiiio  ■  f  the  area  intercepted  by  the  cowl  lip  if  the  peripheral  stream- 
li  to  the  cowl  lip  Is  In  the  undisturbed  free  supersonic  stream, 
since  for  this  case  there  are  constant  momentum  and  pressure  forces 
actlnt:  within  the  stream  tube  and  the  inlet  stream  thrust  may  be  com¬ 
puted  from  the  freestruam  flow  state: 

=  Pq  Pa  ^1  '^O  "  ^0  ^1^^  (absolute-  basis).  (97) 

This  may  be  cori'a  .  ted  to  coefficient  form  by  dividing  by  the  product 
.f  tlio  frees tream  dynamic  pressure  and  a  reference  area  : 

A,  (1  +  y  wh 

C  =  -T - 7 — -  (absolute  basis),  (38) 

or.  aj  expressed  on  a  gauge  basis  by  subtracting  the  quotient  o.t 
Aj)  and  (q^^  A^),  an-d  substituting  q/P  =  (1,'2)  y  , 

C,^  -  2(Aj/.*^)  (gauge  basis).  (99) 

'  i 

-iiiei  natively ,  in  ferm“  of  weight  flow, 

-  2,382  y'f^~  0(M)g  =  2.382  A^  w  0(M)^  Lb 

(absolute  b.-isis),  (lOrt) 
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and 
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(absolute  basis).  viOl) 


The  strea.ir  thrust  of  the  air  entering  the  diffuser  is  alLeied  if  the 
diffuser  ir.let  spills  over  and  intercepts  a  stream  tube  of  smaller 
eross-sec t  ion  than  the  inlet  area  In  this  case  the  streamline 

leading  up  to  the  cowl  lip  is  not  parallel  but  is  divergent  by  an 
.Twount  that  Indicates  the  increase  in  stream  area  reodccT  tc  diffuoc 
the  air  to  the  Mach  number  then  existing  at  the  inlet.  The  inlet  force 
is  the-  the  r-troam  thrust  originally  associated  with  the  flow  in  pass¬ 
ing  through  area  in  its  freestream  state  plus  the  added  force 
applied  over  the  projected  area  (Aj^  -  A^),  i.e., 


(102) 


The  diffusion  process  between  Stations  0  and  1  may  thus  be  either  sub¬ 
sonic  or  supersonic  depending  on  the  manner  in  which  the  flow  Is  de¬ 
flected.  If  this  is  initiated  by  a  center  body  protruding  froir  the 
diffuser  inlet,  the  flow  deflection  may  be  gentle  enough  for  it  to  re¬ 
main  supersonic  (or  with  supercritical  spillover),  whereas  if  the  de¬ 
flection  is  such  as  to  correspond  to  a  low  subsonic  flow  kach  number 
at  the  diffuser  inlet,  flow  deceleration  takes  place  abruptly  through 
a  normal  shock.  Equation  (102)  applies  in  either  instance,  but  the 
pressure  relationship  between  0  and  1  will  of  course  be  different  in 
the  two  cases . 

The  inlet  air  flow  capture  and  the  Inlet  drag  can  be  determined 
conveniently  for  the  supersonic  spillover  cone  flow  case  by  a  method 
presented  by  C.  Sears  (of  the  Marquardt  Aircraft  Corporation)  and  W.  J. 
Oriin  (cf  ARO,  Inc.)^.  The  diffuser  inlet  air  flow  and  the  inlet  drag 
may  he  Gained  by  Integrating  the  flow  and  the  force  over  the  conical 

4  Unpublished  paper. 
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surfaco  defined  bv  in*  vertex  of  tte  cone  and  the  Inlet  lip  of  the 
cowl.  Using  the  angle  u  (Fig.  291  subtended  by  this  roninil  siir- 
tace,  the  freestream  Uacii  nuiubei  and  the  cone  angle,  and  making  use  of 


(n 


Fig.  29  Supercritical  Spillover  and  Additive  Drag 

Zdenek's  cone  flow  tables  (Ref.  13),  one  may  obtain  the  flow  Macb  num'- 
her  and  the  flow  angle  0  at  this  surface.  Ope  may  ?.lso  determine 
the  total  pressure  ratio  (P.  /P.  )  across  the  shock.  Since  flow  prop- 

erties  ^.re  constant  along  the  conical  surfaces  referred  to,  one  may 


d4 


simply  m,.! '  l''>  V  t  h.-  tt.tal  area  by  llio  normal  compoii<»nt  of  the 

flow  to  obtain  the  maos  flo*.  The  momentum  component  of  the  Inlet 
dra^  force  is  similarly  obtained  by  multiplying  the  momentum  per  unit 
area  of  the  Incoming  stream  (which  passes  through  the  conical  surface) 
by  the  surface  area  and  the  cosine  cf  the  angle  between  their  respec¬ 
tive  directions,  •'"H  then  nrolecting  this  force  in  the  direction  of 
diffuser  axis.  The  total  inlet  drag  force  on  an  absolute  basis 
consists  of  this  force  plus  the  product  of  the  inlet  area  into  the 
static  pressure  at  the  conical  reference  surface.  The  resulting 
equations  are : 


(P/P^) 


(cos  8)(1 


tan  Ox 


=  A  P  LX\t- 
0\  tn  t 


-]  ^  • 
Vr  L 


M^(cos^  o)(i  - 

tan  u 


The  equation  for  the  inlet  air  flow  may  be  equated  to  the  correspond¬ 
ing  relation  in  the  undisturbed  flow  to  obtain  the  value  for  (Aq/Aj), 
namely , 


Aq/Ai 


,  \  [(P/P,)  8]_^ 


-(cos  0)  (1  -  ---  -) 


Equation  (104)  may  be  expressed  as  an  inlet  force  coefficient  on  an 
absolute  basis: 


’0 

I  /A )  f^t  \  r  1 


(absolute  basis). 
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t  t  P* •  r  f  mi  ma 


Op.  a  baj'Jrt. 


(gauge  basis).  (107) 


The  force  exerted  by  the  air  streaa!  upon  the  diffuser  inlet  may  be 

vipriri  as  the  total  moEuintum  (stream  thrust)  originally  possessed 
by  the  freestream  plus  the  forces  exerted  upon  it  by  the  surrounding 
air  siream  between  Stations  0  and  1.  Looked  at  in  this  fashion,  the 
inlet  drag  ''onslsts  of  two  terms,  one  of  which  is  due  to  the  initial 
momentum  of  the  captured  air,  the  other  being  the  additive  drag  term; 
this  latter  reflects  the  <ii-ag  caused  by  the  Increase  in  the  momentum 
of  the  gas  between  its  free  state  at  Station  0  and  its  captured  state 
at  Station  1.  This  is  due  to  the  force  exerted  along  the  streamline 
from  0  to  1.  The  additive  drag  lor  the  case  of  conical  flow  with 
supercritical  spillover  la  obtained  by  subtracting  2(qQ  Ajj)/(qjj  A^) 
from  the  right-hand  side  of  Eq.  (107),  yielding 


(108) 


Figures  30  through  35  present  solutions  of  Eqs.  (105)  and  (108)  for 

cone  half-angles  of  20°,  25°,  and  30°,  and  oblique  shock-on-rim  Mach 

numbers  of  about  1.9  through  3.0,  for  freestream  Mach  numbers  ranging 

from  low  values  (corresponding  to  sonic  velocity  on. the  cone  surface) 

*0 

up  to  3.0.  Data  are  plotted  in  terms  of  (C^  so  that  the 

‘^add  *n 

drag  is  referred  to  the  inlet  area  rather  than  to  an  arbitrary  engine 
area . 

The  calculation  of  air  flow  capture  and  additive  drag  or  an  in¬ 
let  force  coefficient  is  much  more  complicated  for  mixed  supersonic 
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and  subsonic  inlet  'low.  Approxl&ate  methods  which  are  appllc<ible  to 
conical  dlffusjers  are  presented  lii  Ref,  14.  In  computlnR  Inlet  addi¬ 
tive  ciras  for  oiher  cases,  it  is  iiecussai y  to  use  an  approxiaf.te  vaHw 
tor  lue  static  pressure  along  streamlines  between  0  to  1  and  then  to 
Integrate  this  between  corresponding  areas.  It  Is  clear  thai  if  the 
frccstream  capture  area  is  equal  to  the  inlet  area  A^,  the  additive 

drag  will  be  zero,  and  the  Inlet  force  cooff iciest  m^y  be  ccap’.*ad 
simply  from  A^ ,  A^^  and  the  freestream  Mach  number.  Data  secured  in 
experimental  tests  on  a  Ferri-type  25°  half-angle  conical  diffuser  are 
shown  in  Figs.  36,  37  and  38  as  examples  of  typical  diffuser  operation. 
Figure  36  r- esents  the  relation  hetweeu  diffuser  capture  area  ratio 
and  tota:..  r.  .sotire  recovery.  The  inlet  force  coefficient  (on  a  gauge 
basis)  is  presented  in  Fig.  37  and  a  correlation  of  exit  Mach  number 
with  air  flow  capture  is  shown  In  Fig.  38.  The  inlet  behavior  can 
be  evaluated  at  a  fixed  combustor  heat  release  by  matching  values  of 
at  the  diffuser  exit  and  at  the  combustion  chamber  Inlet,  and  then 
reading  off  the  air  flow  capture  ratio  from  Fig.  38  at  the  desired 
flight  Mach  number.  These  values  of  air  flow  capture  and  flight  Mach 
number  may  then  be  inserted  in  Fig.  37  to  obtain  the  diffuser  inlet 
force  coefficient.  The  latter  can  then  be  used  together  W7.tb  an  exit 
coefficient  obtained  from  combustor  performance  curves  '.o  determine 
the  engine  net  thrust  coefficient.  Similar  curves  apply  to  more  com¬ 
plex  diffusers,  bearing  In  mind  that  performance  is  then  laore  sensi¬ 
tive  to  air  flow  capture  ratio  and  freestream  Mach  number. 

The  use  of  such  graphical  correlations  has  an  advantage  over 
point -by-polnt  algebraic  calculation  of  performance  if  the  diffusers 
are  to  be  operated  over  a  wide  range  of  flight  Mach  numbers  and  air 
flow  capture  ratios.  If  the  engine  is  to  be  operated  at  only  one  con¬ 
dition  of  flight  Mach  number  and  capture  ratio,  the  algebraic  calcula¬ 
tion  of  inlet  behavior  will  probably  be  more  convenient  than  the  for- 
awr  procedure. 
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0  7  oa  ot  ).o 

>ni«t  Air  Coptur«  Rotio  (Aq/A,) 


Fig.  38  T/picol  Corfeiotion  of  Diffuser  Exit  Moch  Number  and  Inlet 
Air  Capture  Rotio 


4.J  The  Combustor 

A  useful  understanding  of  the  combustion  of  fuel-air  mixtures 
in  a  ramjet  engine  may  be  gained  by  breaking  down  the  flow  processes 
into  the  following  steps: 

1.  The  compressed  air  flow  enters  the  combustion  chamber  from 
the  diffuser  at  subsonic  speed  at  a  total  temperature  equal 
to  the  freestreara  value  and  at  a  rate  of  flow  equal  to  that 

at  the  diffuser  inlet  less  that  diverted  to  power  auxiliaries. 

2.  Fuel  is  added  to  the  air  stream,  atomized,  mixed  with  the 
air,  and  in  part  evaporated. 

3.  The  St -ear  v»locity  is  further  decreaserl  to  the  low  value 
suitable  for  the  efficient  combustion  of  the  mixture. 

4.  The  mixture  is  passed  through  a  flame  stabilization  device 
to  initiate  combustion. 


I 

I 

i 

I 

I 

I 

f 

1 

j 
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C- .  Tlie  nlvturo  Is  burned  In  a  cylindrical  coir.buatiuu  chaaber 
ac'-ompaniec  ny  the  followiin;  pbenouena  : 

n.  Transformat ion  of  the  chemical  enthalpy  of  the  fuel 
into  thermal  enthalpy  of  the  gas,  as  narked  by  a 
rise  ill  total  temperature  (or  increase  of  S^^); 

b.  Reduction  in  gas  density  caused  by  the  Increase  of 

temperature,  with  accompanying  acceleration  of  the 
flow;  and 

c.  Reduction  in  total  pressure  of  the  hot  gas  as  a  result 

of  sudden  acceleration  on  being  heated. 

6.  The  hot  burned  gas  is  passed  into  an  exit  nozzle  to  trans¬ 
form  its  thermal  and  pressure  energy  into  kinetic  energy, 
thereby  developing  thrust  by  Jst  reaction. 

The  processes  listed  under  Steps  1,  2,  3,  4,  5b  and  5c  above 
are  amenable  to  solution  by  the  gas  dynamic  equations,  provided  the 
flow  can  be  treated  as  one-dimensional.  The  amount  of  heat  addition 
in  Step  5a  is  a  quantity  which  depends  on  the  chemical  reaction  rate 
of  fuel  and  air  at  the  particular  flow  condition.  It  is  convenient 
to  compute  the  effect  of  the  change  from  chemical  to  thermal  enthalpy 
In  terms  of  the  tc*r.p€r  ature  function  (S^^/  V^)>  since  this  latter  ex¬ 
pression  Includes  such  parameters  of  state  as  y,  R,  and  . 

Once  the  conbustor  flow  conditions  and  the  amount  of  heat  re¬ 
lease  are  set,  c.'nsideratlor.  can  be  given  to  determining  the  space 
required  to  accomplish  this  heat  release.  Flow  velocities  in  the 
comb’ s tor  are  within  the  range  of  one  to  sev'.,ral  hu>  dred  feet  per 
second,  whereas  normal  flame  velocities  of  the  commoner  fuels  are  only 
a  few  feet  per  second,  so  that  it  is  not  possible  to  effect  significantly 
complete  combustion  of  the  fuel-air  mixture  by  means  of  a  flan,-''  spi-ead- 
ing  from  a  single  point  source  within  the  length  of  a  practical  com¬ 
bustion  chamber.  An  Ignitor  I5;  therefore  provided  to  initiate  the  com¬ 
bustion  of  the  mixture  at  a  ''ufficient  number  of  points  so  that  by 
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1  I  ■ 


■  .);isi,inn.  U'.(  whon-  of  the 


i-  ,i  ]  -  .  !  ■  i:..'  It  ‘  ,  ..-pj.  v-^o- Id  into  thi‘  fuel-air  strc-rtffl  to 
ni-  ii!'  .  ,,  i.  !•  ;;.i  ,  t  ,  1  i-  *  t.ike  ill  Mujeh  .1  flam?  tan  stabilize  as  a 
•  :  I  I  I ,  I.  ,  ill  1..^  1  liteit  l)v  r'--c  1  rcu  1  a  t  iiig  combustion  pro- 

!U(  i  .  Mr  ; ii.ssiliif  10  stabilize  a  flame  aevotlvnamically  b;  the  use 
■'!  I  ipi' .>>  i  I  ,  itas  I,,  lot  tills  pi-offiliiri'  lias  not  yel  been  useJ  in  full 

srale  eiiitir.'  s.?  Such  a  st;:b»li7<.i  flame  i{.nitcs  the  combustible  gases 
tl  lAirt;  .  I'ounil  Use  haffi-.-  n.tl  rapiJly  effects  comDusiion  by  virme  nt 
tile  l.ui;t  value  of  the  ratio  of  ignited  sftar.i  perimeter  to  unbui  ned 
I  low  area,  in  1  m.iniuu’  analogous  to  that  of  a  Meker  burner.  In  a  com- 
liustion  .h.iiiiboi'  oi  i  nsuf  f  ie  lent  length,  the  flame  cannot  thus  consume 
all  the  ri.xtut'i'  within  the  contines  Oj.  the  chamber  and  combustion  will 
.u  eordinttly  lie  incomplete.  A  second  type  of  Dame  stabilizer  commonly 
used,  the  can  lomhustor ,  consists  of  a  perforated  tubular-shaped  struc¬ 
ture  closed  ai  tie  upsire.uu  end  io  provide  a  quiescent  piloting  region. 
The  '-ombust  ib  U  .s  arc  ignited  stagewise  through  the  perforated  holes, 

A  correlatioi.  of  sor’c  of  the  parameters  affecting  can  combustor  design 
has  been  made  by  Zclinski  (Ref.  15).  Since  combustion  chamber  flow  is 
everywhere  subsonic,  conditions  downslream  at  the  nozzle  entrance  may 
lie  propagated  upstreum  to  the  diffuser  after  being  modified  by  heat 
addition  and  drag.  Combustion  and  flow  phenomena  are  thus  coupled, 
and  under  some  circumstances  such  a  system  can  become  unstable,  lead¬ 
ing  to  burning  under  conditions  of  violent  pressure  and  flow  oscilla¬ 
tion.  One  result  oi  this  interrelation  is  that  the  combustor  inlet 
flow  Mach  number  is  dependent  upon  combustion  efficiency  and  vice 
versa.  For  purposes  of  analysis,  however,  Ihis  coupling  effect  may 
be  ignored,  and  analyses  performed  merely  on  the  basis  of  an  assumed 
enthalpy  change  in  the  combustion  chamber;  this  enthalpy  value  may  be 
fstimated  by  considering  the  combustion  of  as  much  fuel  as  is  required 
at  a  given  combustion  efficiency  to  obtain  the  required  heat  release. 
Combustion  efficiency  is  thus  introduced  into  engine  performance  cal¬ 
culations  after  the  gas  dynamic  force  balances  across  the  engire  have 
been  completed. 


Calculation  of  the  flow  conditions  between  the  exit  and  the  in¬ 
let  of  the  combustion  ch.amber  now  proceeds  as  outlined  in  Section  2, 
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r^jr.bu.-i 1 1  :)ii  1'  h'-ii.tIiv  t  •  Wf'ii  t-.  .1  (  '-'jHSt.int  siSE  flew  and  cross- 

1  iciiai  aif.i  ,  In  i.i  ocr  Lo  permit  of  the  initial  and  final  stream 
•..'.vs  b-  -  *,  .  on  the  basis  ot  energy  change  at  constant 

stream  thrust.  Ihe  energy  change  Is  described  In  terms  of  the  in¬ 
crease  of  S„  flora  .an  Initl.al  value  related  to  the  Inlet  total  tempera- 
ture  to  .a  fin.il  value  ci.irresponding  to  the  required  hii.tt  release. 

By  Eq .  (67}  the  burned  gas  stream  thrust  Is 


'b  =  *a  Sa  0(‘‘)b- 


(109) 


the  subscript  (n)  now  being  dropped. 
But 


**2  ^2  ^2 


whence 


'h  ’  ^2  ^2  *2 


V% 


0(M)v 


(110) 


The  exit  stream  thrust  of  a  parallel-walled  combustion  chamber 
must  equal  the  Inlet  stream  thrust  less  the  combustor  drag,  since  wall 
force-s  (except  friction)  have  no  component  in  the  direction  of  flow. 
The  inlet  stream  thrust  and  the  drag  may  be  written  respectively  as 


Agd  y  M‘)  , 


(111) 


and 


^2  *2  ^D. 


y  Mf 


(112) 


whence 


=  -^2  -  °b  > 
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or 


o 

0(M).  P,  A„  8, 
y  T  u  i  A  „ 

'■2 


Pj  A,  1 


(C  -  2)  y  m2 

b 

2 


(113) 


Collecting  together  the  functions  Involving  snd  ellHlnatlng  coamon 
ternvs ,  one  has 


1  - 


(114) 


Unfortunately  It  Is  not  possible  to  solve  this  latter  relation  explic¬ 
itly  for  either  M^  or  M^,  even  for  a  fixed  value  of  the  drag,  so  that 
It  Is  generally  necessary  to  fix  M^^  at  a  value  corresponding  to  the 
•selected  exit  nozzle  area  ratio,  and  then  to  use  dJlferent  values  of 

In  succession  to  find  corresponding  values  vii  ^  K )  for  a 

2 

fixed  geometry,  variable  flow  engine.  Otherwise  may  be  fixed  and 

various  successive  values  of  M^^  chosen  so  that  correspondirif;  values  of 

(®a/  V Tt  )  ®»y  *’6  found  for  a  fixed  flow,  variable  geometry  engine. 

2 


In  the  case  of  a  practical  engine,  it  Is  more  convenient  to  determine 

the  relation  between  (S^^/  VX^  ) ,  Mj^  and  Mg  by  experiment.  At  fl.TSt 

2 

sight  It  may  appear  ;.aat  the  correlation  between  Mg,  M^^,  and  (S^/VT^  ) 

2 


Is  a  variable  one  by  virtue  of  the  change  of  Cp  with  operating  condl- 

b 

tlons.  In  many  practical  cases  will  Indeed  be  found  to  alter,  but 

b 

t’lls  change  can  be  expressed  as  a  function  of  Mg,  from  which  It  may  be 
seen  that  the  left-hand  side  of  Eq.  (114)  will  still  be  solely  a  func¬ 
tion  of  Mg,  whereas  the  right-hand  side  will  remain  a  function  of 

)  and  “b""*  V4..K  are  already  fixed  by  the  engine  exit  nozzle 

geometry?  Thus  even  where  Cj.  is  variable,  and  £q.  (114)  is  incapable 

% 

of  being  solved  explicitly,  a  unique  experimental  relation  may  be 
established  between  (S„/ )  and  M.,  for  a  given  exit  nozzle  size. 

ft  Vft  £ 
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If  this  relation  Is  plotted,  values  may  be  taken  from  the  curve  for 
u<?c  In  subsequent  calcu  ■  atlons ,  and  the  curve  may  also  be  adjusted  for 
use  with  other  sizes  of  exit  nozzles.  As  the  product  of  (S  /  }/’t.  ) 

*  tj 


and  appears  on  the  right-hand  slue  of  Eq.  (114),  a  new  value  of 

(S^  /  )  may  be  found  such  that  the  new  0(1!).  ,  corresponding  to  the 

£l  V  ^  D 

new  noziile  size  times  the  new  (S  /  v'T.  ) ,  equals  the  old  product  of 

(S  /  v'T*  )  and  0 (!<>.,  and  therefore  corresponds  to  the  same  In 

a  '  t2  b  z 

this  manner  one  may  form  a  new  correlation  between  values  of  (S  /  VfT ) 

a  tj 


and  for  a  different  exit  nozzle  area.  Within  the  limits  of  experi¬ 
mental  error,  data  currently  available  support  such  correlations. 


The  relations  IXsLed  Immediately  abov  have  now  defined  condi¬ 
tions  upstream  of  the  exit  from  Step  6  back  tc.  Step  4  (see  page  93), 
but  the  flow  at  the  entrance  to  the  combustor  still  needs  to  be  related 
to  that  at  the  diffuser  exit,  with  some  attention  paid  to  the  effect  of 
fuel  addition.  If  the  latter  Is  Ignored,  one  may  equate  the  flow  at 
Station  <j  to  that  of  Station  2  on  a  basis  corresponding  to  continuity 
at  the  same  total  temperature : 


A 


2 


(115) 


may  be  related  to  II2  soon  as  values  are  assigned  to  A2  and  A^, 

-nd  the  relation  of  P  to  P  known.  This  ratio,  (P.  /P  ),  expresses 

^2  ^d  ^2  ^d 

the  subsonic  diffuser  efficiency  In  the  form  of  the  total  pressure  re¬ 
covery  from  Stations  d  to  2.  !f  luel  Is  added  either  contra-stream  or 
co-stream  at  less  than  stream  velocity,  the  stream  total  pressure  will 
be  reduced  slightly,  as  the  force  needed  to  accelerate  the  fuel  must 
be  balanced  by  a  pressure  gradient  set  up  In  the  air  stream.  (The 
magnitude  of  this  total  pressure  loss  Is  small,  amounting  to  only  0.7 
per  cent  for  a  volatile  fuel  Injected  cross-stream  r.t  a  fuel  to  ..ir 
ratio  of  0,06  Into  an  air  stream  whose  Kacb  number  is  0.4.) 


y  / 


itif  M.,  VI  r^uf.  i;urre lat  lull  Is  co-  r.cd  with  that  Of 
vi.'isus  M,  nnd  t  hi-  t  (’mpo’aturi-  paramoter  <S  y'V  ),  a  single  relation 

V  Sir* 


will  ,'.0!i  I'xt.st  between  liie  u^juibui  I  chizber  Ir.lct  wOr"!  *  * '"tc  inti 
the  exit  coidltions  In  terms  of  (S  /  I'T .  ) .  If  required,  the  total 

pressure  ratio  across  tne  couibusioi,  (V  /1‘  ),  ma/  be  obtaineu  from 

M  ,  M,. ,  and  (S  1 ).  The  exit  stream  thrust  ‘t,  may  also  be  ex- 
d  b  a  '  t  D 

pressed  either  ’  n  te^'iss  of  total  pressure,  a  Uach  numher  funutluii 

(f  P  )  ,  aid  area,  or  as  a  product  of  aaBS  '"low,  air  specific  Impulse, 
t  n 

and  the  Marh  nujiib..r  tuiiction  0(M>.  The  air  flow  Is  in  turn  expressible 
in  terms  of  inlet  total  pressure,  inlet  total  temperature,  inlet  area, 
and  the  Mach  number  functior  ((P/P^)Sl,  thus: 

A  P  tf/P  )  --  P  (S  /  KT")  (P/P.)  Sl  A  0(M)  ,  (116) 

^b  ^  b  '-d  ^  ^2  L  ^  Jd  ^ 


whence 


(Pt  /Pt  ) 


(A^/A^)  |(P/P,)  V'T;^) 


(£/P.  )  “  ^b  ’ 

b 


(S^/  Vf^)(A^j/A^)  [(P/P^)  «]  (A/A* )^]/(  1.253), 


Values  of  (P  /P  )  may  be  plotted  as  parametric  contours  on  curves  of 
'b  ’^d 

the  earlier  (S^/  )  and  relations  if  ihe  value  o£  is  kept 

fixed. 
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Figures  39  through  41  show  the  characteristic  curves  of  BOite 
typical  can-type  conbuetors.  In  Fig.  39,  Is  plotted  against 
(Sa/  )  0(ll)jj;  here  varies  over  a  considerable  range  (0.45  to  ■> 
mean  valSe  of  about  0.7),  corresponding  to  a  change  of  (S^/  )  from 

about  3.8  to  5.5,  and  of  from  1.269  to  1.347.  These  data  are  ob¬ 

served  to  correlate  onto  what  Is  nearly  a  single  curve,  even  though  the 
nozzle  throat  area  has  been  changed  from  65  to  70  per  cent  of  the  com¬ 
bustion  chamber  area  corresponding  to  the  above  change  of  0(M) ,  the 
combustion  chamber  temperature  level  has  been  more  than  doubled,  and 
the  engine  pressure  level  has  been  varied  by  a  factor  of  nearly  ten. 


Fig.  39  Corrtlalion  of  Combustor  Initt  and  Exit  Mach  Numbor  for  Con  Combustor 
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DiffUMr  Mach  Number, 


Fig.  40  Con  Combustor  Gos  Dynomic  Flow  Correlotion  for  Exit  Nozzit  Throat 
Area  Equol  to  65  Per  Cent  of  Cor.ibustian  Chomber  Area 
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Fig.  41  Con  Cambuslar  Gas  Dynomic  Flow  Corrtiolion  for  Exit  Noxzit  Throat 
Arto  Equal  to  70  Per  Cent  of  Combustion  Chomber  Areo 
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Tiu'  acrviracy  nt  tins  t  orrelat  ion  sliuws  how  well  one-dimensional  theory 
raav  b*'  applied  to  describe  the  phenomena.  For  some  purposes  It  Is  more 
convenient  to  set  down  the  data  so  that  the  diffuser  Mach  number  Is 
plotted  aijainst  the  temperature  parameter  (S  /  l/T  ),  with  nozzle  area 

fixed,  it  IS  also  possible  to  cross-plot  values  of  the  total  piessure 
r.Ttio  ai-ross  the  e,-mbu!3tor  alonjt  lines  of  equivalent  combustor  drag. 
Figures  40  a!>d  41  illustrate  the  behavior  of  the  combustor,  again  with 
nozzle  area  ratios  of  65  and  70  per  cent.  The  smaller  exit  nozzle  will 
be  observed  to  be  associated  with  lower  dlXfuscr  exit  Mach  numbers  for 
a  given  value  of  the  temperature  parameter,  together  with  a  lower  total 
pressure  loss  during  combustion.  The  total  pressure  losses  Indicated 
in  these  figures  are  made  up  of  the  following: 

1.  Friction  and  turbulence  losses  during  subsonic  diffusion 
between  the  plane  where  the  combustion  chamber  Joins  the 
diffuser  and  that  representing  the  full  combustion  chamber 
area , 

2.  Friction  and  drag  losses  caused  by  the  presence  of  flame- 
holders  and  fuel  injectors,  and 

3.  Momentum  losses  caused  by  heat  addition. 

In  the  cases  Illustrated,  about  9  per  cent  of  the  available  total  pres¬ 
sure  is  considered  to  have  been  lost  as  a  result  of  heat  addition,  and 
2  or  3  per  cent  because  of  subsonic  diffusion,  the  rest  (about  12  per 
cent;  resulting  from  flameholder  drag.  Combustion  system  total  pres¬ 
sure  losses  are  generally  lower  if  the  air  stream  is  first  diffused  to 
a  low  Mach  number  and  use  is  then  made  of  a  baffle  combustor  of  low 
drag.  The  experimental  correlation  of  diffuser  exit  Mach  numbers  and 
temperature  parameters  for  a  baffle  combustor  of  average  drag,  with  a 
nozzle  area  ratio  of  70  per  cent,  indicates  a  minimum  total  pressure 
loss  of  only  about  18  per  cent,  Instead  of  the  20  per  cent  minimum  for 
the  earlier  can  burner  results,  with  a  burner  drag  corresponding  to  be¬ 
tween  6  and  7  velocity  heads,  instead  of  to  between  S  and  10.  This 
lower  pressure  loss  Is  of  great  significance  In  engines  required  to 
operate  at  low  flight  Mach  numbers  (subsonic  values  to  about  1.7).  and 


Ram.lct  Component  Performauc 


of  lessor  i  >;nifi<'ancc  up  *o  a  Mach  luuBDur  oi  ationt  9  5.  At  higher 
flit^ht  Mach  numbers  aerodynamic  total  pressure  losses  of  25  to  50  per 
cent  in  the  diffusei  so  Jar  exceed  those  in  the  combustor  that  only  a 
minor  Improvement  In  performance  is  obtair.cd  by  reducing  the  burner 
dra^ . 


Consideration  of  combustor  operation  is  here  limited  to  the 
aspects  affecting  over-aii  euttlnc  design,  as  combustor  operation  Is 
more  fully  covered  in  'iG  370-9  and  10.  Theoretical  tbermocbemlcally- 
predlcted  combuslion  total  temperatures  and  values  of  are  shown  In 
Figs,  42  and  43.  In  the  first  figure,  the  effect  of  dissociation 
(which  has  not  been  taken  Into  account)  will  be  to  bring  about  a  smooth 
transition  between  the  portions  of  the  curves  on  either  side  of  the 
stoichiometric  fuel-to-alr  ratio.  The  only  further  aspect  of  combus¬ 
tion  to  be  covered  here  Is  the  utilization  of  either  the  results  of  a 
given  test  or  of  theoretical  data  on  combustor  operation  In  the  pre¬ 
diction  of  over-all  engine  performance.  Combustion  performance  enters 
Into  over-all  calculations  In  two  ways,  (a)  by  setting  the  fuel  flow 
needed  lu  effect  a  desired  value  of  combustion  temperature  as  defined 
by  S^,  and  (b)  by  Halting  the  maximum  temperature  ratio  obtainable 
from  a  power  plant  through  a  limited  total  beat  release.  The  first 
limitation  affects  the  efficiency  at  less  than  maximum  engine  power 
output;  tbe  second  Is  equivalent  to  limiting  the  maximum  available 
power  (tbruet),  whatever  tbe  fuel  consumption.  Both  theory  and  experi¬ 
ment  show  tbe  combustion  behavior  to  be  a  function  of  flow  velocity, 
Inlet  temperature,  combustion  pressure,  and  equivalence  ratio  for  a 
given  fuel.  This  behavior  varies  for  fuels  of  different  types,  but 
Is  relatively  Insensitive  to  minor  variations  within  a  given  type. 

The  paraffin  hydrocarbons  from  propane  through  decane  exhibit  a  very 
sir  liar  behavior  for  equivalent  mixture  states.  II  is  possible  by 
means  of  cither  tbe  theory  outlined  In  TC  370-9  and  10  or  experiment 
to  obtain  a  correlation  between  combustion  beat  release  (as  expressed 
in  terms  of  a  combustion  specific  Impulse)  end  enul valence  ratios  for 
various  Inlet  temperatures,  velocities,  and  pressures.  The  flow  ve¬ 
locity  or  Mach  number  has  been  shown  to  be  dependent  on  both  tbe  area 
ratio  between  the  combustion  chamber  and  tbe  exit  nozzle,  and  on  the 
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Equivalence  Ratio,  (E.R.) 

Fig.  43  Thtarttical  Air  Specific  Impt/Isc,  S,,  versus  Equivalence  Ratio, 
(E.R.)  for  Various  Inlet  Temperatures,  T| 
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liUio  final  and  initial  t^as  temperature  (or  )  states.  One  can 
thorclcn-  cilmina’c  me  inint  f  iow  v«‘liuMty  a-;  a  vai  table  for  an  engine 
uslni;  .T  fixed  are;,  exit  le^ziilc,  because  this  velocity  has  a  particular 
value  for  each  operatlui;  puiiit  and  only  this  value  need  be  considered. 
The  pressure  and  temperature  of  the  olream  of  given  velocity  are 
uniquely  ncterminea  by  the  quantities  of  mass  flow  and  total  tempera¬ 
ture  which  are  more  easily  observed  and  controlled.  The  final  (and 
most  convenient)  correlations  are  then  effected;  these  arc  based  on  a 
fixed  combustion  chamber  exit  area,  and  relate  combustion  efficiency  to 
equivalence  ratio  at  a  lixt  value  of  total  temperature.  It  Is  possi¬ 
ble  to  draw  auxiliary  scales  superimposed  on  such  a  plot  In  order  to 
relate  (S^/  )to  and  equivalence  ratio  ,  Calculated  or  experi¬ 
mental  data  are  then  plotted  in  terms  of  engine  air  now.  Pressure 
appears  in  the  mass  flow  parameter,  and  the  effect  of  temperature  is 
shown  by  a  comparison  of  separate  figures  plotted  for  various  appro¬ 
priate  values  of  total  temperature. 

Figures  44  through  46  show  such  experimental  tsst  data  for  a 
piloted  can-type  burner  operated  at  Inlet  total  temperatures  of  150°, 
300°  and  450°F  respectively,  for  a  wide  range  of  air  flow  conditions. 
The  required  values  of  (S^^/  y^T^  ) ,  engine  air  flow,  and  total  tsmpera- 

ture  can  be  entered  In  these  curves.  Interpolating  where  necessary  to 
determine  the  equivalence  ratio  needed  to  produce  the  required  combus¬ 
tion  temperature  (measured  as  .  Total  fuel  flow  Is  computed  from 
the  product  of  air  flow  Into  equivalence  ratio  and  from  the  stoichio¬ 
metric  fuel-to-air  ratio.  In  cases  where  either  test  or  theoretical 
data  are  lacking  for  the  purpose  of  correlating  the  equivalence  ratio 
with  an  actual  operating  condition,  a  reasonable  approximation  io 
likely  performance  may  be  secured  by  selecting  the  equivalence  ratio 
required  for  complete  combustion,  rnd  then  Increasing  the  fuel  flow  by 
dividing  by  a  combustion  efficiency  based  upon  similar  engines  operated 
under  comparable  conditions.  To  do  this  It  Is  necessary  to  employ  the 
theoretical  curve  relating  to  equivalence  ratio,  as  presented  In 
Fig.  43, using  a  value  of  determined  from  the  require*'  engine 

value  of  (S  /  yT  )  and  Inlet  total  temperature,  to  obtain  an  equlva- 
a  tj 
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recent  survey  paper  by  Dugger  (Ref.  16)  summarizes  many  of  the  coB' 
bustion  problems  oi  present-day  ramjet  engines. 


Fi9'  Typicol  Combustion  Efficiency  Correlation  for  Ti 


I  I  Til!  K  \  1  t  ]  e 


it.  '  1.M  1  IS  fii  111  r.i  1 1 V  of  the  de  Laval  tyf  C  ,  although  at 

tin.i'-  iiiil\  itic  l  oiivei^ioii l  pari  ol  the  no7.7.1e  need  be  used  ii  the  exit 
1  low  M.u  h  iiunhoi'  is  only  required  to  be  some  or  sJbsoriAe .  It  is 
cssirit  la  1  i  y  a  matehi.ig  device  fitted  to  the  discharge  er.d  of  the  eom- 
Inist  io!i  cli.vn.lier  aad  serves  to  choke  the  latter  by  Lcildiiig  up  pressure 
to  match  the  available  I'if laser  pressure  at  the  flow  eonditions  of  eap- 
liircd  air  flow  and  eiiKiiie  exit  temperature.  The  divergent  side  of  the 
n')./-7.  h  to  expand  ti."'  combustion  gases  to  a  supersonic  velocity 

iind  correspondingly  low  temperaiure  and  pressure,  and  under  such  con- 
dltioi's  the  output  of  the  combustor  is  independent  of  downstream  at¬ 
mospheric  conditions.  In  the  limiting  case,  the  throat  aiea  of  the 
exit  ni'z/le  mav  approach  that  of  the  combustioi.  chamber;  such  an 
arrangei!’  it  is  known  as  "straight-tailed”  or  "unconstricted."  (This 
could  co.i'espond  to  the  setting  up  of  a  thermal  throat.)  FIom'  through 
the  no7.7.1o  is  nearly  isentropic  if  it  is  uncooled.  For  one-dimensional 
flow  the  i.nlot  and  exit  M.ach  numbers  are  simply  determined  by  the.  fol¬ 
lowing  Mach  number  relationships; 

A  A,  - - —  (subsonic  portion),  (110) 

“  (A/a 

and 

A  A  -  - —  (supersonic  portion).  fl'i'd) 

(a/a  )^ 

^  value  has  already  boon  assigned  to  and  thus  to  0(M)j^  tEq.(114)]  , 
and  this  value  of  ,  when  used  in  Eq .  (Ill,'),  will  define  the  are.a  ratio 

ol  the  convei’giiit;  section  of  liio  exit  nozzle. 


Coropniieiit  Porf  jrmj 


'  '  1  •  1  ;  -.ilii!'  t  1  h*'  p.' niii'ii  ;iini  underoxpa  nded 

.  ■■  ■  1:  ^  ;■  i  ■  V  1  M- 1  y  t  iiiu-i  lies  in  a  eumbinaiion 

;■■■■  •  ‘1  .ei.iitiiii  wiilkiii  1  i.e  ramjet  (which  affects 

•...  ;r.  ip-'vii!’.  '1  me  i  x  i  I  .'Uul  flight  altitude  (which 

:  ur*' )  ,  M;!,';.'.lc  sli.ipcs  f''r  parallel  !?;«-- 

ir  '|.:c  1111*  mai.  .n  i-enpino  hi-  varuiiis  riethods,  iiR'luding  the  "method 

■  ;  I  ii.u  ii' I  el  ■  i  n  .  ihe  resulting  dcsigiis  afford  thooretiral  (ideal) 
\  ,1  1  ic.  s  ef  exii  taiu'.t.  apart  Irnni  the  momentum  loss  uue  to  skin  tric- 
1  1  in  ,  whivli  Is  ;  ;;pr  ■>>;  imately  1  per  ■•ent  of  the  ‘otai,  but  the  minimum 
l.  ngtli  oi  .UK'h  an  lOeal  n(.)zzle  tends  to  come  out  cxctssivc  sn  terms  of 
ivailalilc  sjiacc  .  Since  an  ideal  diveigfcnl  design  is  therefore  longer 
than  IS  generally  acceptable,  the  subsonic  convergence  is  often  made 
more  rapid  thai:  the  ideal  This  may  cause  a  separation  oi  flow  at  the 
thiii.it,  leading  to  a  reduced  eoe*ficient  of  discharge  with  ."iOiiunlforra 
f'n*  ill  the  s  inersonlc  s*'  lion.  If  there  were  no  need  for  further  ex¬ 
pansion,  a  sonic  jet  would  he  perfectly  efficient,  but  this  is  modi- 
lioU  by  supersonic  expansion  and  may  be  less  eiiicient  than  the  ideal 
because  of  nnnun t f ormi t y  of  the  entering  flow.  Figure  47  illustrates 
the  theoretical  reduction  oi  eiticicncy  which  occurs  when  a  family  of 
supersonic  tscntropic  nozzles  are  shortened  to  les.s  than  the  lengths 
required  to  produce  parallel  flow  at  the  exit. 


Since,  exit  nozzles  do  not  in  practice  give  rise  to  isentropic 
one-dimensional  flow,  the  coefficient  of  discharge — which,  depends  to 
■n  considerable  extent  on  the  inlet  geometry--is  always  less  than  unity, 
so  that  the  realizable  exit  thrust  will  be  below  the  theoretical  value, 
■file  principal  loss  arises  from  the  existence  of  a  radial  component  in 
the  exit  velocity,  smaller  losses  occurring  as  :<  result  of  nonuniformity 
of  t.ic  exit  velocity  and  as  an  effect  of  skin  frictiO!':  This  drop  of 
actual  exit  jet  thrust  below  the  theoretical  value  is  manifest  in  a 
drop  of  the  static  pressure  acting  on  the  diverging  nozzle  are.i.  The 
etiectivenesa  of  a  nozzle  may  thus  be  determined  by  '’ithcr  a  direct 
experimental  measurement  of  the  thrust  or  of  the  pre.xsure  distribution, 
or  l)y  a  calculation  of  the  internal  flow  field,  either  within  the  noz¬ 
zle  or  across  the  exit  jet  Figures  4h  and  40  Illustrate  the  effect 
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Fig.  48  Actuol  Variation  of  Thrust  with 
Length  of  Conical  Nozzle  of 
Fixed  Area  Rotio 
Doti  taken  from  R«f.  iti. 


Fig.  49  Actual  Voriation  of  Thrust  of 
Conical  Nozzle  with 
Divergence  Angle 
Data  taken  from  Ref.  18. 


on  actual  norzie  thrust  efficiency  of  shortening  the  nozzle  and  vary¬ 
ing  the  divergenct  angle.  Further  data  on  the  efficiency  of  practical 
nozzles  will  be  found  in  TG  370-13. 


The  magnitude  of  the  expansion  area  ratio  (A^/A^)  to  yield  maxi¬ 
mum  thrust  thus  involves  use  of  an  optimum  expansion  angle  which 
enables  a  compromise  to  be  struck  between  excessive  friction  and  ex¬ 
cessive  expansion  losses.  In  an  over-all  force  balance,  the  gross 
forward  engine  thrust  will  be  made  up  of  the  rate  of  flow  of  momentum 
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thi  '.'uuh  )ic  t  s  1 1  ii/./iii-  ihi  iai  together  with  the  Icrward  acting  com- 
ol  the  wall  gauge  pressure  summed  over  the  expanding  area  of 
the  exit  nez/.le.  This  latter  force  attains  a  maximum  when  the  expan¬ 
sion  no/,/, io  area  is  extended  until  the  gauge  pressure  acting  on  the 
added  area  leduees  t;)  zero.  External  drag  considerations,  however, 
modify  the  optimum  size  of  the  exit  nozzle  since,  for  example,  a  gain 
of  thrust  at  the  expense  of  an  equal  increase  of  drag  obviously  pro¬ 
vides  no  not  improvement.  The  external  drag  will  be  a  minim-'r'  if  the 


nozzle  exit  area  Is  the  same  as  the  combustion  chamber,  thereby  avoid¬ 
ing  the  setting  up  of  base  or  wave  drag.  In  this  case,  (A„/A„)  - 
(A^A^),  It  is  impossible  to  satisfy  this  requirement  for  all  noBei- 
ble  flight  altitudes,  unless  a  variable  geometry  nozzle  is  employed, 
so  that  expansion  to  an  arbitrary  area  different  from  A^^  results  in 
the  exit  static  pressure  differing  in  general  from  the  ambient  pres¬ 
sure.  The  simplest,  type  of  variable  geometry  design  is  an  inverted 
nozzle,  having  a  movable  center  bouy,  such  a  design  is  associated  with 
difficulties  in  the  way  both  of  cooling  and  of  rigidity.  If  an  ex¬ 
pansion  to  ambient  pressure  is  required,  the  exit  total  pressure  must 
first  be  obtained  by  following  the  total  pressure  change  through  the 
diffuser  and  the  combustor;  the  ratio  of  the  exit  static  to  the  total 
pressure  than  determines  the  exit  Mach  number,  which  in  turn  fixes  the 
expansion  ratio  from  nozzle  throat  to  exit: 


(P/P.)  =  (Pn/Pf  ) 

'■  e  "  c 

and 

(A^/A^)  =  (A/A*)g  . 


The  thrust  of  the  exit  jet  r’ay  thti.  be  tivnputcQ  Xio’u  either  Eq. 
(39)  or  Eq.  (67),  as  given  in  Section  2,  using  Eq,  (59)  for  expansion 
to  the  ambient  static  pressure,  and  Eq.  (67)  for  expansion  to  the  full 
area  value. 


Ill 
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TiiL'sc  thrust  values  are  expressed  on  an  absolute  basis.  If  they 
reiiuired  on  a  gaa^e  basis,  a  force  (A^  ^0^  ’  product  of 

i.'  .I'lbient  static  pressure  into  the  exit  area,  must  be  deducted  from 
i  e  l  i^at  i-.aiid  Side.  Engine  net  thrust  on  a  gauge  basis  may  be  found 
o  i  a  ditterence  of  the  exit  s.nd  inlet  forces, 

ilic  additional  loss  oi  t-ngine  thrust  ai  ising  from  nonone-diraen- 
.i.rr.lity  oi  thermal  cr,er,,y  (t(>tal  temperature)  in  the  exit  stream  has 
1. 1'.i  discussed  by  Bader  (lief,  .19),  There  are  likely  to  be  regions  of 
.  jyii;.;  fuel-air  latio  or  combii.st  ion  efficiency  in  the  burnt  gas  stream 
lo'rii'.j:  irom  the  corabustoi .  K  stream  which  ha.s  burned  at  ,a  leaner  than 
-o  t  'a;,*-  iuc’l-air  mixture  or  which  has  less  than  average  combustion  ef- 
:  eiencv  will  thus  have  a  lower  value  of  S,  or  of  total  temperature. 
S.i.."ii.lavly,  regions  of  irrcatcr  heat  relenae  will  have  a  higher  or 
t  eiiiperat  ure.  An  approx ima to  calculation  of  this  effect  may  be 
de  by  consjdering  a  set  of  concentric  streams  differing  stepwise 
,  :  (ur.  l  acii  other  in  respect  of  total  temperatur®  or  of  S  ,  The  loss  of 

A 

u.rust  m.ay  easi.'y  be  demonstrated  qualitatively  for  the  case  when  the 
e:;it  jet  is  consideic-d  to  consist  of  a  hot  and  a  cold  stream  both  at 
,s;r,".c'  Mach  number  as  tlie  inlet  flow.  The  exit  gauge  thrust  caused 
i.y  tl',c.  cold  gas  will  be  equal  to  or  less  than  its  inlet  thrust  and  it 
,ill  •nakc  no  net  contribution  to  engine  thrust  other  than  drag  (corre¬ 
sponding  to  diffuser  losses  and  burner  drag).  The  engine  will  then 
opciate  with  less  tnan  the  optimum  air  flow  and  at  a  higher  equiva¬ 
lence  ratio  and  will  accordingly  be  less  efflcie’'t.  The  over-all  fuel 
.specific  impulse  of  the  engine  will  approximate,  but  will  be  lower  than. 


,  .  t  .i-.i  111  iioiiolu'-il  iPK  iis  tonal  i  t  y 

;  .  ii  .  1  :;t  ’  :.'a  I  ';-,  i!  ;  -^imj  It^'  t<i  Lif.it  it  arj  a  icUuctioii 
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.  pip-.- ii'i  M.iili  oumhi'i'  may  bo  iit'*ermined  b>  t  indinn 

I  ■,iim  ,  I  till.'  I  I  I' t  !•  1  !)ii ;  ■ ' 'IIS  I  t  hi  stream  tubes  at  constant  local 
t'l.i;  l  ■•tniM' I'.u  ui'e  .  ■*liiii.  e.iu.iis  i  itt  iiroiiuc  t  ot  the  a'/oraije  total  tem- 
•  ■r.it'.u'i  ni-a  v  i.M'.uii'.  ol  the  tutal  flow: 
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uiiit  lime  m;iy  be  expressed  in  terms  of  Mach  number, 


pre5sur<',  ai'.u  total  temperature  as  follous: 
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A,(i  rx;^) 


T!i'‘  thrust  of  the  aiv  stream  may  bo  obtained  by  summing  the  thrust 
coi’.tr  iljut  ions  oi  tlic  stream  lubes: 


A  0(M)  1-  A.  ‘  r.  m,  (S  )  0(M) 

••i  i  ^  i  ^  “  i  ' 


1.  p  ^(s  V  r, )  0(M)^ 

P  0(M)  8  i-  A  (1/  KtT ) 

1  ’ 


V)ui  (S^  \'^t^  ■’  nade  up  larfcly  of  gas  constants  since 
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This  expression  may  be  simplified  if  Y,  c^,  R,  and  y  are  considered 
constant : 


z  k AI/  fTT) 

i  ^  ^1 


(125) 


If  Eq.  (124)  is  divided  by  VW^,  there  follows 


r  A, 

L  A.  (1/  fT') 
1  ^  ^1 


If  this  latter  equation  is  now  divided  Into  Eq.  (125),  there  results 
the  ratio  of  the  actual  thrust  to  that  produced  by  a  one-dimensional 
stream  of  the  same  energy: 


perfi.’i  rri."r.; 
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(126) 
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f. 


U',.  .;.i\  :  :i<  i.tl.-i  I  Dll  ,1  tempoi'a  lure  profile  when,  ex- 

■;! . '  !  i  he  e\!i  jet  ai't'.i  IS  at  twico  the  mean  temperature 

Hill!,  till  iitl.e!  half  is  at  half  t  lie  mean  temperature,  tnus  satisfying 
1.  ;.  t  !  1' i  )  II.  1  I  i,aril  to  eonstaiK  y  of  heat  release.  E.itci  ing  these 

t  11.;,.  i;,|,  (12-1). 
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Thas  there  is  approximately  a  6  per  cent  loss  of  thrust  in  this  simple 
'■.IS-.  Calculation  of  the  thrust  loss  for  actual  combustion  systems  is 
further  ror.plicatee.  by  dissociation  phenomena  that  occur  during  burn- 
iiiii  at  elevated  temperatures  and  by  changes  in  the  values  of  Y,  Cp, 
anh  R  with  temperature  and  ga:>  composition,  as  indicated  earlier. 

The  equations  presented  so  far  are  useful  only  when  wotal  tem¬ 
perature  distributions  are  known.  If  it  is  desired  to  calculate  the 
effective  thrust  available  from  burning  known  amounts  of  fuel  in  vari¬ 
ous  parts  of  the  air  stream,  the  following  procedure  needs  to  be  fol¬ 
lowed  i.i  order  to  deduce  the  stream  temperatures  and  stream  thrusts 
from  the  local  equivaleiice  ratios.  It  is  assumed  that  the  cembustor 
has  been  designed  on  a  une-dimrnsl ona 1  flow  basis  and  that  all  flow 
areas  and  Mach  numbers  and  the  equivalence  ratio  are  fixed.  The  thrust 
having  been  determined  for  one-dimensional  liow,  it  is  desired  to  cal- 
<  ulale  the  thrust  loss  that  would  result  If  the  hot  burned  gases  were 
not  mixed  with  the  remaining  throughput.  Nominal  one-dimensional 

values  exist  for  T  ,  n  ,  (E.R.),  M  ,  and  o  (the  fraction  of  the 
'•Q  ' 

mroughput  initially  heated);  the  nominal  value  of  S_  corresponds 

1 

to  i7^(E.R.).  The  heat,  however,  is  added  to  only  a  fraction  a  of  the 
flow  and  this  fraction  burns  at  an  equivalence  ratio  of  (E.R.)q/o  . 
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The  ratio  of  exit  thruBf)  nay  now  be  tound  iroa 


■frni  “  *. 

uni  a  • 
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whence 
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This  reduction  of  effective  thrust  can  be  regarded  as  a  reduc¬ 
tion  In  effective  combustion  efficiency  for  a  particular  chosen  engine 
operating  condition.  Figure  50  shows  this  loss  of  effective  coabus- 
ttcn  efficiency  In  a  combustion  system  burning  kerosene  at  an  inlet 
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Fiy.  50  Ptreantagv  Lest  of  Apporent 
Cembustien  Efficiency  Due  to 
Nentnixed  Exit  Jet 


total  temperature  of  400'’f.  The  effective  thrust  loss  is  greatest  for 
engines  operating  at  low  equivalence  ratios  where  only  a  small  fraction 
of  the  over-all  flow,  numerically  equal  to  the  equivalence  ratio,  is 
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■  R  f.i  -  I  ]st  un' .  Thf  rt'sidue  of  the  air  not  involved  In  conbus- 

lioii  ill  !!■.  i  ;r.a i  1  y  '.'v-pass-od  aroupo  the  combustion  space  and  Mixed  *’lth 
the  products  ''f  comlnis  1 1  on ,  the  process  being  slMl!ltir  to  that  occurring 
In  a  turbojet  lombust ion  system.  Ii  may  be  noted  that  this  reduction 
of  comhusiion  cfilciencv  Is  a  factor  to  be  applied  to  a  calcu  .ted  fuel 
specific  impulGC  or  ctr.hustion  efficiency.  For  lean  mixtures,  in  spite 
of  .t  certain  .tmount  of  loss  due  to  nonmlxlng.  the  fuel  specific  i»- 
pulsc'-  usually  optimize  .at  values  of  y  less  than  rtclehiometrl.-  ao 
that  an  optimum  engine  design  will  still  be  similar  to  one  with  uni¬ 
form  c.\it  flow. 

No  compensation  factor  Is  needed  for  nonone-dlnenslonallty  of 
the  exit  stream  when  the  exit  jet  thrust  is  Measured  directly  for  the 
determination  of  combustion  efficiency.  In  the  case  of  experimental 
combustion  efficiency  evaluation  from  thrust,  the  thrust  losses  from 
jet  iior.one-dir..-r.3i..r,alit ics  ef  all  kinds  are  Included  in  the  combus¬ 
tion  efficicnev. 

A  further  problem  in  the  design  of  the  exit  nozzle  results  from 
the  fact  that  the  high  temperature  of  the  gases  leaving  the  combustion 
chamber  exposes  the  nozzle  to  a  considerable  flux  of  heat  from  both  con- 
uuciion  and  radiation,  which  may  cause  Jt  to  fail  structurally.  Consid¬ 
erable  effort  has  therefore  been  expended  in  an  endeavor  to  develop 
suitable  nozzle  materials  to  withstand  the  high  temperatures  of  the 
gases  issuing  from  the  combuotlon  chamber.  For  example,  while  a 
nickel-cobalt  steel  alloy  may  melt  in  the  range  2500-2600°F,  while 
retaining  its-  strength  up  to  2000°F,  the  melting  point  of  molybdenum 
approxim.itely  4750°F,  with  an  uppei-  stress  limit  at  about  3000°r. 
Another  approach  is  to  Increase  the  reflectivity  of  the  subsonic  and 
throat  sections  by  coating  them  with  such  metals  as  gold  or  platinum. 
.Alternatively,  a  combustion  chamber  cooling  shroud  may  be  designed  to 
afford  relief  by  spilling  an  annulus  of  cooler  air  onto  the  nozzle, 
but  this  arrangement  mav  result  in  a  reduction  of  the  effective  throat. 
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V  ramjet  performance  calculations 


ihe  Dt'jof't  of  this  section  is  to  show  how  the  basic  therwo-^y- 
numic  calculations  of  a  ramjet  engine  are  made.  The  various  engine 
arr-c.gcncr.lc  v,hlch  c.an  be  chosen  for  specific  requirements  are  first 
considered,  and  from  these  the  one  having  the  beat  riauiC  of  meiit.  is 
pti'ked.  Some  arrani^ementx  may  be  found  which  will  produce  high  fuel 
economy  only  at  the  design  point,  otuers  will  provide  an  acceptable 
performance  over  a  wide  range  of  speeds  or  power  outputs.  The  main 
requirements  and  Indices  of  merit  for  engine  performance  will  arise 
from  considerations  of  engine  application  presented  In  a  later  chapter. 
Typical  requirements  will  be  stated  here  without  extensive  discussion, 
.IS  there  are  a  number  of  Inevitable  compromises  that  arise  In  the  de¬ 
sign  ot  engines  for  any  purpose.  These  arise  from  requirements  of 
>>s»ge,  from  problems  of  Installation  and  from  the  conflicting  require¬ 
ments  for  high  performance  presented  by  the  various  major  components. 


5.1  Thrust  Performance 

The  useful  oulput  cf  a  ramjet  engine  is  the  net  difference  be¬ 
tween  the  forward-acting  reaction  on  the  engine  set  up  by  the  ejection 
of  the  exhaust  gases  and  the  drag  associated  with  the  air  flow  capture 
at  Inlet.  This  result  Is  expressed  either  as  a  net  thrust  or  as  a 
dimensionless  coefficient  Cj,  ,  the  main  Inlet  and  exit  flow  relation* 
having  been  developed  In  Sections  2  and  3.  Application  of  these  rela¬ 
tions  requires  the  selection  of  values  for  the  two  Independent  gas 
dynamic  variables,  flight  Mac  i  number  and  temperature  parameter 

(S  /-i/t.  ),  In  addition  to  the  need  to  detail  the  engine  geometry, 
a  r  iQ 

These  relations,  togethc*-  with  the  ambient  condllions,  determine  con¬ 
ditions  within  the  combustor,  and  allow  the  selected  temperature  parame¬ 
ter  to  be  translated  Into  an  operating  fucl-to-alj  ratio  by  ae^r.s  of 


r.ir  1 1  r ii i>i  j;.,  Iiurmr  which  lolate  the  combusticr.  effi- 


u  ..I  1  t '  i  h 
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Tile  nt'i’hanics  ,)i  calculation  i-roceeu  in  a  straightforward  man- 
r.(  >'  1  r<iin  the  iiiiii  t  inn  of  the  difluser  and  combustor,  using  theoretical 
icl.itiJns  or  lie'll  •  lulls  such  as  those  si  nwn  in  Figs.  36  to  for 
tl.i  ;i  i  i  1  ti''cr  .  ;uiii  -lo,  •!  1  and  l-l  to  46  for  the  combustor.  A  re- 

l:it  lot. ship  su.  !■  dial  of  Eq .  (114),  or  test  results  zuch  as  those 
shown  It,  Fii;-  dll.  will  yield  the  required  relationship  oetween  luo 
I'tiiosii-  c;;:!  M.u'h  rv'mtior  and  the  temperature  parameter.  One  may 
then  ueicrmino  tlic  air  flrr  capture  by  means  of  Eq.  (96',  or  from  test 
insults  Such  as  those  shown  in  Fig.  38.  The  air  capture  ietermlnes 
the  inlet  drac.  as  itiven  by  Eq .  (98)  or  by  experimental  results,  as 
snown  ill  Fit;.  37.  The  xlt  force  is  determined  from  the  selected 
value  ol  (S^  '  ^  T^. ) ,  the  nir  capture  is  found  for  the  diffuser,  and  the 
exit  noz/.ic  performance  is  then  found  from  Eq.  '120).  The  exit  and 
inlet  forces  may  be  converted  to  a  gauge  basis  by  subtraction  of  the 
ambient  static  pressure  force  on  the  projected  surfaces,  or  expressed 
as  coefficients  by  division  by  the  product  of  a  reference  area  into 
the  dynamic  pressure;  they  are  thus  functions  of  the  flight  Mach  num¬ 
ber,  (S^/  external  pressure  level,  but  normalizing  them 

by  dividing  by  q^  results  in  values  which  are  independent  of  the  pres¬ 
sure  level.  Different  ways  of  plotting  these  three  variables  of  flight 
Mach  number,  thrust  coefficient  and  air  capture  ratio, are  shown  in  Figs 
cl  ?‘'d  c2.  Each  figure  thus  completely  represents  the  gas  dynamic  be¬ 
havior  of  the  engine  under  all  conditions,  but  neither  includes  any  in¬ 
formation  as  to  how  the  exit  gas  temperature  is  to  be  obtained. 


It  is  clear  from  an  inspection  of  Fig.  51  that  the  thrust  co¬ 
efficient  varies  weakly  with  the  freestream  Mach  number  and  strongly 

with  the  value  of  (S  /  -i/T  ).  Ir.  addition,  the  lines  of  constant 
a  r  tQ 

(S,^/ )  relating  to  Cj,  show  three  types  of  engine  inlet  aerody¬ 


namic  behavior.  Above  The  line  (A,,/A, ) 

U  I  mux 


i  the  thrust  rises 


only  slightly  because  of  the  reduction  of  engine  air  flow  due  to  inlet 
spillovei'  with  an  increase  in  (S  /\'T.  ).  Below  this  curve  the 

a  r  to 
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Fig  $1  EHtel  of  Fr*«ttr*om  Moch  Numbtr  ond  Engin*  T*mp*ratui« 
Poram*t*r  on  Thrutt  Coofficitnr 


FriMlrcom  Moch  Number, Mg 

Fig.  32  Effect  of  Freestreain  Mach  Number  and  Engine  Tt.r.psrature 
Parameter  an  Thrust  Gtefficient  and  Airflovr  Captura  Areo 
Ratio 
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(S  )  t:urves  c,{hibit  two  types  of  behavior,  depending  on  whether 

or  not  the  flight  speed  is  above  or  below  the  value  for  which  the  ob¬ 
lique  snocK  fri'm  the  spiki  Just  falls  within  the  inlet.  As  long 

as  tne  obfiqu-;  shock  is  ?  •.  >f  the  inlet,  an  Increase  in  the  Mach 

number  reduces  spillover  between  the  oblique  shock  and  the  inlet  lip, 
so  that  the  air  flc."  captured  increases  more  rapidly  than  does  the 
speed.  This  increase  is  reflected  in  the  rise  of  exit  thrust  or 
thrust  coefficient  with  Mach  number  until  full  air  flow  capture  is 
obtained.  Above  this  speed  no  supersonic  spillover  occurs,  a,  \  the 
inlet  capture  area  becomes  fixed,  leading  to  only  a  linear  ji^ase  in 
the  air  flow  with  increasing  speed.  As  the  dynamic  pressure  Increases 
as  the  square  of  the  speed,  the  net  thrust  coefficient  drops  off  In 
the  ratio  of  the  increase  of  air  flow  with  dynamic  pressure,  or  in¬ 
versely  as  the  Mach  number.  This  simplified  picture  Is  somewhat 
complicated  by  the  relationship  between  Inlet  drag  and  mass  flow  cap¬ 
ture  so  that  as  one  moves  above  the  (An/Ai)___  -  1  line  the  air 

u  ^  max 

flow  and  inlet  force  coefficient  both  fall  off.  For  most  diffusers, 
this  fall  in  the  inlet  force  coefficient  is  enough  to  offset  the  thrust 
reduction  caused  by  lessened  air  capture,  and  a  small,  but  at  times 
useful,  Increase  in  thrust  is  obtained  by  subcrltlcal  dlffusor  opera¬ 
tion.  Figure  52  thus  shows  the  aerodynamic  aspects  of  the  engine  be¬ 
havior  better  than  Fig.  51,  as  it  Is  easier  to  sL  w  the  variation  of 
subsonic  inlet  flow  with  the  temperature  paraneter  (S./^  T*  )  which 

largely  determines  (An/A,)  ;  this  results  in  linen  of  constant  thrust 

coefficient  becoming  lines  sloping  In  general  down  to  the  lower  right 
from  the  upper  left,  with  breakf’  Ip  slope  and  spacing  on  passing  through 
regluac.  of  different  rates  of  air  flow  change.  Figure  51  Is  simpler  to 
Intel pret  from  the  thrust  output  point  of  view,  while  Fig.  52  repre¬ 
sents  the  aerouynaslc  and  thermodynamic  picture  more  directly. 

Pressure  and  temperature  (In  addition  to  the  temperature  parame¬ 
ter)  must  be  fixed  in  order  to  set  the  fuel  flow  or  the  equivalence 
ratio  required  to  produce  a  given  value  of  •  I*  combustion 

were,  always  perfect  with  no  changes  due  to  dlssoclatl'''n ,  one  could  ob¬ 
tain  a  correlation  between  S  aiiU  the  equivalence  ratio  uc  shewn  in 

a 
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Fit^;.  1},  ;.r.  1  v.o'.;Ut  ihus  teed  to  fix  only  T  ,  ,  and  (S  /  i/t.  )  to  ob- 

0  U  a  '  tQ 

tain  the  ru*'l-air  i at io  ne'ded  to  produce  a  given  thrust  coefficient. 
Sineo  ceiTtbustio:.  •:  f  f  if  ie i  ■  f  i"  fact  so.newhat  pressure-sensitive, 
one  must  .ilso  fix  the  ambient  pressure  .  which  thus  enables  the  com- 
buatiin  pressure  to  be  doterminec  from  the  relationship  between  ambi¬ 
ent  pressur'' .  flight  Mach  number,  and  the  diffuser  recovery.  If,  for 
instance,  one  were  to  set  the  ambient  pressure  and  temperature  to  that 
corresponding  to  NACA  Standard  Day  conditions  at  an  altitude  of  30,000 
feet,  one  could  combine  the  relationship  between  YTq  and  the 

erjulvalence  ratio  shown  in  Figs.  44,  45,  and  id  with  Figs.  51  or  52  to 
obtain  the  tlirust  coefficient  in  terms  of  equivalence  ratio.  If  the 

given  aiuc-o  for  (S  /'  )  .  T.  ,  and  the  air  weight  flow  are  entered 

a  V  ty  ty 

into  Fig.s.  M .  45,  or  46,  the  corresponding  equivalence  ratio  may  be 
derived,  as  shown  in  Fig.  53  for  the  given  conditions.  These  equiva¬ 
lence  rati<(  cu» ves  could  also  have  been  plotted  on  Fig.  52  but  for  the 


1 
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resulting:  multiplicity  oi  Grid  lines.  If  the  temperature  level  at  the 

;:ivor.  iii”'"  altituat-  weitf  changed,  engine  thrust  would  be  affected 

because  of  the  effect  of  inlti.  t‘*:fipe’*’ture  on  for  burning  at  a 

fixed  equivalence  ratio,  resulting  from  the  change  of  combustion 

efficiency  and  the  increase  of  T  in  (S  / 4/T  ).  The  general  result 

to  ^  f 

is  that  while  S  increases  slightly  as  the  inlet  total  temperature  Is 
varied,  the  thrust  fall.s  because  of  the  correspondingly  greater  in¬ 
crease  of  the  root  of  the  inlet  total  temperature  in  tue  denominator 
of  the  temperature  parameter  term.  Figure  51  shows  the  over-all  effect 


Fig.  54  Effect  of  Ambient  Temperature  on  Thrust  C  ifficient  of 
Constant  Altitude  of  30,000  Feet  and  Equivalence  Ratio 
of  0.8 


of  ambient  temperature  on  the  engine  described  in  Figs,  51  through  53, 
as  obtained  by  contrasting  the  thrust  coefficient  for  a  constant  equi¬ 
valence  ratio  of  0.8  under  Standard  Hot,  Standard  Cold,  and  Standard 
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Day  conditions  at  30,000  feet.  This  effect  is  large  and  can  be  quite 
serious  fu.  marginal  powc-i-  plants.  When  the  engine  operating  altitude 
is  changed,  two  effects  ensue,  caused  respectively  by  pressure  level 
and  ambient  temperature.  At  low  altitude,  higher  temperatures  gener¬ 
ally  occur  together  with  htgb> r  pressures;  combustion  at  low  altitude 
is  also  generally  more  efficient  because  of  the  high  pressure  level, 
hut  thi;  value  of  (S^/  )  is  lower  because  of  the  greater  propor¬ 


tional  increase  of  over  that  of  which  generally  leads  to  a 

G 

duction  of  the  thrust  '.•oef f iclent  at  eca  level  as  compared  to  the  value 
at  the  base  of  the  stratosphere.  An  increase  of  altitude  above  the 
base  of  the  stratosphere  results  in  lower  combustion  pressures  attended 
by  reduced  values  of  combustion  efficiency  and  S,  without  any  reduction 
in  T  so  that  the  thrust  again  falls.  These  further  effects  are  shown 


in  Fig.  55  for  equivalence  ratios  of  0.3  and  0,8.  Of  note  is  the  rapid 
decrease  in  thrust  coefficient  for  reduced  flight  speed  at  an  altitude 


Fig.  55  Effect  of  Altitude  on  Thrust  Coefficient  for  NACA  Stondord 
Doy,  ot  Constont  Equivolence  Rotio 
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nf  5u,000  feet  and  an  equivalence  ratio  of  0.8,  due  to  the  fact  that 
the  fall  in  combustion  efficiency  is  thus  large  for  a  change  cf  mass 
flow  or  combustion  pressure.  The  change  of  combustion  efficiency  with 
altitude  is  less  pronounced  at  .an  equivalence  ratio  of  0.3  and  the 
slopes  of  tie  corresponding  Cp  curves  are  more  nearly  the  sane. 

One  may  compute  ramjet  engine  air  flows  at  suitable  points,  and 
also  values  of  fuel  flow  from  the  air  flow  and  the  equivalence  ratio, 
as  indicated  by  curves  such  as  those  ot  Fig.  53.  Th=  .tesultlua  rela¬ 
tions  between  fuel  flow,  flight  Mach  number  and  thrust  may  be  plotted 
as  in  Fig.  5b  to  give  the  final  engine  performance  relationships  re- 


Fin.  56 

Thrust  Coefficient  iot  Various 
Fuel  Flows  Wf  at  Constont  Altitude 


quirod.  Associated  with  each  point  on  Fig.  56  are  a  thrust  and  a  fuel 
rate.  From  the  quotient  of  these  parameters  a  figure  of  merit,  1^,, 
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rr-  tijci  Specif  ic  •  mptilse ,  may  be  obtained  and  plotted  as  shown  in 
I'll,,  f.7  iQ  -.r.c!  ic«r.i.  t!'c  efficiency  of  the  engine  by  tha  oae 

oj.  V  i;/nt,  >v..‘ 7)  of  I,  Sriper lapcsed  t-n  the  basic  grid  of  Fig.  53.  'Tie  kinks 
ill  the  cvu'Vi.’P  i.  t  •’■2.2  have  Ic-re  been  smoothed  out. 


FutI  Spttihc  Imriils*  for  Typical  Enfant 
•A'itt-.  Fixtd  tnl*t  ct.d  Ex't  Gtomotry 


Optimum  engine  thrust  efficiency  corresponds  to  optimum  1^, 
(’*’hich  for  the  particular  case  shown  occurs  at  a  Mach  number  of  2. IS 
and  i  thrust  coefficiert  of  0.535).  The  locus  of  this  optimum  point 
set  by  the  shape  of  the  combustion  efficiency  and  the  diffuser  air 
flow  curves.  As  the  particular  engine  design  was  set  originally  by 
considerations  of  thrust  rather  than  of  efficiency,  the  diffuser  and 
combustor  optimum  design  points  do  not  coincide.  A  higher  Impulse 
peaK  value  could  have  been  obtained  by  arranging  tn^  engine  geoiaactry 
so  that  the  diffuser  operated  at  peak  pressure  recovery  while  combus- 


tier)  took  p.’iace  at  maximum  efficiency.  The  choice,  however,  warn  dic¬ 
tated  by  the  possibility  of  obtaining  greater  thrust,  by  operating  the 
combustor  at  an  equivalence  ratio  corresponding  to  higher  than  optlaua 
combustion  efficiency.  Maximum  thrust  is  seemed  by  uatchlng  the  dif¬ 
fuser  inlet  and  combustor  exit  areas  to  the  maximum  available  heal  re¬ 
lease,  rather  than  by  achieving  the  most  efficient  combustion.  A 
whole  series  of  engine  designs  may  thus  be  assembled  from  a  combustor 
and  a  diffuser  of  givc.n  performance  by  adjusting  inlet  and  exit  areas 
to  optimize  economy  at  various  points.  Figure  58  shows  optimum  fuel 


Fig.  5i'  Fuel  Specific  Impuisc  for  Typicoi  Engine 
with  Vcriobl*  Inlet  ond  Exit  Geometry 


specific  Impulses  associated  with  various  values  of  design  thrusts  for 
engines  having  the  diffuser  and  combustor  perfprmance  shown  in  the 
preceding  figures,  but  with  inlet  and  exit  areas  adjusted  for  optimum 
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efficiency  at  each  point.  This  figure  thus  represents  the  performance 
oi  a  variaole  geometry  engine— If  such  could  be  built.  A  comparison 
of  Figs,  hi  and  58  shows  In  effect  the  penalty  paid  by  a  fixed  geome¬ 
try  engine  for  off-deslgn  perf ornaiice .  The  fuel  specific  Impulses  for 
the  fixed  geometry  engine  approach  those  for  the  variable  geometry  one 
at  the  optimum  point,  but  fall  20  to  30  per  cent  below  them  at  the  ex¬ 
treme  values  of  thrust  coefficient  and  Mach  number. 


5.2  Optimum  Fuel  Specific  Impulse 

To  find  the  optimum  of  the  parameter  of  fuel  specific  Impulse, 
it  is  first  necessary  to  maximize  the  expression 


^net/*f 


by  specifying  the  thrust  value  ,  the  fuel  flow,  Is  thus  re¬ 

garded  as  the  variable.  Since  Is  la  turn  a  product  of  the  engine 
air  flow  and  the  fuol-to-alr  ratio  y, 


I 


f 


The  two  terms  In  the  denominator  must  be  related  to  each  other  through 
the  medium  of  the  thrust  coefficient  equation  and  the  internal  flow  re¬ 
lationships  for  particular  engine  designs  as  characterized  by  their 
geometrical  features.  These  are  implicit  relations  so  that  the  char- 
act'.'rlstlc  design  features  must  be  such  that  values  can  be  arbitrarily 
assigned  to  them  In  order  lo  evaluate  the  required  functions.  The  deri¬ 
vations  given  below  will  show  that  these  features  are  the  diffuser  exit 
Mach  number,  and  the  ratio  of  nozzle  exit  area  tu  combustion  uhamber 

area,  (A„/A„).  These  are  related  to  the  equations  for  thrust  coeffl- 
dent,  diffuser  flow,  combustion  chamber  momentum  balance,  and  exit  noz¬ 
zle  flow,  as  follows; 
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«,/i  A,  A„ 

- IL^-  —±  r  0(Ji)  -  2ji - S -  .  (128) 

(q/P)^  2  AgCq/P)^ 


Ai  [(P/P.)  S]^  P.^  -  Aj  [(P/P.)  8]^  P,^  , 


1  -  (Cp  -  2)(q/P)2 

- V -  ■  <V  "<“>6  ■ 


One  «sy  evaluate  (A^/Ag)  from  Eq.  (129)  in  order  to  express  it 
in  terms  of  functions  of  11^,  and 


*(  't,  f''”*’  *]p 

^2  ■  K,  [(P/P.)  s]^ 


■  ’>d 


(P/Pt)  * 

(P/P^)  ^ 


Solving  for  (S  /  V'T.  )  in  Eq.  (130), 
a  tQ 

1  -  (Cn  -  2)(q/P) 

<vht;-)  -  (jtrij;) - V - ^ 


Entering  this  value  into  Eq.  (128),  along  with  Eq.  (131),  there 


follows 


Cy(q/P)o 


«o  ’Id  b 


(P/Pt)  «]^ 


»2  -I 
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-  2(q/P)^ 


i(P/V. )  S 

i-  *■  -^2 

[(pTiTsT 


M- 


This  may  be  solved  for  0(M)g/0(M)jj  to  obtain 


0(M)b 


(P/P^) 


’’e  <P/Pt>2  ['  -  -  2)(q/P)J 

[(P/P^)  Sj 


(q/p>o 


2  ’Jd  1 - 

°  [(p/p^)  m  ^  0 

L  ^  J  0 


(i;i2) 


since  fx,  Cp  ,  Cp,  Mq,  Vl^,  and  have  been  fixed,  the  right 
b 

hand  side  of  Eq .  (132)  may  be  evaluated  by  defining  the  particular  re¬ 
lationship  between  0(U)g  o-nd  that  must  be  met  if  the  engine  Is 

to  produce  the  required  Cp  at  these  conditions.  This  relationship 
implicitly  determines  the  throat  area  ratio  (Aj/A^.)  because  particular 
values  for  0(M)  and  0(M).  correspond  to  each  value  of  (A  /A  )  and 
(Ag/A^).  (Either  one  ratio  or  the  other  will  have  been  fixed  so  that 
they  both  change  together  as  A^  is  varied.)  As  A^  is  reduced  below  A^ , 

the  exit  stream  becomes  supersonic,  and  0(U)  also  increases  with  11., 

©  © 

but  Its  rate  of  Increase  is  slow.  The  subsonic  value  of  0(H)|,  like¬ 
wise  Increases  with  decreasing  but  does  so  more  rapidly;  the  ratio 

of  0(M)  /0(1I).  is  therefore  unity  for  A.  -  k  -  A„,  and  decreases  as 

6  D  ©  C  <w 

A^  diminishes.  Figure  59  shows  this  ratio  pl-oLted  against  the  combus- 
tio.i  chamber  to  exit  nozzle  area  ratio,  as  obtained  from  Eq.  (132). 

From  the  abscissa  of  this  figure  may  be  determined  the  exit  nozzle  con¬ 
traction  ratio  for  the  chosen  value  of  (A^/Ag).  The  value  of  which 
corresponds  to  this  value  of  (A/A  )  may  then  be  found  from  Mach  number 
tables  in  order  lu  evaluate  0(M)^,  which  is  then  inserted  into  Eq.  (130) 

to  determine  the  value  of  (S^/  ^T^  )  associated  with  the  operating 
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value  of  Mg  to  give  the  prescribed  Cj..  A  set  of  such  calculations 
must  be  performed  for  various  values  of  Mg  to  obtain  corresponding 
valuer,  of  the  temperature  parameter.  for  each  engine  design 

are  computed  from  the  respective  values  of  M2  by  means  of  Eq.  (129), 
and  i:he  corresponding  fuel-lo-alr  ratios  are  deduced  from  the  values 
of  {ii  /  JY~) ,  as  obtained  above.  The  value  of  the  fuel  specific  im- 

a  f  ^2 

puis.=‘  is  then  computed  for-  eaf  "  design,  and  the  values  plotted  agrirst 

diffuser  Mach  number,  inlet  area  capture  ratio,  equivalence  ratio  or 

exit  nozzle  throat  area  ratio,  in  order  to  find  the  relation  of  any 

one  of  these  to  the  optimum  fuel  specific  impulse.  Other  values  of 

the  engine  design  parameters  may  also  be  determined  from  Eqs. (128) 

through  (130)  Such  an  optimum  value  of  if  valid  only  for  the 

selected  value  of  (A  /A„);  if,  therefore,  (A  /A„>  is  also  a  variable 
0  «  0  2 

for  optimization  purposes,  the  process  needs  to  be  repeated  for  a 
range  of  values  to  select  a  further  optimum  as  a  function  of  this  ratic. 
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Piectuliif;  set  tiuns  l-.avo  troalo<i  '.lie  general  operation  of  engine 
compoiu'iUs  tin  <,■•)  aa  t ;  11.0  flew  state'  across  each  com- 

ooaoni.  !i  is  t  lijiv  »>.’ f  p^ssi^lc,  in  principle  at  leact,  to  combine 
various  dtllusiTs,  cumhu^lui  s ,  and  exit  'lozzics  to  produce  t.  wide 
vuM'  l;,  ;  wlm  h  may  bo  operated  under  various  conditions, 

siiive  the  i:as  d;.  .i.iniio  aspoels  of  engine  design  arc  [tore  amenable  to 
ni'.nlvsis  tha-i  are  tlu'  Structural  aspects.  Although  all  design  is  in  a 
sense  ipl 1 raiza t ion .  the  procedure'  discussed  here  Involves  a  comparison 
iictwcoe  vario'is  compone,‘ii‘.  designs,  each  optimized  with  respect  to  the 
demands  upon  it,  in  oiuci'  to  .select  the  ucst  compj'omise  lor  the  engine 
as  a  whole'. 


The  t!\ree  ba.sic  vari.ible.s  affecting  the  output  of  a  heat  engine 
art  an  intensive  factor,  the  energy  added  per  unit  of  vorklng  fluid, 
an  extensive  factor,  the  quantity  of  working  fluid  handled, and  the 
effii. iency  of  the  cycle.  Once  the  choice  of  a  ramlet  cycle  hus  been 
made,  various  combinations  of  these  factors  are  possible.  The  practical 
variables  cerresponding  to  these  are  the  temperature  parameter  (S^/yx^) 
and  the  diffuser  area  ratio  (A„/A^).  There  are  other  parameters  which 
also  need  to  be  considered  when  a  design  is  being  calculated,  but  as 
these  are  either  environmental  (such  as  flight  Mach  number,  pressure, 
and  temperature)  or  experimental  (such  as  burner  drag),  they  are  ir.l- 
llaily  determined  constants  rather  than  dependent  variables.  The  dis¬ 
charge  area  of  the  exit  nozzle  occupies  a  peculiar  place  in  the  calcu¬ 
lation,  as  does  diffuser  cowl  drag,  in  that  except  for  the  considera¬ 
tion  of  external  aerodynamic  drag,  It  Is  at  an  obvious  maximum  when 
the  gases  have  expanded  t.o  ambient  pressure.  If  the  external  drag  is 
not  to  he  considered  in  the  optimization  of  an  engine,  all  designs 
will  i.'ivolve  an  expansion  to  ambient  pressure.  On  the  other  hand, 
si.nce  the  rale  of  change  of  thrust  with  area  is  zero  when  the  ambient 
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(  ovi’oust  i!>  i  •  rti.ioiK'v  1-  t»io  c'j'jrsc  of  opt  imlz.it  ion  procedure 
Hii‘'  * !-  jmr.iif;  I  >ir  ilesigncd  foi  optiniur.  combustion  efficiency 

:  tb'-  ^ucl  to  :iir  r.'.rio  Irdicated  by  the  analysis.)  The  likely  range 
of  diffuser  exit  Mach  numbers  aill  lange  from  0.250  down  to  0.100  for 
'  eni'ln*'  Jesigns  .lepomling  on  the  tenperr.ture  ratio  across  the 

eoiuiMisior .  Tne  pleasure  ratio  across  the  combustor  will  also  be  such 
as  to  require  e-xpansior.  of  the  gases  to  suhstaiitially  the  full  engine 
area  to  bring  them  close  to  ambient  pressure  at  exit.  It  will  he  con¬ 
venient  io  consider  six  values  of  diffuser  exit  Mach  number  from  0.100 
to  0.250,  at  intervals  of  d.025,  and  to  evaluate  the  exit  nozzle  throat 
area  and  the  temperature  parameter  (S^/  from  Eq.  (132)  for  the  re¬ 

quired  value  of  thrust.  Combustor  drag  values  may  he  chosen  to  he  2, 

6.  and  12  combustion  cnamher  inlet  velocity  heads.  These  values  corre¬ 
spond  respectively  to  the  drag  range  associated  with  either  a  huffle- 
lype  burner  at  the  diffuser  exit,  a  can-type  combustor  placed  wholly 
in  the  combustion  chamber,  or  with  the  latter  mounted  wlthun  the  dif¬ 
fuser  near  its  exit. 

Six  sets  of  calculations  are  then  performed  using  the  equation 
referred  to,  introducing  in  sequence  the  six  values  of  II2 ;  values  of 
(,K^/K^)  and  for  orag  coefficients  of  2,  6,  and  12  are  then  deter¬ 
mined.  Since  it  is  assumed  that  the  exhaust  gases  are  to  bo  reexpanded 
to  full  engi.ne  area,  may  he  arct  equal  to  Ag.  The  exit  forces  are 
next  computed,  except  for  the  term  which  is  determined  when  the 
oicit  is  matched  with  an  inlet  corresponding  to  the  assumed  values  of 
diffuser  exit  Mach  number  M2  from  Eq .  (129)  for  each  of  the  values 
0.7,  0.8,  and  0.9  for  tj^.  There  are  thus  nine  famlMes  of  six  seta  of 
engine  designs  requii  ing  various  values  for  the  teii.pt  .-ature  parameter 
(S^/'  yr^)  to  be  determined  fi-om  values  of  thrust  coefficient,  diffuser 
efii'  ieiicy,  diffuser  exit  Mach  number, and  combustor  drag.  One  may  now 
select  hy  interpolation  from  each  sot  the  engine  design  yielding  the 
desired  thrust  coefficient,  and  calculate  the  fuel  flow  and  fuel  spe¬ 
cific  impulse  for  the  selected  thrust  coef f iclei.t  aa  explained  in 
Section  2,  page  31,  and  Section  5,  page  132,  and  then  plot  these  re¬ 
sults  to  determine  the  optimum  fuel  specific  impula’o  for  each  family. 
The  resultant  fuel  specific  Impulse  curves  are  shown  in  Figs.  60 


! 


M;r< 'Uiii!  bi:.  ait  '..omputod  100  per  cent  combustion 

Iti  u  tn  >  iHil  vv.*\  t.isi  ■  b’  :^!!;uslt‘d  to  Itjwcr  values  by  nultipiica- 
t  ioi\  '  f  t.'.e  trojaati  aiui  divlsio!:  ui  ine  abscissa  by  any  appropriate 
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Fij.  62  Efftct  of  Burntr  Drag  and  Equivaltnct  Ratio  on  FutI 
Sptcific  Impulst,  rij  ^  0.7 


The  design  of  an  engine  for  optimum  fuel  specific  Impulse  with 
expansion  of  exit  gases  to  ambient  pressure  Instead  of  to  the  fixed 
nominal  engine  area  wou.d  differ  only  In  regard  to  the  selection  of 
the  proper  value  of  the  area  ratio  (A^/Ag)  which  fixes  the  exit  Mach 
number.  Since  ^  Is  needed  to  obtain  the  ratio  of  the  throat  to  com¬ 
bustor  area,  an  approximate  value  must  first  be  employed  to  obtain  a 
valve  for  the  latter,  since  changes  of  ^  affect  the  design  less  than 
changes  in  (A^/Ag).  One  may  now  proceed  to  determine  other  features 
of  the  engine  design  that  depend  on  th>  rlenlred  thrust.  For  example, 
the  exit  Mach  nujnbei-  may  he  obtained  irom 

(A/A*)^  =  M  . 


(133) 


anu  tljc  o>.it  pi'os.b‘ir('  froa 


(A-/A  ) 

•— i— 2-(s,  /  )  w 

n  “  '  '-o  0 


0(M)^(P/P^) 


Tho  quantities  (A^  'A. ).  and  ( ''2^ /"^t  ^  w:j''n  the  opti¬ 

mum  desl^jn  for  selected  valueti  of  Cp,  ,  and  ^  w.is  beini?  established. 
The  value  of  is  found  from  (A/A  bjf  means  01  £q,  (J3i»),  anu  it  is 
.,t,  a  < (^2  tlie  functions  of  in  Eq  (134).  The  value  of 
(P^  P,j)  is  not  likely  to  turn  out  to  be  exactly  1.0,  corresponding  to 
the  case  P^ :  if  ^  P„,  the  further  expansion  required  points  to 

the  need  to  select  a  iarg-er  value  of  p  ;  likewise,  if  P^  <  P^^,  a  smaller 
v.iliu'  (li  j,  is  required.  A  second  more  accurate  trial  value  of  4, 

fi' ,  may  then  be  derived  by  computing  the  exit  total  pressure  P^  ,  using 

'"e 

values  found  for  ,  0(M)^,  and  {K^/k^)\ 


<'‘l"‘2)<Vt"qh  »<“>e 


Having  found  P^  from  this  equation,  one  may  next  c:omputfe  (Pq/P.  )  to 
e  ,  e 

yield  the  exit  Macn  number  M  for  nozzle  expansion  down  to  atmospheric 

I  ®  341  ' 

pressure.  This  value  of  M  also  yields  (A/A  )  ,  the  new  value  for 

6  I  0 

(A^/A^),  and  the  corresponding  value  of  4  is  then  derived  from  this 
by  multiplying  by  that 

4'  -  (A/A*)g(A^/A2)  .  (136) 

I 

'.lese  calculations  are  then  repeated,  using  4  to  obtain  a  second  ap¬ 
proximation  for  (A^/Ag) ,  (Aj/Ag),  etc.  The  calculations  implicit  in 
Eqs .  (133)  and  (134)  are  lepeated  to  determine  whether  the  exit  static 
pressure  then  matches  the  ambient  value. 

A  question  arises  concerning  the  basis  tor  a  first  estimate  for 
4,  as  it  varies  widely,  ranging  in  value  from  a  quantity  of  the  order 
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of  unity  for  high  thrust  ooof f icient  designs  down  to  much  tower  values 
fui'  low  thrust  designs.  The  wisdom  ot  using  exit  areas  larger  than 
the  combustion  chamber  is  questionable  since  the  gain  in  thrust  on  the 
downstream  side  of  the  nozzle  is  offset  by  the  Increased  frontal  drag 
produced  by  increacing  the  engine  diameter  from  A„  to  A  .  Accordingly, 
the  assumption  of  ii  =  1  is  a  reasonable  first  approximation. 

Determination  of  an  engine  design  to  provide  maximum  thrust  is 
in  practice  more  straightforward  than  that  shown  in  the  preceding  cal¬ 
culations  since  only  the  maximum  value  of  (S  '  i/tT)  need  be  considered. 

In  cases  of  low  burner  drag,  the  engine  air  flow  will  be  at  a  maximum, 
corresponding  to  sonic  discharge  from  the  combustion  chamber  without 
an  exit  nozzle.  The  calculation  proceeds  from  the  momentum  balance 
across  the  combustion  chamber,  starting  with  -  1 ,  and  (S^^/  )  at 

a  maximum;  inserting  these  values  into  Eq .  (135)  enables  one  to  deter¬ 
mine  Mg.  The  inlet  area  ratio  is  then  found  from  M^,  the  diffuser 
total  pressure  efficiency,  and  the  weight  flow  relationship,  Eq.  (129). 

Exit  and  inlet  fcrce  coefficients  may  now  be  computed,  and  their  dif¬ 
ference  is  then  the  desired  thrust  coefficient.  If  the  combustor  drag 
is  high,  it  will  be  found  that  substantially  as  much  thrust  will  be 
obtained  using  a  slightly  convergent-divergent  exit  nozzle  “ 

0.94, say)  as  with  an  unconstricted  combustion  chamber  —  in  addition  to 
better  fuel  economy. 

6.2  Typical  Results  of  Engine  Optimization  Calculations 

The  procedure  outlined  in  the  preceding  section  has  been  carried 
out  for  some  engine  design  cases  of  general  interest.  The  approach 
chosen  has  been  to  invertlgate  first  the  dependence  of  Mach  number  on 
the  engine  fuel  specific  impulse.  This  has  then  provided  a  basis  for 
the  selection  of  the  design  flight  speed  when  it  was  combined  with  ex¬ 
perimental  test  data  on  the  parameter  of  diffuser  efficiency  as  a  func¬ 
tion  of  flight  Mach  number.  Other  engine  performance  parameters  such 
as  thrust  coefficient,  burner  drag,  combustion  efficiency,  and  exit 
nozzle  efficiency  may  be  considered  fixed  since  they  are  not  dependent 
on  the  flight  Mach  number.  The  next  stage  of  optimization  is  then 
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earned  out  at  the  .selected  flight  speed  and  the  corresponding  availa- 
tilc  diffus_-r  pressure  lecovery  in  order  to  determine  the  best  com¬ 
promise  hc'tween  combu.stion  efficiency,  burner  drag,  and  exit  nozzle 
ef  f  ic  tenc^y .  Tlie  study  oi  the  compromise  between  burner  drag  and  exit 
nozzle  offic  tenrv  «ers  c  onioii  St  ion  efficiency  may  be  carried  out  most 

simply  if  one  notes  that  the  eifect  of  combustion  efficiency  is  to 
alter  the  amount  of  fuel  to  be  burned — in  an  inverse  relationship. 
Accordingly  the  values  of  fuel  specific  impulse  calculated  for  »  com- 
bustio'’  efficiency  of  100  per  cent  may  be  adjusted  to  any  lower  com- 
h.ustion  ciliciency  by  increasing  the  fuel  flow  inversely  with  the  de- 
Ci-ease  in  combustion  efficiency.  This  leads  to  a  new  value  for  the 
fuel  specific  impulse  equal  to  the  product  of  the  fuel  specific  im¬ 
pulse  for  100  per  cent  combustion  efficiency  and  the  actual  combustion 
efficiency,  and  to  a  new  optimum  equivalence  ratio  equal  to  the  quo¬ 
tient  of  the  ideal  equivalence  ratio  Into  the  combustion  efficiency. 

A  relalionship  between  engine  geometry  and  the  fuel  specific  im¬ 
pulse  is  shown  in  Figs.  63,  64,  and  65  for  flight  Mach  numbers  of  1.50, 
2.25,  and  3.50.  (Figure  6i  is  here  repeated  for  convenient  reference 
as  Fig,  64.)  Equivalence  ratio  is  plotted  as  the  abscissa,  and  there 
Is  a  simultaneous  change  of  inlet  area  ratio  (A, /A  )  and  exit  nozzle 
throat  size  to  maintain  a  constant  thrust  coef  f  li’ >wr.t .  The  maximum  oi 


i 

Fig.  63  Effect  of  Burner  Drug  and  Equivalence  Ratio  on  Fuel 
Specific  Impulse,  .  0.9 
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this  cuivc  do t eruii nc-s  sin.jltaneously  tho  proper  inlet  and  exit  areas 
as  well  as  the  equivalence  ratio  for  optimum  performance.  A  S(.>t  of 
such  optimization  calculations  has  been  made  using  suitably  selected 
values  of  diffuser  recovery,  and  resulting  values  of  optimum  fuel  spe¬ 
cific  impulse  .are  plotted  ago  •  -I  iiio  rii^IiL  nacii  number  as  the  abscissa 
in  Fig.  eb.  This  figure  is  not  complete  without  the  insertion  of  a  eon- 
sicieraiioii  of  ihe  recoveries  actually  achievable  in  practical  diffusers. 


Fig.  66  Effect  of  Flight  Mach  Number  and  Diffuser  Efficiency 
on  Fuel  Specific  im.ptibe 


I 

i 

! 

J 
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Reference  is  made  to  Fig.  67  whlcli  shows  conservative  values  for  dif¬ 
fuser  efficiency;  these  are  incorporated  in  Fig.  66  to  establish  the 
proper  relationship  between  the  optimum  fuel  specific  impulse  and  the 
flight  Mach  number.  The  fuel  specific  impulses  shown  in  Fi^.  S6  were 
computed  on  the  basis  of  an  assumed  exit  nozzle  efficiency  of  100  per 
cent,  a  burner  drag  equal  to  6  velocity  heads,  and  100  per  cent  com¬ 
bustion  efficiency.  An  assumption  of  100  per  cent  nozzle  efficiency 
is,  of  course,  somewhat  optimistic:  the  burner  drag  figure  Is,  however, 
a  realistic  one.  (The  combustion  efficiency  has  no  bearing  on  the 
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I'ig.  67  Expc-rimentol  DiHuser  Reccvery  Ratios 


locution  of  the  optimum  as  long  as  it  can  be  ccnsldorod  coiiotaut . )  On 
these  assumptions  the  flight  Mach  number  for  optimum  fuel  specific  im¬ 
pulse  will  he  found  to  occur  at  a  value  of  about  2.6.  An  analysis 
similar  to  the  above  may  be  performed  using  the  more  pessimistic  value 
for  the  exit  nozzle  efficiency  tj^  oi  U.67  in  order  to  obtain  an  in¬ 
crease  of  ex:.t  nozzle  ;i.rust  above  that  available  from  a  sonic  nozzle; 
the  results  of  this  calculation  are  presented  in  Fig.  66.  A  comparison 
of  Figs.  66  and  66  shows  that  a  loss  of  exit  nozzle  efficiency  has  a 
large  detrimental  effect  on  the  magnitude  of  the  available  fuel  spe¬ 
cific  impulse  for  a  given  diffuser  efficiency  and  flight  Mach  number, 
but  it  does  not  alter  appreciably  the  flight  Mach  number  at  which  the 
optimum  Impulse  is  obtained. 

The  diffuser  efficiency  strongly  affects  both  the  optimum  geome¬ 
try  and  the  operating  equivalence  ratio,  as  shown  in  Figs.  69  and  70, 

As  the  diffuser  efficiency  rails,  the  fuel  specific  impulse  is  opti¬ 
mized  by  reducing  the  diffuser  air  flow,  as  shown  in  Fig.  69,  at  the 
expense  of  greater  combustor  heat  addition  (Fig.  70)  and  increased 
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F'n.  70  Voriatian  of  Optimv>n  Equivoltnct  Rotia  >vith  Frtottrtom 
Mach  Numbor  and  Dilfuttr  Efficiency 


combustion  losses.  The  effect  of  on  fuel  specific  impulse  has  hnen 

seen  In  Figs.  66  and  68.  The  burner  drag  losses  have  an  effect  similar 
to  that  of  diffuser  pressure  losses  except  that  they  are  related  to  the 
local  combustion  chamber  Inlet  Hach  number  Instead  of  to  the  freestream 
Mach  number.  Figure  71  shows  the  reduction  of  optimum  fuel  specific  Im¬ 
pulse  with  varying  burner  drag  coefficient.  The  effect  will  be  observed 
to  be  large  at  low  fllgut  Mach  numbers  (1.5  to  2.0)  where  the  diffuser 
exit  velocity  bead  Is  a  significant  part  of  the  total  piessurc  loss, 
but  at  higher  flight  Mach  numbers  from  2.25  to  3. 40, the  diffuser  total 
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Fig-  71  Effect  of  Burner  Drog  on  Fuel  Specific  Impulse 


pressure  losses  Inr'rense  considerably  and  so  swamp  the  effect  of  burner 
drat;.  As  the  flight  Mach  number  increases,  the  burner  Mach  number  di¬ 
minishes  because  of  the  greater  density  of  the  captured  air,  and  inter¬ 
nal  drag  losses  are  thereby  reduced.  Because  of  the  relationship  of 
interiial  drag  losses  to  the  diffuser  exit  Mach  number,  there  is  a  marked 
shift  of  the  optimum  diffuser  Mach  number  with  combustor  drag  coeffi- 
rlont,  as  shown  in  Fig.  72.  Tit-  ovtr-'ll  efficiency  in  burners  of  Uig'a 
drag  is  improved  if  a  lower  diffuser  Mach  number  is  employed,  i.e.,  the 
distribution  of  cor.busto,  aerodynamic  and  thermodynamic  losses  is  read¬ 
justed  to  minimize  the  sum  of  these  quantities. 

The  diffuser  exit  Mach  nu  '  r  for  optimum  fuel  specific  Impulse 
is  relatively  insensitive  to  chai  in  diffuser  efficiency,  exit  noz¬ 
zle  efficiency,  and  thrust  coefficient.  This  inseusltivity  arises 
from  a  self-adjusting  feature  associated  with  the  diffuser  and  combus¬ 
tor  geometry.  A  reduction  of  diffuser  efficiency,  for  Instance,  calls 
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Fig.  72  Voiiatien  of  Diffuttt  Meek  Number 
wirk  Buner  DfOg  for  Optimum 
Fuel  Specific  Impulse 


for  a  lessened  capture  area  ratio  which,  through  the  reaultant  Increase 
In  the  ratio  of  diffuser  exit  to  inlet  area,  offsets  the  lower  pressure 
recovery  and  brings  the  air  essentially  to  a  constant  Mach  nus'jer  In 
the  dlffuani-.  hlallarly,  variations  In  exit  nozzle  efficiency  and  de¬ 
sign  thrust  coefficient  principally  affect  the  design  of  the  exit  noz¬ 
zle  and  combustor,  together  with  the  value  of  the  over-all  specific 
Impulse.  This  variation  ol  the  lapulse  with  exit  nozzle  efficiency  Is 
shown  for  over-all  nozzle  efficiencies  of  0.975  and  0.950  in  Fig.  73 
for  engines  of  optimized  georietry.  The  ratio  of  the  impulse  loss  to 
the  loss  of  exit  momentum  Is  nearly  constant  with  flight  Hach  number, 
but  the  proportion  of  the  former  loss  becomes  gr**ter  at  the  lower  spe¬ 
cific  Impulses  obtained  at  lower  flight  speeds.  The  fuel  flow  and  the 
equivalence  ratio  Increase  inversely  with  the  change  in  fuel  specific 
Impulse.  £xlt  nozzle  efficiency,  as  explained  earlier,  is  a  function 
of  the  divergence  angle  and  the  exit  to  throat  area  ratio.  Subsonic 
and  sonic  nozzles  achieve  nearly  theoretical  values  of  exit  momentum 
because  of  their  relative  uniformity  of  flow  and  ti.o  absence  of  flow 
divergence.  Since  a  supersonic  nozzle  of  fixed  length  will  approach 
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Fig.  73  Effect  of  E«H  Noiile  Efficiency  on  Optimum  Fuel 
Specific  Impulse 


the  ideal  as  the  area  ratio  approaches  unity,  engine  design  uptlBlza- 
tion  calculations  Involving  such  exit  nozzles  should  enploy  a  nozzle 
efficiency  varying  with  the  expansion  ratio.  A  fair  approximation  to 
this  situation  was  obtained  in  the  calculations  for  Flga.  66  and  68 
by  assiimtnj  than  th**  tncreas-'  of  exit  soscr.tus  above  the  sonic  value 
for  the  exit  stream  was  a  constant  fraction  of  the  Ideal.  Since  the 
exit  momentum  Is  ^a*a®^**^eff  sonic  value  Is  the  approxl- 

natlon  takes  the  foi-m 

0(II)^ff  -  1  +  |0{M)ideal  -  l]  ^e  •  <137) 

An  increase  in  the  value  of  the  design  tnrnst  coefficient  leads 
i.,j  aii  iiiciease  in  the  combustion  temperature  and  exit  nozzle  area,  and 
a  fall  i.i  the  over-all  fuel  specific  impulse  wi'.h  oniy  a  slight  effect 
on  the  optimum  diffuser  exit  Mach  number  which  therefore  increases 
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propc  rt  idii.i  n y  at  iiie  rale  ol  about  only  one-te.ith  of  that  of  the 
thrust  i.  han^f  The  i  'tcct  of  incroaslnR  the  thrust  coefficient  on 
fuel  specific  impulse  and  equivalence  ratio  for  three  different  fl.if,^' 
Math  numbers.  1.8,  2 .  “i ,  and  S-.l,  is  shown  in  FJg.  74;  these  curves 
cover  likely  roqions  of  liirust  and  flight  speed  operation  with  likely 
values  of  diffuser  recovery,  burner  drag,  and  exit  nor.zle  efficiency. 
Similar  effects  occur  at  other  flight  speeds. 


Fig.  74  Effect  of  Thrust  Cocfficitnt  on  Fuel  Specific  ir.ipulse  and 
Optimum  Equivoltnce  Rotio 


In  conclusion,  It  may  said  that  tne  most  straightforward  way 
of  arriving  at  an  ontimi.in  engine  design  configuration  iu  first  to  assign 
prohable  values  to  the  various  component  performance  parameters  for  use 
in  the  indicated  calculations.  The  latter  are  then  repeated  using  upper 
and  lower  llsit  values  in  the  case  of  components  of  uncertain  perform¬ 
ance  in  order  to  determine  the  variation  of  optimum  pi'uportlons  and 
performance  within  this  possible  component  variation.  By  this  means, 
one  is  able  to  weigh  the  net  performance  gains  to  be  made  oy  optiBlza- 
tior,  of  one  component  at  the  expense  of  another . 
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Tabulation  of  Mach  Number  Functions 

The  flow  relationships  involving  Mach  number  also  Involve  the 
value  of  y  (which  decrea"  u  with  temperature  fo  any  given  gas,  at 
least  up  to  5000°F,  at  aii,  pressures  except  the  lowest),  but  a  change 
in  y  from,  say,  7/5  to  9/7  produces  a  change  of  only  a  few  per  cent  in 
most  Mach  number  functions  at  any  particular  value  of  the  Mach  number. 
(The  area  function  (A/A*)  Is  particularly  sensitive  to  a  varlatioii  in 
y  at  the  higher  values  of  M  however.)  The  Mach  number  functions  tabu¬ 
lated  below  are  presented  for  comparison  at  these  two  values,  since 

2y 

they  give  simple  whole  number  values  for  such  exponents  as  - — ■ — ^  , 


_ 2_  y  +  1 

y  -  T’  ^'(y  -  1 


,  etc.,  as  Indicated  In  Section  2;  these  values  are 


respectively  9,  7,  and  4  for  y  ■  9/7.  The  value  7/5  corresponds  to 
air  at  low  temperatures  (0-500®F),  while  y  -  9/7  corresponds  to  the 
products  of  combustion  of  fuel  aud  air  mixtures  at  more  elevated  tem¬ 
peratures  (about  -WOO'F).  The  set  of  Mach  number  functions  for  V  -  7/5 
may  thus  be  used  to  describe  flow  into  a  combustion  system,  while  those 
for y =  9/7  describe  flow  out  of  a  combustion  system.  For  low  heat  re¬ 
leases,  howe  '-'r,  it  may  be  more  appropriate  to  use  the  7/5  value  for 
both  such  streams. 


The  functions  listed  below  are  Intended  to  be  uniform  with  Refs. 
22  and  23  and  the  present  text.  Other  tabulations  may  be  found  In  Ref. 
i*!  and  in  some  texts  listed  la  the  Bibliography.  (Quantities  are  all 
dimensionless  unless  otherwise  stated.  Figures  A-l  through  A-17  illus- 
trp.te  the  variation  with  Mach  number  of  the  various  functions  for 
values  of  the  latter  up  to  5.  (It  has  not  seemed  worthwhile  to  extend 
the  rei at innshipv  beyond  this  value  owing  to  the  need  to  incorporate 
real  gas  effects  at  higher  Mach  numbers.)  To  avoid  confusion  with  the 
curve  of  (Tj^/T^;,  the  function  (p*/p^)  has  been  plotted  separately. 
Constant  specific  heat  values  are  assumea  in  the  derivations. 
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